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Abstract 



We calculate the leading power corrections to the decay rates, distributions and 
hadronic spectral moments in rare inclusive B Xgi^l^ decays in the standard 
model, using heavy quark expansion (HQE) in (l/mb) and a phenomenological model 
implementing the Fermi motion effects of the 5-quark bound in the -B-hadron. We 
include next-to-leading order perturbative QCD corrections and work out the depen- 
dences of the spectra, decay rates and hadronic moments on the model parameters in 
either HQE and the Fermi motion model. In the latter, we take into account long- 
distance effects via B ^ Xg + {J/tp, tp', ...) — > Xs£'^£~ with a vector meson domi- 
nance ansatz and study the influence of kinematical cuts in the dilepton and hadronic 
invariant masses on branching ratios, hadron spectra and hadronic moments. 
We present leading logarithmic QCD corrections to the b — > S77 amplitude. The 
QCD perturbative improved Bg — > 77 branching ratio is given in the standard 
model including our estimate of long-distance effects via Bs ^ <fry ^ 77 and 
Bg — > (pip ^ (/)7 ^ 77 decays. The uncertainties due to the renormalization scale 
and the parameters of the HQE inspired bound state model are worked out. 

Zusammenfassung 

Wir berechnen im Standardmodell mit der Entwicklung schwerer Quarks (ESQ) Kor- 
rekturen in (l/mb) zu Zerfallsbreiten, Verteilungen und hadronischen Momenten 
des seltenen, inklusiven Zerfalls B — > . Dieselbe Analyse haben wir 

im phanomenologischen Fermibewegungsmodell durchgeftihrt, welches die Effekte 
des im B-Meson gebundenen 6-Quarks beschreibt. In beiden Methoden zeigen wir 
unter Einbeziehung nicht fuhrender QCD-Korrekturen die Abhangigkeit der Spek- 
tren, Breiten und Momente von den ModeUparametem. Mit einem Vektormeson- 
dominanzansatz modeUieren wir im FermibewegungsmodeU die langreichweitigen 
Beitrage der Charmoniumresonanzen, die durch B ^ Xs + {J/i>, tp', •••) — > Xsi~^i~ 
entstehen. Wir studieren den EinfluB kinematischer Schnitte in der invarianten Masse 
des Leptonpaares und des hadronischen Endzustandes auf Verzweigungsverhaltnisse, 
hadronische Verteilungen und Momente des Zerfalls B Xgf.^l^ ■ 
Ferner prasentieren wir fiihrende logarithmische QCD-Korrekturen in der b 577 
Amplitude. Wir geben das QCD-verbesserte Verzweigungsverhaltnis von Bg — > 77 
im Standardmodell unter Beriicksichtigung unserer Abschatzung langreichweitiger 
Beitrage durch die Zerfallsketten Bg ^ (fn ^ 77 und Bg ^ (j)tp ^ (pj 77 
an. Wir schatzen die Unsicherheiten durch die Renormierungsskala und des Bg 
Bindungszustandsmodelles, welches der Theorie schwerer Quarks entlehnt ist, ab. 
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Chapter 1 



Introduction 



Raie B decays probe the flavour sector of the standard model (SM) ^ and perturbative and non-perturbative 
aspects of QCD. Since the first measurement of rare radiative B K*j decays by the CLEO collabo- 
ration [§] this area of particle physics has received much experimental |^] and theoretical |Q] attention. 
Flavour changing neutral current (FCNC) induced decays are of particular interest, since in the SM they are 
governed by loop effects and depend sensitively on virtual particles like, e.g., the top quark and Cabibbo 
Kobayashi Maskawa (CKM) matrix elements [Q]. 

The theory of the effective Hamiltonian ||6|-|l0|] (see section 22 for a discussion) Tleff ~ J^CiOi 
enables a description of low energy weak processes in terms of short-distance (Wilson) coefficients Cj, 
which can be calculated perturbatively. The (new) vertices 0,, which are absent in the full Lagrangian, are 
obtained by integrating out the heavy particles (W, t, (j) in the SM ) from the full theory. Their coupling 
strength is given by the Ci, which characterize the short-distance dynamics of the underlying theory. 

The Wilson coefficient Cf^^ corresponds to the effective 657 vertex. Its modulus \C^^^\ is constrained 
by the measurement of the inclusive B Xg^ branching ratio at CLEO [11]. The CLEO result is in 
agreement with the present theoretical SM prediction in B Xg^ decay and moreover, can exclude large 
parameter spaces of non-standard models. 

The transition b —>■ sl'^l" with i = e, fi,T involves besides the electromagnetic penguin b — > 57* 
s£~^i~ also electroweak penguins b sZ^* —^ si^£^ and boxes. They give rise to two additional Wil- 
son coefficients in the semileptonic decay B Xgi'^i^ , Cg and Ciq. A model independent fit of the 
short-distance coefficients Cg, Cio and Cy^ is possible from the following three observables [12]: The 
B — > Xs7 branching ratio B{B Xsj), the (partly integrated) invariant dilepton mass spectrum and the 
Forward-Backward (FB) asymmetry [13] in i? — > Xgl^t^ decays. They involve independent combina- 
tions of the Wilson coefficients, which allows the determination of sign and magnitude of C7 , Cg and 
Cio from data. 

The presence of charmonium resonances in the decay B —>■ Xgl^l" complicates this analysis. The 
cc states appear via B ^ Xs + {J/i^, i^' , •••) — > Xgi'^i'^ and can be taken into account in a phenomeno- 
logical vector meson dominance (VMD) ansatz, which is assumed to hold near the J/ipjip', . . . peaks 
[13]. Hence, kinematical cuts on the dilepton mass have to be imposed to remove the dominant res- 
onance contributions and to disentangle the short-distance information from the long-distant one, e.g.. 
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< rn?ji^ — Here 6 is an experimental cut-off parameter and typically of order a few hundred 
MeV. The restriction to certain regions of phase space explains the use of partly integrated spectra. The 
net effect of the resonances is an additional error on the distributions and hinders the determination of 
the short-distance coefficients Cj from experimental data. Hence, they must be evaluated by taking into 



account theoretical dispersion by using different phenomenological models [ 14 ] or further experimental 
input. 

The bound state nature of heavy hadrons can be explored within the S-system. Non-perturbative 
power corrections in {^Qco/mh)^ are systematically obtained with the heavy quark expansion (HQE) 



technique [15], parametrizing distributions, decay rates etc. in terms of higher order matrix elements. 



denoted by Ai, A2 for n = 2. We recall the HQE relation niB = mi, + A — (Ai + 3A2) /2mi, + 0(l/m^) 
between the physical i?-meson and the 6-quark mass, where A accounts for the "binding energy". The 
HQE method has been applied to semileptonic charged current B —>■ X£i>£ and FCNC B — > Xg^ 



[17] decays, and it is known that inclusive spectra are not entirely calculable with the HQE approach [16], 



[18-21]. S-meson wave function effects have been estimated in i? ^ c^^'^ [22,23] and B Xsj 



[^,Q decays with a phenomenological Fermi motion model (FM) [22]. 

We present here a detailed analysis of inclusive B — > Xs£~^i^ decays in the SM with £ = e, fi, (since 
we neglected lepton masses in our calculation the result cannot be used in the r case), following similar 
studies for charged current B —^ Xiv^ and radiative B —>■ Xs^ decays. 



This thesis contains the following points [18,25,26,14]: 



We calculate the I /ml power corrections using HQE techniques in B —^ Xgi^i decays in the 
dilepton invariant mass distribution. This corrects an earlier calculation [111] and has been confirmed 



recently ||27|]. We find that the HQE breaks down near the high q end-point, hence the spectrum 



cannot be used in this region. 

Alternatively, we study B-meson wave function effects with a Gaussian FM. We present the dilepton 
invariant mass distribution and the FB asymmetry and investigate the dependence on the FM model 
parameters. 

We include cc resonance effects in i? — > Xsi~^£~ spectra with the help of a VMD model and present 
the distributions including next-to-leading order perturbative QCD corrections. 

We present l/mf, power corrections in the HQE approach and 0{as) perturbative corrections in the 
hadron spectra and hadronic spectral moments in B ^ Xsi~^i~ decays. 

The explicit dependences of the lowest moments of the hadronic energy Eh and the hadronic in- 
variant mass Sh, {Ejj), {Sjj) for n = 1,2 on the non-perturbative HQE parameters are worked 
out. We find that the first two moments of the hadronic invariant mass in B ^ Xgi^i^ decay are 
sensitive to A and Ai . 

We complement this profile of hadron spectra and moments in B —^ Xs£^i^ decays by an anal- 
ysis in the FM. The hadronic energy spectrum is found to be stable against a variation of the FM 
parameters, however, the hadronic invariant mass distribution depends significantly on them. 
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We incorporate the charmonium resonances by means of the VMD model into the analysis of hadron 
spectra and moments in B Xsi~^£~ decays. The broadening of the resonances in the FM in the 
hadron spectra is worked out. Also the cc resonances turn out to be important in the moments. 

We investigate the resulting uncertainties in spectra and moments in B ^ Xgi^i^ decays from 



different parametrizations [ ]28|J29| , |18[ ] of the resonant and non-resonant cc contributions. 



We work out hadron spectra, spectral moments and branching ratios in the FM with kinematical 



cuts, as used in the CLEO analysis in their search for the decays B Xsi^i^ , i = e, fj, [3Q]. 
They imposed two kinds of cuts: one in the dilepton invariant mass to exclude the main bulk of 
the J/V', V'') • • • resonances and another one in the hadronic invariant mass Sh < 3.24 GeV^, to 
suppress the BB background. Our study of the B — > Xs£~^i~ branching ratios with cuts is of direct 
use to estimate the efficiency of the remaining signal. 

To summarize, we shall present spectra, decay rates and hadronic spectral moments in inclusive rare 
B —^ Xsi^i~ decays in the standard model. The dilepton mass spectra and the FB asymmetry presented 
here are of use to extract the short-distance coefficients. We show that the uncertainties lying in different 
parametrizations of the charmonium resonance effects are not the dominant ones. Hadron spectra and 
moments in B Xsi^£~ decays can be used to test HQE and FM and/or to determine their parameters. 
We point out that the HQE, where it is valid, and the FM show very similar behaviour Moreover, some of 
their parameters can be related. In particular, the moments of the hadronic invariant mass inB —^ Xgl'^l" 
decays provide a complementary constraint on the non-perturbative HQE parameters A, Ai to the one in 



charged current semileptonic B —>■ Xlvi decays [|31|]. These decays can be used for a precise determi 



nation of these parameters. A related issue is the question of universality of the FM/HQE parameters for 
b q transitions with final quarks q = u,d,s or q = c. This remains to be tested. Finally, relating 
partly integrated hadron specti^a of B —^ Xsi^i^ to semileptonic B —^ X„£z/£ decays, we expect the 
cancelation of uncertainties resulting from bound state effects, thus this offers the possibility of a precise 
determination of Vub- 

Further, we study the exclusive mode Bs — > 77, which in the SM has a branching ratio in reach 



of future i?-facilities [^. We improve earlier analyses [32-34] by including leading logarithmic QCD 
corrections to the short-distance b — > 577 amplitudes. We use the same effective Hamiltonian as for 
b —>■ 57, which, as we will show, contains a complete set of operators for both decays. The Bs —>■ 77 
decay rate gets enhanced under renormalization, like the B — > Xsj decay rate [^. Likewise, we obtain 
a strong dependence on the renormalization scale of the Bs —>■ 77 observables, the branching ratio and 
the CP ratio [^, the latter resulting from CP-odd and CP-even parts in the FCNC 2-photon amplitude. 
Moreover, modeling Bs bound state effects in a HQE inspired approach, we avoid the constituent quark 
mass value nis ~ rrix, as used in previous analyses. 

In our analysis of Bs —>■ 77 decay we take into account long-distance (LD) effects via intermediate 
(neutral) vector mesons. Especially the contribution due to Bs 07 and subsequent decay — > 7 is 
estimated and found to sizeably reduce the Bs — > 77 branching ratio. We use QCD sum rules to evaluate 
the Bs 07 form factor and include the contribution from the gluon condensate. The VMD model is 
employed for the 0-meson photon conversion. 
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Finally, we perform a VMD based calculation of the decay Bg —>■ (pj ^ 77, where we abbre- 

viate = {J/tp,il^' , ■■■)■ We compare the LD-contribution to Bg —>■ 77 decay resulting from intermediate 



^p production with the one obtained by the interaction of the virtual charm loop with soft gluons [36]. 
We find that both amplitudes are in good agreement within the accuracy of the calculation. The contri- 
bution due to intermediate charmonium resonances changes the Bg —>■ 77 branching ratio including the 
intermediate contribution by less than 1%. 



Our work in Bg —>■ 77 decays [37,3a] can be summarized as follows: 



• We present leading logarithmic QCD corrections to the short-distance amplitude in 6 — > S77 decay. 

• We model Bg bound state effects in a HQE inspired approach, in contrast to the constituent quark 
model, which is used in the literature. 

• We estimate the LD effect due to the decay chain Bg —>■ cfry —>■ 77, using QCD sum rules and VMD. 

• We estimate the LD contribution through Bg — > (pip — > </)7 — > 77 using VMD. 

We present the branching ratio and the CP ratio for Bg — > 77 decay in the SM, taking into account 
improved perturbative 0{ag) contributions and long-distance effects via intermediate vector mesons. Un- 
certainties resulting from the renormalization scale and the bound state parameters are worked out. 

Organization of the work 

An introduction to rare B decays and the methods used is given in Chapter ^, where we discuss the 
effective Hamiltonian theory, rare radiative B —>■ Xg-y decays and long-distance methods in inclusive B 
decays (HQET, FM and VMD). Chapter ||is based on refs. [18,25,|^Q]. Here we investigate inclusive 



B Xgt^i decays in the SM. We present branching ratios and various spectra, furthermore hadronic 



spectral moments are estimated. In chapter ^ the exclusive channel Bg — > 77 is analysed [ p7| , [38| ]. SM 
based branching ratios and CP ratios in Bg 77 decay are given and their uncertainties are worked out. 
Finally, chapter |5| contains a summary and an outlook. Input parameters, Feynman rules and utilities are 
collected in appendix The power corrections to the structure functions of the hadron tensor in S — > 
Xgi'^l~ decays are given in appendix ^, together with auxiliary functions and the FM double differential 
Dalitz distribution in the context of the dilepton mass spectrum and the FB asymmetry. Appendix ^ 
contains analytic expressions used in the derivation of hadron spectra and hadronic spectral moment in 
B Xgi+i- decays. 



Chapter 2 



Rare B Decays: Motivation and Methods 



In this chapter we outUne the flavour structure of the standard model (SM). We discuss the CKM mixing 
matrix and motivate the importance of studying flavour changing neutral current (FCNC) 6 — > s tran- 
sitions. We introduce the necessary tool to include QCD perturbative corrections in weak decays, the 
effective Hamiltonian theory. As an application of the former we discuss B Xg'y decay as the most 
prominent example of a rare B decay. Finally non-perturbative methods like the heavy quark expansion 
technique, the Fermi motion model and vector meson dominance are sketched. 



2.1 The Flavour Sector in the Standard Model 



In the quark sector of the SM, there are six quarks organized in 3 families. The left-handed quarks are put 
into weak isospin SU{2)l doublets 



Qup 



(2.1) 

i=l,2,3 \" / L \" / L \" / L 

and the corresponding right-handed fields transform as singlets under SU {2)l. Under the weak interaction 
an up-quark (with Qu = 2/3e) can decay into a down-quark (with Qd = — l/3e) and a boson. This 
charged current is given as 

( d\ 



J. 



cc 



\/2 sin 9w 



CKM 



(2.2) 



where the subscript L = (l — 75)72 denotes the left-handed projector and reflects the V — A structure of 
J^*^ in the SM. Here the weak mixing (Weinberg-)angle 9w is a parameter of the SM, which is measured 
with high accuracy [|39|]. The so-called Cabibbo Kobayashi Maskawa (CKM) matrix Vckm [§I describes 
the mixing between different quark flavours. It contains the angles describing the rotation between the 
eigen vectors of the weak interaction {q') and the mass eigen states (g) 

( d'\ ( d\ 



s' 



Vckm 



s 

\b J 



(2.3) 
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Symbolically, Vckm can be written as 



Vr 



CKM 



^ Vud 

Vcd 



Vus 
Vts 



Vub ^ 
Vcb 

Vtb ) 



(2.4) 



In general all the entries are complex numbers, only restricted by unitarity VckmVckm ~ ^- They are 
parameters of the SM and can only be obtained from an experiment. Note that only three independent real 
parameters and one phase are left after imposing the unitarity condition. Some parametrizations of Vckm 
can be seen in ref. [|^]. 

A useful parametrization of the CKM matrix has been proposed by Wolfenstein [ pO| ] 

-A 
V ^A3(1 - 



^'Wolfenstein 



AX^{p - irj) \ 



iA2 



iri) 



2 



AX^ 
1 



+ 0(A^ 



(2.5) 



/ 



The parameters A, A, p and the phase r] are real numbers. A is related to the Cabibbo angle through 
A = sm9c [^], which describes the quark mixing with 4 quark flavours. Since A ~ 0.221, the relative 



sizes of the matrix elements in eq. ( |2.4| ) can be read off from eq. ([2.51). As can be seen, the diagonal entries 
are close to unity and the more off-diagonal they are, the smaller is the value of their corresponding matrix 
elements. The parameter A has been determined from the decays b ciu£ and B D*£i'g, yielding 
A = 0.81 ±0.07. The measurement of the ratio \ Vub/Vcb\ = 0.08 ±0.02 yields ^/p^Trf = 0.36 ±0.09. 
Likewise the mass difference AM^ = M{B^p) — M{Bj^^) ~ 0.46 (ps)"^ constrains the combination 



— p)^ + r?^. The observed CP-asymmetry parameter ex = 2.26 • 10"'^ constrains p and r]. The 
precise determination of the parameters p and t] is a. high and important goal, since it corresponds to two 
important questions: 

• Does CP hold in the SM ?? A non zero phase 77 7^ in the CKM matrix directly leads to CP 
violating effects. 

• The unitarity of the CKM matrix can be used to write down relations between its elements J2^=i ^ij^jk 
= 6ik, i,k = 1, 2, 3. There are 6 orthogonality equations possible (i ^ k), and each can be rep- 
resented graphically as a triangle, a so-called unitarity triangle (UT) pT]]. The sides and angles of 
such an UT can be constrained by different types of experiments. For the UT given by the relation 



VudV:, + VusV* + vM = o 



(2.6) 



there are 3 scenarios possible, which at present are not excluded experimentally and are a sign for 
new physics: 1. the UT does not close, i. e., J2i=i o^i 7^ 0' where Oj denotes the three angles of 
the triangle. 2. X]?=i ct* = 0' but the values of the are outside of their SM ranges determined by 
another type of experiment 3. J2i=i = 0' but the values of the angles are inconsistent with the 
measured sides of the triangle. 

In the literature special unitarity triangles are discussed. A recent review over the present status on the 
CKM matrix and the unitarity triangle is given by [H]. 
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2.1.1 Flavour changing neutral currents and why do we look for them in the 5 -system ? 

In the SM, the neutral current mediated through the gauge bosons Z*^, 7, g does not change flavour. There- 
fore, the so-called Flavour changing neutral currents (FCNC) do not appear at tree level and are described 
by loop effects. The subject of the present work is an analysis of such rare (FCNC mediated) B decays 
in the SM. The quarks are grouped into light (n, d, s) and heavy (c, 6, t) ones in the sense, that the mass 
of a heavy quark is much larger than the typical scale of the strong interaction, Aqcd 200 MeV. The 
sixth quark, the top, is too heavy to build bound states because it decays too fast. The special role of the 
6-quark is that it is the heaviest one building hadrons. We will not discuss the "double" heavy Be and con- 
centrate on i? = (bq) meson transitions with light q = u,d, s. Since the 6-quark is heavy, the i?-system 
is well suited for a clean extraction of the underlying short-distance dynamics. Unlike the iT-system, 
long-distance effects are expected to play a subdominant role in B decays except where such effects are 
present in a resonant form. 

The motivation to investigate b s{d) transitions is to improve the knowledge of the CKM matrix 
elements and to study loop effects. For the latter the interest is large, since there is no tree level FCNC 
decay possible in the SM. The leading loops give the leading contribution and they are sensitive to the 
masses and other properties of the internal virtual particles like e. g. the top. They can be heavy and 
therefore can be studied in a rare B decay at energies which are much lower than the ones necessary for 
a direct production of such particles. The idea is to compare the SM based prediction for a rare B decay 
with an experiment. A possible deviation gives a hint not only for the existence, but also for the structure 
of the "new physics" beyond the SM. 

Further the S-system can be used as a testing ground for QCD, to check perturbative and non- 
perturbative methods. One example is the decay B —>■ Xsj, which can be described in the lowest or- 
der at parton level through b — > 57. As a 2-body decay, the photon energy in the 6-quark rest frame is 
fixed: = (m^ — m'^)/2mb for an on-shell 7. A possible non trivial spectrum can result from gluon 
bremsstrahlung 6 —>■ s^g and/or a non-perturbative mechanism, which is responsible for the motion of 
the 6-quark inside the meson thus boosting the distribution. Such a bound state effect can be incorporated 



with e.g., the Fermi motion model, see section [2.4.2| for a brief discussion. 



Some rare B decays have already been detected. The channel B K*^ has been measured from 
the CLEO collaboration some time ago however the most prominent example of a rare B decay is the 
inclusive B — > ^^7 [pi]], where Xs is any hadronic state with strangeness s = 1 and B is a mixture of 
B^ and 5° (5°). The branching ratios are found to be 

B{B K*7)CLEO ^ 4 2 ± 0.8 ± 0.6 • 10"^ , (2.7) 
13{B Xsjf^^^ = 2.32 ± 0.57 ± 0.35 • 10"^ . (2.8) 



Also the ALEPH collaboration has presented a preliminary result [42] 



B{hb X,7)^^^^^ = 3.29 ± 0.71 ± 0.68 • 10"^ , (2.9) 

where hi, is any 6 flavoured hadron originating from decays, h^X. 

The calculation of the exclusive mode introduces large theoretical uncertainties due to the hadronic 
matrix elements. The inclusive decay is under better control, leading to the following result in the recently 



10 



CHAPTER 2. RARE B DECAYS: MOTIVATION AND METHODS 



completed NLO calculation [||], Q 



B{B Xs-i)NLO = 3.50 ± 0.33 • 10"^ . (2.10) 



Comparing eq. ( [2 .81 ) and eq. ( |2. 10| ), the CLEO measurement is found to be 2a away from the theory, 
but the SM cannot be ruled out. One has to look for other decay modes, since improving the theoretical 
accuracy in S — > decay seems not at hand. After displaying the methods developed in i? ^ ^s7, 
two other rare B decays, B — > with i = e, ^ and the exclusive decay Bg 77 will be discussed 

in chapter || and chapter ^ of this work, respectively. Both candidates have branching ratios 10~^ 
which are in reach of future B experiments. The aim of this thesis is to analyse these decays within the 
framework of the SM and to present up to date predictions for measurable quantities (branching ratios, 
distributions, asymmetries, etc) as accurately as possible with the present available techniques. 



/ \ 
/ \ 
I t 



b c,t s 

Figure 2.1: A FCNC b ^ s diagram. 



A typical diagram for 6 — > s is displayed in Fig. 2.1 from where the CKM couplings can be directly 



read off. The amplitude T is the sum over all internal up-quarks 

r= ^ ; Xi = V^bV*. (2.11) 



i=u,c,t 



Using the CKM unitarity J2i=u c t -^j = ^^'^ the smallness of Vub yielding ^ A^, we arrive at 

T = Xt{Tt - Te) + Xu{Tu - T,) ~ Xt{Tt - T,) (2.12) 



for a 6 ^ s amplitude in the SM. In the Li-system the FCNC transition rates (c — > u) are much more 
suppressed due to an inbuilt GIM mechanism [|4^. Here we have 

i=d,s,b 

= V,bV:,{n-Ts) + V,dVM-Ts) , (2.13) 

in which the first term is CKM suppressed and the second one GIM suppressed since — <C m^r. 
The SM rates in the charm sector for decays such as —>■ 77, —>■ aie out of reach for present 
experiments. If one nevertheless finds something in the rare charm sector, it would be a direct hint for the 
desired physics beyond the SM. 
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2.2 The Effective Hamiltonian Theory 

As a weak decay under the presence of the strong mteraction, rare B decays require special techniques, 
to be treated economically. The main tool to calculate such rare B decays is the effective Hamiltonian 
theory. It is a two step program, starting with an operator product expansion (OPE) and performing a 
renormalization group equation (RGE) analysis afterwards. The necessary machinery has been developed 
over the last years, see [j^-p^, [ |4^ ] and references therein. 

The derivation starts as follows: If the kinematics of the decay are of the kind that the masses of the 
internal particles rrii are much larger than the external momenta p mf ^ p^, then the heavy particles can 
be integrated out. This concept takes a concrete form with the functional integral formalism. It means 
that the heavy particles are removed as dynamical degrees of freedom from the theory, hence their fields 
do not appear in the (effective) Lagrangian anymore. Their residual effect lies in the generated effective 
vertices. In this way an effective low energy theory can be constructed from a full theory like the SM. 
A well known example is the four-Fermi interaction, where the -boson propagator is made local for 

<C roy^ (q denotes the momentum transfer through the W): 

-i (^ + + , (2.14) 

where the ellipses denote terms of higher order in 1/ mw- ' Performing an OPE for QCD and electroweak 
interactions, the effective Hamiltonian for a FCNC b ^ transition in the SM can be obtained by inte- 
grating out W, 0. Up to 0{^j^) it is given as: 

Heffib ^ S7) = -4^At Ci{^l)0^{f^) , (2.15) 

where the weak couplings gw = sing^ collected in the Fermi constant G p 

Gf _ 9w 



(2.16) 



w 

-5 r^^^T-2 



Gf = 1.16639 • 10"^ GeV"^ . (2.17) 
The on-shell operator basis is chosen to be 



01 = {sLalfjbLa){cLpl^CLf3), 

02 = isLalf^bLf3)icLi37''cLa), 

q=u,d,s,c,b 
q=u,d,s,c,b 
q=u,d,s,c,b 



' We remark here that the original way was reversed: The main historical step was to extrapolate the observed low energy 
4-Fermi theory in nuclear /3-decay to a dynamical theory of the weak interaction with heavy particle exchange. 
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q=u,d,s,c,b 

e 

16^2 



^7 = -:r^So,cr^^{mbR + msL)bo,F'"', 



Os = r^SaT^pa^,u{mi,R + msL)hpG''^"' , (2.18) 
where L{R) = 1/2(1 =F 75). '^^lv = ■^[l^l^lu\ and a, /3 are 5*^7(3) colour indices. T°, a = 1 ... 8 are 



the generators of QCD, some of their identities can be seen appendix |A.2| . Here F'"^, 0""^^^ denote the 
electromagnetic and chromomagnetic field strength tensor, respectively. As can be seen from the operator 
basis, only degrees of freedom which are light compared to the heavy integrated out fields {W,t,(f)), remain 
in the theory. The basis given above contains four-quark operators O1...6, which differ by colour and 
left-right structure. Among them, the current-current operators Oi and O2 are the dominant four-Fermi 



operators. A typical diagram generating the so-called gluonic penguins O3...6 is displayed in Fig. 2.2. The 



operators O7 and 0% are effective h — > 57, h ^ sg vertices, respectively. All operators have dimension 



6. For h —>■ sv£ transitions the basis eq. ( |2.18[ ) should be complemented by two additional operators 
containing dileptons. They are discussed together with their corresponding Wilson coefficients in chapter 

s 

/ \ 




Figure 2.2: A gluonic penguin diagram. 

The coupling strength of the introduced effective vertices Oi is given by the (c-numbers) Wilson 
coefficients Ci{ii). Their values at a large scale fi = mw are obtained from a "matching" of the effective 



with the full theory. In the SM, the Ci{m\Y) read as follows 1 134] 



Ci,z.Amw) = 0, (2.19) 

C2{mw) = 1 , (2.20) 
^, , 3a;3-2a;2 _8a;3 _ 53^2 _^ 7^ 

CiM = 1(^31)4 1^"+ 24(x-l)3 ' ^'-'^^ 

CsM = i(^3T)4l-^+ 8(x-l)3 ' ^'-"'^ 

with X = m^/m^r. It is convenient to define effective coefficients C^^g{fi) of the operators Oj and 
Og- They contain renormalization scheme dependent contributions from the four-quark operators Oi,..6 in 
Ti-eff to the effective vertices in 6 — > 57 and b — > sg, respectively. In the NDR scheme , which will be 
^We recall that in the naive dimensional regularization (NDR) scheme the 75 matrix is total anti-commuting, i. e. {75 , 7/j} = 
0, thus L-jfj, = 7^7?. 
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used throughout this work, they are written as 

Cfifi) = Crifi) + QdC^ifi) + QdN^Ceifi) , (2.23) 
Cfip) = Csifi) + C^ip) . (2.24) 

Here Nc denotes the number of colours, Nc = 3 for QCD. The above expressions can be found from 




Figure 2.3: The diagram contributing to the one loop b 57, b ^ sg matrix element, respectively. 

evaluating the diagram shown in Fig. Contributions from Oi,,,4, which correspond to an 7^L 7^L 
like insertion, vanish for an on-shell photon, gluon, respectively. The Feynman rules consistent with these 



definitions are given in appendix A.3 



2.2.1 QCD improved corrections 

Our aim is now to include perturbative QCD corrections in the framework of the effective Hamiltonian 

theory. This can be done by writing down the renormalization group equation for the Wilson coefficients 

3 

At^C'i(/x) = -fjiCjin) , (2.25) 

where 7 denotes the anomalous dimension matrix, i.e., in general the operators mix under renormalization. 
Solving this equation yields the running of the couplings Cj(//) under QCD from the large matching scale 



(here = mw) down to the low scale fi « m^, which is the relevant one for 6-decays. Eq. ( 2.25 ) can be 
solved in perturbation theory = Anas : 

-Y,-,- = ^7^°) + (^)M^) + (2 26) 

aifi) = Q(^)W + J^a(^)« + ... . (2.27) 
The initial values of the above RGE are the Ci{mw), which in the lowest order in the SM are given in 



eq. ( pJ9[p:22D . 

Let us for the moment concentrate on the special case that 7 is a number. Then the lowest order 
solution is given by 

C{fi) = 77 2/3o C{mw) , (2.28) 



'with Ci = Ci{fi, g) we have equivalently t^-^Ci = (^i^ + Mg^ ^)Ci' 
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Q!s(/i) 



(2.29) 



which can be easily checked by substituting it into eq. ( 2.25| ). In the derivation we have used the QCD /3 
function, which describes the running of the strong coupling: 



with its lowest order solution 



d 



167r- 



167r2 



as{mw) 



l + /3o^ln( 



4n 



(2.30) 



(2.31) 



We see that our obtained solution eq. ( 2.28 ) contains all powers of as(^) ln(^^). It is called leading 
logarithmic (LLog) approximation and is an improvement of the conventional perturbation theory. In 
general such a QCD improved solution contains all large logarithms like n = 0, 1, . . . (here with /i = m^) 



a:(m,) ln™(^) 



(2.32) 



where m = n corresponds to LLog. A calculation including the next to lowest order terms is called next 
to leading order (NLO) and would result in a summation of all terms with m = n — 1 and so on. In the 
following we use the 2-loop expression for as(^) which can be always cast into the form 



47r 



1 



/?ilnln(/.VA^^^) 



"-QCDJ 

With Nf = 5 active flavours (note that we integrated out the top) and SU{Nc 
coefficients of the /? function are 



(2.33) 

3) the values of the 
(2.34) 

They are given in appendix A.2 for arbitrary N^, and Nf. The strong scale parameter Aqcd = ^qcd^^ 
chosen to reproduce the measured value of as{fJ,) at the pole. For as{mz) = 0.112, 0.117, 0.122 we 
have Ag^^ = 160, 214, 280 MeV, corresponding to the values of the input parameters listed in appendix 



/3o-y , 



A.l 



We recall that in LLog the calculation of the anomalous dimension and the matching conditions at 
lowest order, 7^'^^ cf^\'m\Y) is required. In NLO a further evaluation of higher order diagrams yielding 
7^^\ c\^\mw/) is necessary and in addition the hadronic matrix elements (Oj) have also to be known in 



In a general theory and also in the one relevant for rare radiative 6 decays given in eq. (2.15), the 
operators mix and the matrix 7 has to diagonalized. In the latter case the (8 x 8) matrix 7^^^ has been 
obtained by [|8[^ and the running of the Ci(/u) in LLog approximation cannot be given analytically. The 
LLog solution for the Wilson coefficients ready for numerical analysis can be taken from [p7|]. We display 
the Ci for different values of the scale ^ in Table 2.1. As can be seen, there is a strong dependence on 
the renormalization scale /i, especially for Ci and C^^. Other sources of uncertainty in the short-distance 
coefficients Q are the top mass and the value of Qs(m^). We keep them fixed to their central values given 



in appendix A. 1 
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C^(M) 


/i = m^y 


/i = 10 GeV 


= 5 GeV 


^l = 2.5 GeV 


^1 





—0.161 


—0.240 


—0.347 


Co 


1 


1.064 


1.103 


1.161 







0.007 


0.011 


0.016 







—0.017 


-0.025 


—0.035 


^5 





0.005 


0.007 


0.010 







-0.019 


-0.030 


-0.046 




-0.196 


-0.277 


-0.311 


-0.353 


^8 


-0.098 


-0.134 


-0.148 


-0.164 



Table 2.1: Leading order Wilson coefficients in the Standard Model as a function of the renormalization 
scale (1. 



Here a comment about power counting in our effective theory is in order: As can be seen from Fig. g3 
with an external photon, the insertion of four-Fermi operators generates a contribution to 6 ^ 57, which 
is also called a "penguin". It is a 1-loop diagram, but urdike "normal" perturbation theory, of order q^. To 
get the a], contribution, one has to perform already 2 loops and so on. That means, the calculation of the 
LO(NLO) anomalous dimension matrix was a 2(3)-loop task. 



A comprehensive discussion of weak decays beyond leading logarithms can be seen in ref. [|46|]. The 
main results of the NLO calculation in — > decay will be given in section 2.3, 



The advantages of the effective theory compared to the full theory can be summarized as follows: 

• The effective theory is the more appropriate way to include QCD corrections. Large logarithms like 
ln(^/mvi/) multiplied by powers of the strong coupling as{n), which spoil the perturbation series 
in the full theory, can be resummed with the help of the RGE. 

• On the level of a generic amplitude A = (Heff) ~ Z]j Ci{^){Oi){fi) the problem can be factorized 
into two parts: The short-distance (SD) information, which can be calculated perturbatively, is en- 
coded in the Ci, and it is independent of the external states, i.e. quarks or hadrons. The long-distance 
(LD) contribution lies in the hadronic matrix elements. Both are separated by the renormalization 
scale fjL. Of course the complete physical answer should not depend on the scale /x, truncating the 
perturbation series causes such a remaining dependence, which can be reduced only after including 
higher order terms or a full resummation of the theory. 

• As long as the basis is complete, the effective Hamiltonian theory enables one to write down a 
model independent analysis in terms of the SD coefficients Cj. This is true for SM near extensions 
hke the two Higgs doublet model (2HDM) and the minimal supersymmetric model (MSSM). Here 



one can try to fit the Cj from the data []12|]. However, new physics scenarios like, e.g., the left-right 
symmetric model (LRM) require an extended operator set [48,4^|^]. Wilson coefficients in the 
2HDM and in supersymmetry (SUSY) can be seen in ref. [51] and ref. [|52|], respectively. 
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2.3 B Xsj in the Effective Hamiltonian Theory 

The effective Hamiltonian theory displayed in the previous section is applied to i? ^ Xs^ decay. Several 
groups have worked on the completion of the LLog calculation [||,^. The anomalous dimension matrix 
at leading order 7*^'') and the lowest order matching conditions (eq. ( 2.19^2.2^ )) govern the running of 



the LLog Wilson coefficients, denoted in this and only this section by to separate them from the 

NLO coefficients. We discuss the improvement of the theory in i? — > obtained from NLO analysis. 
In the remainder of this work we treat the Wilson coefficients Q, i = 1, . . . 8 in LLog approximation. 
In the spectator model, the branching ratio for B —>■ Xgj in LLog approximation can be written as 

B(B ^ X.,) = B.^^l^ = j^-^^C-r^ipf . (2.35) 

T{b ^ ceue) \Vcb\ TTJimc) 

where a normalization to the semileptonic decay B — > Xcive to reduce the uncertainty in the fe-quark 
mass has been performed. Here Bsi denotes the measured semileptonic branching ratio and the phase 



space factor /(rn-c) with rhc = rric/mb for T{B Xc^-ff) can be seen in eq. ( 3.30 ). 

'7 



As the branching ratio for B —>■ X^"/ is mainly driven by c!f''^^^{fj.), several effects can be deduced: 
Including LLog QCD corrections enhance the branching ratio for B Xg^ about a factor 2 — 3, 



as can be seen in Table |2jJ (here denoted by Ci{fi)) and changing the scale from /i = niw down to 
/X ~ rub. 

• While the sign of c!f^^^^ is fixed within the SM, i.e. negative, it can be plus or minus in possible 
extensions of the SM. A measurement of B{B —^ Xsj) alone is not sufficient to determine the sign 

, or in general, the sign of Cj^^ resulting from possible higher order calculations. 

• The strong scale dependence of the value of df^^^^dj,) causes serious problems in the accuracy of 
the LLog prediction. Varying the scale between ^ < < 2mb , results in an error in the branching 
ratio of ±25% [|3]l, 0. 

Because of the last point the NLO calculation was required. Several steps have been necessary for a 
complete NLO analysis. Let us illustrate how the individual pieces look like: At NLO, the matrix element 
for 6^57 renormalized around fi = nib can be written as |0] : 

M{b ^ S7) = -A^XtD{07{mb))tree , (2.36) 

with 

D = Cf ff(M) + ^ E (cP'"M-,P In ^ + . (237) 

The rij are computed in ref. [43|. They contain the bremsstrahlung corrections [|6|], [ ^ ] b —>■ sjg and 



virtual corrections to the O7 matrix element [43]. Especially the latter with an O2 operator insertion 



demands an involved 2-loop calculation, see Figs. 1-4 in [43], where the corresponding diagrams are 
shown. It is consistent to keep the pieces in the parentheses in eq. (2.37), which are multiplied by as{mb), 
in LLog approximation. 
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Now Cj^^{fi) has be be known at NLO precision, 

Cf%) = + ^4^C«^^%) . (2.38) 

As this job consists out of two parts, the work has been done by two groups: The 0(a1) anomalous 
dimension matrix was obtained in ref. [p4|], which required the calculation of the residue of a large number 
of 3-loop diagrams, describing the mixing between the four-Fermi operators O1...6 and Oj^s- The second 



part, the NLO matching at /i = niw has been done in ref. [ p^ ] and confirmed in ref. [p6|]. The NLO 
calculation reduces the fi = 0{mh) scale uncertainty in varying fi in the range ^ < ji < 2mi, drastically 



to ±4.3% [57] and suggests for B —>■ Xg'^ a scale /x = ^ as an "effective" NLO calculation through 



FiB ^ Xsj)Loif^ = ^ ^ ^s7)nlo ■ (2.39) 

As a final remark on scale uncertainties it should be noted that in the foregoing the top quark and 
the W have been integrated out at the same scale jj, = mw, which is an approximation to be tested. It is 
justified by the fact that the difference between {my/) and {nit) is much smaller than the one between 
as(miy) and as(mf,). ^ The authors of [57] analysed the dependencies on both the W matching scale 



= 0{mw) and the one at which the running top mass is defined: rht{fJ-t) and jj-w / fit. Similar to 
the rub scale they allowed for fiy/^nt a possible range: ^ < /U^ < 2mx where x = W,t. Their findings 
are that the l^t uncertainty is much smaller (namely ±1.1%, ±0.4% at ^ ~ nib in NLO, respectively) 
than the uncertainty in the scale around nib and therefore negligible. 



2.4 Long-Distance Effects in Inclusive B Decays 

In this section we sketch the methods to treat the LD effects in inclusive B decays. We have effects due 
to the confinement of the quarks in a bound state and due to resonances. They will be explained more 
detailed in the following chapters when and where necessary. For the evaluation of the exclusive channels 
we refer to chapter ^ in which the rare mode Bg —>■ 77 is discussed, especially section i2. 



There are mainly two different approaches to take into account the effects of the S-meson bound state, 
the heavy quark expansion (HQE) and the phenomenologically motivated Fermi motion model (FM). 
While the former is a field theory in the framework of the heavy quark limit of QCD and has an interest of 
its own, the latter serves as a model of the data and has no intrinsic problems like end-point singularities 
etc. like HQE. Both models have parameters which can be related to each other and as they are used as 
inclusive methods, no form factors appear in the amplitudes. Inclusive decays are good from theoretical 
point of view and a challenge for experimentalists: An inclusive final state X is an average over a suf- 
ficiently high number of exclusive single (resonances) and continuous multi body states with the same 
quantum numbers as X. Inclusive decays involve the calculation of quark level processes. The underlying 
assumption of quark-hadron duality requires a large and dense enough populated phase space. By means 
of a "smearing" procedure, the singular behaviour in a local form is avoided and the differential spectra 
can be measured in a distribution sense. 



Using eq. (2.33) and the input parameters in Table /^.l, we have as{mw) ~ 0.12, Qs(mt) = 0.11 and as{mt) = 0.21. 
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Another type of LD effect beside the bound state effect mentioned above, is due to resonances. A {qq) 
spin 1 state can hadronize from a virtual g-loop. The conversion of such a vector meson into a photon is 
described by the phenomenological vector meson dominance (VMD) model. 



2.4.1 The heavy quark expansion of QCD 

Consider a hadron containing one heavy quark Q in the limit mq — > cxd. The other ingredients, light 
quarks and gluons, are seen as a light cloud around the heavy quark, sometimes also called the "light 
brown muck" [58], which exchanges small momenta of order Aqcd with Q for which a perturbative 
expansion is not useful. The parameter Aqcd characterizes the soft hadronic interaction scale. The heavy 
quark inside the bound state is treated as a static source of gauge charge (colour and electric charge): Q is 
so heavy compared to Aqcd, that it does not recoil as a result of the soft exchanges, it sits at rest in the 
hadron rest frame. This is the heavy quark limit of QCD. 

New symmetries can be explored which are exact in the limit of an infinitely heavy quark: The light 
degrees of freedom are insensitive to the mass, flavour and spin of the heavy quark! This brings an enor- 
mous simplification of certain aspects of QCD, like the calculation of heavy quark matrix elements and 
hadron spectroscopy. However, for a firm phenomenological analysis we need to go from the predictions 
of heavy quark symmetry in the strict limit uiq ^ cxd to a theory which provides a controlled expansion 
around this (academic) case. This can be done with the heavy quark expansion technique (HQET) in the 
limit mq » Aqcd- The necessary technology has been developed over the last decade and can be seen 
in a selection of papers [15 16] and references therein. A nice review on the HQE technique is given in 
ref. pi]. 



l/mh power corrections 

Let us now switch to the system under consideration, the S-mesons. Since the 6-quark is heavy, i. e. 
rrih ~ 4.8GeV » Aqcd ~ 200 MeV, the success of the spectator model in B = (bq) meson decays can 
be understood. Moreover, corrections in inverse powers of mi, to this can be systematically obtained with 
the help of the HQET. 

The light degrees of freedom in the i?-meson give rise to the parameter A which accounts for the 
binding energy of the bound state. In the limit of an infinitely heavy 6-quark, i. e. ^ oo we have 

mB = mb + A, (2.40) 

where ms denotes the i?-meson mass. Corrections to this can be calculated within the following set 
up of HQE: The heavy 6-quark momentum is written as = mi^v + k, where k is a small residual 
momentum of order Aqcd- v denotes the velocity of the meson with momentum P = niBV, which at 
rest is u = (1,0, 0, 0). It follows that the relative movement between the heavy quark and the meson is 
suppressed by powers of k/mi). Performing an operator product expansion up to operators with canonical 
field dimension 5, the HQET mass relation modifies to [ ]60| ] 

niB = mb + A- -^(Ai + 3A2) + . . . , (2.41) 
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Reference 


Method 


A [GeV] Al [GeV2] 


Falk et al. [|6|] 


Hadron Spectrum 


ssO.45 « -0.1 


Gremm et al flvill 
Chemyak ^ 
Gremm, Stewart |68] 


T pntnti Snpctnim 

{B ^ Xii)) 


n39±nii — niQinio 

. (J ty —1— v_f . i X yj , ij _i_ yj , j^yj 

0.28 ±0.04 -0.14 ±0.03 
0.33 ±0.11 -0.17 ±0.10 


Li, Yu [69] 


Photon Spectrum 


0(tK+0A2 _n 71 +0-70 



Table 2.2: Determinations of the parameters A and \\from inclusive decay spectra. 



Ukewise for B* vector mesons 



TUB* 



mb ± A 



(Ai-A2)±... , 



(2.42) 



where the ellipses denote terms higher order in I /mi). In general, the next to leading power corrections in 
HQET are parameterized in terms of these matrix elements of the kinetic energy and the magnetic moment 
operators Ai and A2, respectively. We can get the value of A2 from spectroscopy 



A2 ~ 



0.12 GeV^ . 



(2.43) 



The quantity Ai is subject to a theoretical dispersion. Its value has been determined from QCD sum rules. 



yielding Ai = -(0.52 ± 0.12) GeV^ (Ball and Braun in [51]) and Ai = -(0.10 ± 0.05) GeV^ (Neubert 



2|]). Further, the value for Ai has been extracted from an analysis of data on semileptonic B decays 
(B — > Xlui), yielding Ai = —0.20 GeV^ with a corresponding value A = 0.39 GeV, as the two are 



correlated [31]. For a review on the spread in the present values of these non-perturbative parameters 



extracted from inclusive decay spectra, see Table |2.2|, which is adopted from [p3|]. 

Now there is an intrinsic difficulty in the relation eq. ( 2.41[ ). One can ask for the meaning of the "pole" 
mass rri}, and A ? First of all, they are non-perturbative parameters and they add up in the combination 



given by eq. (2.41) to the physical B-meson mass. However, while the sum is fixed, there is a scheme 
dependent "renormalon" ambiguity of order Aqcd in both nih and A, which cancels out in physically 



measurable quantities [64]. Assuming universaUty, the parameters nib, A determined by one experiment 
can be used to help the analysis of another decay, provided that one uses the same renormalization scheme 
prescription. 

The power corrections in i? — > Xsj decay including I /ml terms have been calculated in ref. [17]. 
The l/mf corrections have been recently reported in ref. []65|]. However, the use of the last calculation is 
limited by the fact, that the matrix elements of the higher dimensional operators are almost unknown. In 
B —^ Xsi^i" decay the 1 / m^ have been first calculated in ref. [|l^] (with massless s-quark) and corrected 



in ref. [18] with full m^. The latter has been confirmed recently in the massless s-quark case in ref. [27] 
Details of the HQET calculation are given in chapter H. 
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l/rric power corrections 



With rUc ~ 1.4GeV » ^qcd the charm is still a heavy quark and also l/rric power corrections are 
subject of present i?-physics. A theoretically interesting structure, an effective hs^g-vextex appears from 



the diagrams displayed in Fig. 2^ [36]. The amplitude of this operator can be expanded in l/m-c. The 
power corrections in Aqcd/ w-c to the B Xsj decay rate have been calculated first in ref. [ [70| ] and ^J%, 
however with the wrong sign. This has been settled now in favor of [72,7^,27]. The resulting correction 
to the decay rate is found to be small: 



6riB ^ X,7) 



2^2 



+0.03 



(2.44) 



V{B ^ X,7) 

The 1 / mc expansion is an alternative description of the virtual charm loop to the traditional vector meson 





Figure 2.4: The hs^g-vertex at lowest order in QCD. 



dominance ansatz, (see section 2.4.3 , Fig. 2^) in regions of momentum transfer far away from the 
resonances. 

Pros and cons of the HQE approach 

• systematic expansion in AQC£i/mc,fe 

• well-defined limit of QCD 

• unknown matrix elements of higher order operators 

• expansion breaks down in some regions of the phase space 

The last point above corresponds to an uncomfortable property of the HQE method: The development of 
end-point singularities in inclusive decay spectra [T^18|, hence only quantities smeared over a sufficiently 
large phase space interval are calculable. Such a smearing is incorporated, e.g., in the Fermi motion model 



[22], discussed in the next section. 



2.4.2 The Fermi motion model 

Another model to handle the effects of the bound state is the phenomenological Fermi motion model (FM) 



[22|. The FM is defined through the requirement that the 6-quark and the spectator quark q four-momenta 
add up to the B = (6g)-meson four-momentum. In the rest frame of the S-meson the quarks fly back to 
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back with momentum p = Pb = —Pq- From energy conservation follows than that either the spectator or 
the h quark has to have a momentum dependent mass. In this work we choose niq to be a parameter of the 
model and the b quark to have a p dependent mass. It reads 



mf,(p) = m_B + rrig - 2mB\Jp'^ + rriq^ ; p=\p\. (2.45) 

The next assumption is that the momentum p obeys the Gaussian distribution function (j){p) weighted with 
the Fermi momentum pp 

4>{p) = exp(^) , (2.46) 

y/T^PF-^ PF^ 

with the normalization /q°° dpp^ cl){p) = 1. The procedure to implement these wave function effects to a 
general parton model distribution obtained in the b quark rest frame is as follows: 

1. replace the b quark mass by mb{p) 

2. boost the distribution into the S-meson rest frame and 

3. fold the result with the wave function given in eq. ( 2.46| ). 



For subsequent use in working out the normalization (decay widths) in the FM model, we also define 
an effective 5-quark mass by 

= ( r dpp^ mb{pf^{p)Y/^ . (2.47) 
Jo 

The two parameters of the FM model, the Fermi momentum pp and the spectator mass niq can be 
fitted from data, however, up to now this procedure has not been very conclusive as still large ranges of 
the parameters are possible. The question appears here whether the FM parameters do depend on flavour, 
i. e., are they universal for B —^ Xj + {'yjeptons) with / = d,u,s,c transitions. Further relations 
between the FM parameters {pp, niq) and the HQET parameters (Ai, A) can be obtained. However, there 
is no analogue of the magnetic moment coupling A2 in the FM. 

We will return to the FM in chapter ^ to model the wave function effects in B ^ Xst^l~ decay. 

2.4.3 Vector meson dominance 

Vector meson dominance (VMD) provides a mechanism to convert a spin 1 meson, here generically de- 



noted hy V = (qq) into a photon [74]. The creation of a (qq) bound state from a virtual q, q pair and its 



subsequent conversion into a photon is displayed in Fig. 0|. The intermediate vector meson propagator 




V 



7 

IPwwwv 



Figure 2.5: The VMD conversion. 
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equals l/(my — — imi/Fy*), where the imaginary parts takes the effect of a finite total width Fy* of 
the V into account. The matrix element of the constituent current is defined as 

< 0|^7^g|y(g,e^) >= fv{q^)mvel , (2.48) 

and the photon interacts with the (neutral) bound state through the electromagnetic current J^*" = eQqqj^q 
as: 

< 0|J^-|y(g,e^) >= eQJviq^)mve^ • (2.49) 

Here q, are the momentum and the polarization vector of the vector meson V, respectively and Qq 
denotes the charge of the quark q in units of e. VMD conversion means that while y — > 7 also e^, 
where is the polarization vector of the photon. The form factor at q"^ = rriy can be obtained from 
measurement of the leptonic width 

r v^e+e- =/y^^| . (2.50) 

What about fv{q^) at q^ 7^ my ? For the calculation of the form factor at < < my, there 
exist different ways, which are neither straightforward nor unique, a) First of all data on photo production 
7A^ ^ yiV at g2 = can be used. They indicate a large suppression of /i/(0) compared to fvi'm'y) 
only if the vector meson is heavy such as y = J/i/jjip'. As in this case = nT^jj^^rn^i to = 
involves a large extrapolation in q^. It is not expected to be so significant for F = p, w, as = etc. 
is not far away from q^ = 0. This is what all the methods listed here share, that fv{q^) decreases with 
decreasing q^. b) Often a single-pole form is assumed to extrapolate the form factor to smaller values of 
the momentum transfer 

fv{q') = , , (2.51) 

1 - 97%oZe 



where mp^ie corresponds to the masses of higher resonances of the V . c) The approach by [ [75| ] is based 
on a dispersion relation calculation. It yields an interpolation formula between /i/(0) and fvi'm'y) where 
both can be fixed by data. The situation for q^ > my is unclear. 

The VMD mechanism has been applied to i? ^ Xg'y decay to estimate the long distance contribution 
through B Xsip' X^j by Here = J/i^,i^', ■ • • V'^"^ are the known six charmonium 

resonances, see p9|]. The -0* ^ 7 conversion requires the knowledge of the form factor at = for an 



on-shell photon. It has been shown in ref. [76], that the methods a) and c) listed above yield a consistent 
suppression at q"^ = 0. The longitudinal degrees of freedom of the have been removed using the 



procedure proposed by |77]. However, assumptions made remain as uncertainties in the calculation. 

We will make extensive use of VMD in chapter ^ and chapter ^ to include long-distance effects from 
intermediate vector mesons in B Xgi^i^ and Bg — > 77 decays, respectively. However, in chapter ^ 
we estimate the resulting uncertainties which emerge from various theoretical approaches in implementing 
the -dependence of the VMD-domrnated amplitude. 



Chapter 3 



Inclusive B Xgi i Decay 



This chapter contains a comprehensive analysis of B ^ Xs£~^i~ decay in the standard model (SM). We 
include QCD improved 0{as) corrections, use heavy quark expansion techniques (HQET) and apply the 
Fermi motion model (FM). Further, the long-distance effects via intermediate J/ip, ip', ... resonances are 
taken into account with a vector meson dominance (VMD) ansatz. 



3.1 Introduction 



Flavour changing neutral current (FCNC) decays B Xgi^i^ and B Xgj are governed in the SM by 
loop effects. They provide a sensitive probe of the flavour sector in the SM and search for physics beyond. 
In the context of rare B decays the radiative mode B — > has been extensively discussed in chapter |2[ 
In this chapter we address inclusive B — > Xgi'^i'^ decay with £ = e, //. Since we are neglecting 
finite lepton masses we cannot apply our results to the r-case. The b — > s£~^i~ transition has been 
studied earUer in the free quark model in refs. 80 1 in the lowest order in the SM context. The NLO 
0{as) improvement in the invariant dilepton mass distribution and the decay rate has been worked out in 
refs. [p|,p7[]. Leading (1/m^) power corrections in the HQET framework [15,16] in the invariant dilepton 
spectrum in i? — > Xsi~^£~ decay have been reported in ref. [18], correcting an earlier calculation [17]. 



This has been recently confirmed in ref. [27] for the massless s-quark case. Another interesting quantity 
in B ^ Xsi^£~ is the FB asymmetry [13], also known at l/ml [p^]. It can be used together with 
the branching ratio of i? — > and the dilepton spectrum in B —>^ Xgl'^l^ for a model independent 



analysis of the short-distance coefficients [ p^ ] in the search for SUSY effects [pT[-83], [|50|]. The 1/ 



imr 



power corrections to the left-right asymmetry [84,85] have been presented in [Bg] correcting an earlier 



calculation of the same [85]. Both the FB asymmetry and the left-right asymmetry are defined in section 
p.l.3| . The longitudinal polarization of the lepton, Pl, in B ^ X^r+r^ at the partonic level has been 



worked out [87]; the other two orthogonal polarization components Pt (the component in the decay 
plane) and Pl (the component normal to the decay plane) were subsequently worked out in [p8[]. The 
0{l/m1) correction to the dilepton invariant mass spectrum in — > Xgl'^l" has also been calculated 



in []88[], however, the result differs in sign from the one in [[27|]. This controversy, which goes back to 
the corresponding one in i? ^ ^^7 decay (see section 2.4.1 for a discussion) has been settled in favor 
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of [^]. It is known that inclusive decay spectra are not entirely calculable with HQET [19-21], [16], 
especially, the expansion in powers of l/nib diverges in the high dilepton mass region in i? ^ Xsi^i" 
decay [111]. An alternative approach to take into account i?-meson bound state effects is the FM model 



[E2p. In the FM a prediction for the entire range for the dilepton mass distribution and FB asymmetry 



in B Xgi I decay has been given ref. [18]. The sensitivity of the distributions on the FM parameters 
is worked out there. Long-distance effects due to intermediate B ^ Xs + {J /'4^, ip' , •••) Xsl'^i^ have 



been discussed in refs. [ p9|j2qj9(]| ] and recently [|T4p. Hadron spectra and hadronic spectral moments are 
presented in refs. [ ^J2^JT4| ] in both the HQE approach and the FM. 

This chapter is divided into two parts. The first one (this section up to and including section based 
on ref. [p^], contains an introduction to i? ^ Xgl^i^ decay, basic definitions and the 0{as) and I /ml 
corrected matrix element for b — > s£~^i~. Is is mainly devoted to the analysis of the dilepton invariant mass 
distribution and the FB asymmetry. In doing that, we derive leading power corrections to the decay rates 
and distributions in the decay B Xst^l~ using heavy quark expansion (HQE) in (l/m^). Further, 
wave function effects of the 6-quark bound in the i?-hadron are studied by us in the phenomenologically 
motivated Gaussian Fermi motion model. Using this model for estimating the non-perturbative effects, 
we include the dominant long-distance (LD) contributions from the decays B ^ Xg + {J /'^l^, ■■■) ~^ 
Xs£~^i~. Further, taking into account the next-to-leading order perturbative QCD corrections in 6 — > 
si~^i~, we present the decay rates and distributions for the inclusive process B Xst^l~ in the SM. 

The second part, starting from section complements the study of i3 — > Xs£'^l~ decay and investi- 
gates hadron spectra and hadronic spectral moments. It is mainly based on refs. [25,26jl4]. We compute 
the leading order (in ctg) perturbative QCD and power (1/m^) corrections to the hadronic invariant mass 
and hadron energy spectra in the decay B Xs£~^i~. The computations are carried out using HQET 
and a perturbative-QCD improved Fermi motion model which takes into account S-meson bound state 
effects. We also present results for the first two hadronic moments (5^) and {E^), n = 1,2, working 
out their sensitivity on the HQET and FM model parameters. In the FM, also the LD effects due to in- 
termediate charmonium resonances are taken into account. We study uncertainties in the parametrization 
of the cc effects. Further, we investigated the effect of the experimental cuts, used recently by the CLEO 
collaboration in searching for the decay B —^ Xsi'^l^ pO|], on the branching ratios, hadron spectra and 
hadronic invariant mass moments using the FM model. 



3.1.1 Kinematics 

We start with the definition of the kinematics of i? — > Xs£~^£~ decay at parton level, 

b{Pb) ^ s{ps){+g{pg))+i+{p+)+e-{p-) , (3.1) 

where g denotes a gluon from the 0{as) correction (see Fig. We define the momentum transfer to 
the lepton pair and the invariant mass of the dilepton system, respectively, as 



q = P++P- 
s = (f' . 



(3.2) 
(3.3) 
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The dimensionless variables with a hat are related to the dimensionful variables by the scale rrib, the 
6-quark mass, e.g., 

s . rus 

s = ^ , ms = — , (3.4) 

nih 

etc.. Further, we define a 4- vector v, which denotes the velocity of both the 6-quark and the i?-meson, 
Pb = rubV and pb = msv. We shall also need the variable u and the scaled variable u = defined as: 

b 

u = -{Pb - P+f + {Pb - P-f , (3.5) 
u = 2v{p+-p^), (3.6) 

and further the kinematical phase factor 

u{s,ms) = \/is- {nib + msY){s - [nib - nisY) . (3.7) 

The scaled variables s and {t in the decay b s£'^£^ are bounded as follows, 

—u{s,ms) <u < +u{s,ms) , 
u{s,ms) = ^[s - (1 + ms)2] [s - (1 - rhs)^] , 

4mf < s < (1 - rhsf . (3.8) 

3.1.2 NLO-corrected amplitude for b — > si~^i~ 

Next, the explicit expressions for the matrix element and (partial) branching ratios in the decays b — > 
s£~^£~ are presented in terms of the Wilson coefficients of the effective Hamiltonian obtained by integrat- 
ing out the top quark and the bosons, 

Heffib ^ s + £+£-) = Heffib ^ s + 7) - ^KVtb [^9(^)09 + CloOlo] , (3.9) 

where Heffib —>■ s + ^) together with the operators Oi.,,8 and their corresponding Wilson coefficients 
Ci{fi) [0j6|] can be seen in section |Z2[ The two additional operators rnvolvrng the dileptons Og and Oio 



are defined as: 

o2 



= -^Sa7^Lbc,£-f^i , 

'^10 = T-^s,7^iM7M75^ • (3.10) 

A usual, CKM unitarity has been used in factoring out the product Vf*Vtb- Note that the chromomagnetic 
operator Og does not contribute to the decay B Xs£^£~ in the approximation which we use here. The 
Wilson coefficients are given in the literature (see, for example, [^^). They depend, in general, on the 
renormalization scale fi, except for Ciq. At leading logarithmic (LLog) approximation, we use the values 



of the Ci given in Table 3.1 



With the help of the effective Hamiltonian in eq. ( |3.9| ) the matrix element for the decay b si^£ 
can be factorized into a leptonic and a hadronic part as, 

M{b^s£+£-) = ^VZVtb[{cf-Cio){-s^^Lb){ij'^L£) 



26 



CHAPTER 3. INCLUSIVE B Xsi+l' DECAY 





C2 


C3 


C4 


C5 


c. 


cf 


Cg 


Cio 


C(o) 


-0.240 


+1.103 


+0.011 


-0.025 


+0.007 


-0.030 


-0.311 


+4.153 


-4.546 


+0.381 



Table 3.1: Values of the Wilson coefficients used in the numerical calculations corresponding to the 

.ejf. 



central values of the parameters given in Table AJ_. Here, C-f = C7 — C5/3 — Cg, and for Cg we use 
the NDR scheme and C(°) = 3Ci + C2 + 3C3 + C4 + 3C5 + Cq. 



-2Cf [sta 



(3.11) 



where we abbreviate cf = cf{s). We have kept the s-quark mass term in the matrix element exphcitly 
and this will be kept consistently in the calculation of power corrections and phase space. The above matrix 
element can be written in a compact form, 



v,:v,, (r^L^^ + r^L^'^) , 



with 



yL/r ^ 



R^^ (cf ^ Cio + 2Cf\\ + 2ms Cf^, |l 



(3.12) 

(3.13) 
(3.14) 



where we have already used massless leptons in substituting —2ia^^q'^ = [7/^)^ by 27^^ in the term 
proportional to cf. 

The effective Wilson coefficient cf{s) receives contributions from various pieces. The resonant cc 
states also contribute to cf{s) and will be discussed in section 3^; hence the contribution given below 
is just the perturbative part: 



cfis) = Cgr^is)+Yis) 



(3.15) 



The function Y{s) represents the one-loop matrix element of the four-Fermi operators [47,|], see Fig. 3.1 




Figure 3.1: The Feynman diagram responsible for the four-Fermi-operator contribution (depicted by the 
blob) to the operator O9. 
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It is written as: 



Y{s) = gimc,s){3Ci + C2 + 3C3 + Ci + 3C5 + Ce) 

-^5(1, s) (4 C3 + 4 C4 + 3 Cs + Ce) - ^9(0, s) (C3 + 3 C4) 

+1 (3C73 + C4 + 3C5 + Ce) - e ^ (3Ci + C2 - C3 - 3C4) 
y y 



(3.16) 



(3.17) 



(NDR) , 
-1 (HV), 

where (NDR) and (HV) correspond to the naive dimensional regularization and the 't Hooft-Veltman 
schemes, respectively. We recall that while Cg is a renormalization scheme-dependent quantity, this de- 
pendence cancels out with the corresponding one in the function Y{s) (the value of ^, see above). The 
function g{z, s) includes the quark-antiquark pair contribution |P,P7||: 



9{z,s) = -9M-) 



e(l-y)(ln 



-Inz H h -y (2 + y)\ \l 

9 27 9^ 9^ I 



y\ 



1 + 
1- 



ivr) + Q{y — l)2arctan 



= 27-9^"^^) 



4, . 4. 

— m s H — ITT , 
9 9 ' 



(3.18) 
(3.19) 



where y = Az^/s. As can be seen from the above equations, internal 6-quarks ~ (7(1,5), c-quarks 
~ g{rhc, s) and light quarks q, (with ruq = for q = u, d, s) ~ g{0, s) contribute to the function Y{s); 
only the charm loop involves the dominant "current-current" operators Oi and 02- 

The 0{as) correction [^] from the one-gluon exchange in the matrix element of Og in the invariant 
dilepton mass s is represented by 



1 H LOis) 



vr 



(3.20) 



where 



a;(s) 



2^4 2 
-gTT - -Li2(s) - -lnsln(l 



2s(l + 5)(l-2s) 
30" 



s)2(l + 2sl 



In s + 



5 + 9s 



5 + 4s 
~ 3(1 + 2s) 
6s2 



ln(l 



(3.21) 



Note that the function uj{s) is given with 



6(l-s)(l + 2s) 
0. The one-gluon correction to Og with respect to the final 



partonic energy and the invariant mass will be presented below in section 3.6 



In the order we are working only Og is subject to corrections since the renormalization group 
improved perturbation series for Cg is 0{l/as) + 0(1) + 0{as) + . . ., due to the large logarithm in Cg 



represented by 0{l/as) [p7|]. The Feynman diagrams, which contribute to the matrix element of Og in 



0{as), corresponding to the virtual one-gluon and bremsstrahlung corrections, are shown in Fig. 3.2 . 
With the help of the above expressions, the differential decay width becomes on using p± = {E± , p±), 



dr 



1 



2mB 27r 



GF'a\ 2 



(27r)32E+ (27r)32E_ 



(3.22) 
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Figure 3.2: Feynman diagrams contributing to the explicit order Us corrections of the operator Og. Curly 
lines denote a gluon. Wave function corrections are not shown. 

J /7-> 

where Wfj_y and L^y are the hadronic and leptonic tensors, respectively. The hadronic tensor W^l is 
related to the discontinuity in the forward scattering amplitude, denoted by T^l , through the relation 



W, 



2ImT^j,. Transforming the integration variables to s, u and v ■ q, one can express the Dalitz 



distribution in 6 ^ si^i (neglecting the lepton masses) as: 



dr 



dudsd{v 



1 Gp^a^ rub^ 
2mB 2 7r2 256 vr^ 



(3.23) 



with 



rpL/R 



d\ e-'^-y (b |t {ri^/^(y), r2^/^(o)} I b) , 



spin 



2 



P+aP-/3 



(3.24) 



(3.25) 



where Ti 



L/Rt 



L/R _ -pL/iJ 



2/1 



, given in eq. (3.14). The Dalitz distribution eq. (3.23) contains the ex- 



plicit 0(as) -improvement, and the distributions in which we are principally interested in can be obtained 
by straight-forward integrations. 

Using Lorentz decomposition, the tensor T^^, can be expanded in terms of three structure functions , 



-Ti g^v + T2 Vy + Ts ie^uais v° , 



(3.26) 



where the structure functions which do not contribute to the amplitude in the hmit of massless leptons 
have been neglected. After contracting the hadronic and leptonic tensors, one finds 



J.L/R j^L/Rl^'^ = m^^ 1 2 s T^^l^ + 



\2 -'- -2 

) --U -s 



(3.27) 



We remark here that the T3 term will contribute to the FB asymmetry but not to the branching ratio or the 
dilepton invariant mass spectrum in the decay B — > Xs£~^£~ ■ 



'We use the convention rr(7'"7''7"7'^75) = -4ie'"'"'', with e"^^^ = -1. 
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Normalization 

It has become customary to express the branching ratio for B Xsi~^£~ in terms of the well-measured 
semileptonic branching ratio Bgi for the decays B —^ {Xc, Xu)£u^ according to 



dB{B ^ Xsl^r) =Bs 



dr{B Xs 



'T{B ^{X„X^)£iye) ' 

This fixes the normalization constant Bq, which will be used throughout the following sections to be, 



(3.28) 



Bn = B 



si 



3a2 \V,lVtb\' 



1 



(3.29) 



Here f{rhc) is the phase space factor for r{B Xciv^) and the function K{mc) accounts for both 
the 0{as) QCD correction to the semileptonic decay width [ p2[ ] and the leading order (l/m^)^ power 
correction [p3[]. They read as: 



l-8mi + 8rh° 



2Amt \nrhr 



and 



where 



K{mc) 



asirrih) h{m. 
1 H giiric) + 



TT 



2ml ' 



g{mc) 
h{rhc) 



firric) 
Ai + 



f{mc 



-9 + 24m2 - 72mt + 72m^ - 15m^ - 72mt In 



rrir 



(3.30) 
(3.31) 

(3.32) 
(3.33) 



and the analytic form of AQ{rhc) can be seen in ref. [p3[]. Note that the frequently used approximation 
g{z) -|((7r^ - ~ ■^)^ + f) holds within 1.4% accuracy in the range 0.2 < z < 0.4. The equation 
g{z) = -1.671 + 2.04(z - 0.3) - 2.l5{z - O.S)^ is accurate for 0.2 < z < 0.4 to better than one per 
mille accuracy. 



3.1.3 Asymmetries in B ^ Xgi^i Decay 

Besides the differential branching ratio, B Xsi^i~ decay offers other distributions (with different 
combinations of Wilson coefficients) to be measured. An interesting quantity is the Forward-Backward 



(FB) asymmetry defined in | |13|J12| ] 

dAis) 



ds 



d^B 

ds dz 



dz 



^ 

-1 dsdz 



dz . 



(3.34) 



where z = cos is the angle of measured w.r.t. the 6-quark direction in the dilepton cm. system. 
From the experimental point of view, a more useful quantity is the normalized FB asymmetry, obtained 
by normahzing dA/ds with the dilepton mass distribution, dB/ds, 
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The asymmetry A, which we recall is defined in the dilepton c.m.s. frame, is identical to the energy 
asymmetry Ae introduced in as shown in ref. [llS]]. It is defined in the B rest frame as 



, _ iNiE^>E+)-NiE+>E.)) 

^ ~ {N{E^ > E+) + N{E+ > S-)) ■ ^ ' 

Here N{E^ > E^) denotes the number of lepton pairs whose negatively charged member is more ener- 
getic than its positive partner, where E± denote the £^ charged lepton energy in the B rest frame. The FB 
asymmetry is odd under charge conjugation in contrary to the differential branching ratio, which is charge 
conjugation even. Both observables contain non overlapping information which together can be used to 
test the SM. 



Another quantity is the left-right-asymmetry [ ]84| , |85| , ^q ], defined as 



(3.37) 



ds ds ds 

with dB^/ds (dB^/ds) denoting the invariant dilepton mass distribution for B Xs£^i~ decay into 
purely left-handed (right-handed) leptons. We can obtain dB^^^ /ds from the dilepton invariant mass 
distribution dB / ds by the replacements 

^eff Cf ^ =F C'lo ^ Cio =F Cg^ |^eff|2 l|^eff|2 iq\ 
^9 2 ' 2 ' 1^7 I ^ 2' ' ■ (-).38) 

Measurement of these asymmetries provides additional information on the underlying short-distance 
physics. 

3.1.4 Leading power {l/rrib) corrections in the decay B Xsi~^i~ 

We start with a discussion of the analyticity properties of the forward scattering amplitude T^^. They 
are determined by cuts, depending on the external states. We consider real particle production in the 
inclusive decay B —>■ Xsi^i^ , thus we have pb = Px + Q, where pb,Px denotes the 4-momentum of 
the S-meson, final hadronic state Xs, respectively. The hadronic invariant mass is in the range 

m\<p\<ml, (3.39) 

and the physical cut runs along the real axis in the complex v ■ q plane in the limits 



999 

Iq^ <vq< -^--1 i£ . (3.40) 

2mB 



The phase space integration, which follows the above cut, is over intermediate physical states and hence, 
depends on long-distance QCD. Moreover, at the upper bound of the cut, where ~ m^, resonances are 
dominating. In order to perform a reliable expansion in perturbative QCD, the contour of the integration 
has to be deformed in such a way that a) it encloses the cut and b) stays away from it by a distance large 
compared to Aqcd (see Fig. 1 in [ p^ ] for the contour of integration). The expansion is valid except in 
the corner of the Dalitz plot, where the hadronic invariant mass of the final state is small (the s-quark 



in Fig. 3.3 is almost on-shell.) However, results of perturbative QCD are expected to be recovered after 



suitable smearing. 
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Operator product expansion 



The next task is to expand the forward scattering amplitude T^p in eq. (3.26) in inverse powers of rri},. 



We employ HQE techniques which have been already sketched in section y.AA\ The leading term in this 
expansion, i.e., 0{m^) reproduces the parton model result. Let us describe how to get next to leading 
power corrections. First we write the momentum of the heavy 6-quark as pi,^ = mi^v^ + A;^, fix the four- 
velocity of the external 6-quark field to be and treat the components of the "residual momentum" of 
order Aqcd- We obtain the condition = ml + 2mhV.k + k'^, which yields p^ = ml + 0{k/mb). The 
heavy quark remains almost on-shell under soft gluon exchange with the light degrees of freedom, there 
is no anti-quark generated and the total 6-number is conserved. 
It is customary to define a field h with fixed velocity v through 



h{x) = e'""'''-''P+b{x) 



(3.41) 



with inversion 



b{x) 



1_ 

2mf) 



+ ... 



h{x) , 



(3.42) 



and the projection operators P± = (1 it y)/2. For the Dirac field b{x) the following identities hold 
P+b{x) = b{x) and P^b{x) = 0, which are corrected by terms of 0{l/mb)- Inserting this into the usual 
QCD Lagrangian £ = b{ilf) — mq)b we get the one in the HQET: 



C-HQET = hiv.Dh + S£ , 



(3.43) 



with 6C containing the corrections in l/mh'. 



6C 



1 



■h{iDfh 



1 



■Zi{ij)h{iv.Dfh + 



-^Z2{ij)h^a^'^G^,h + 0{l/ml) . (3.44) 



2mb 2mh 

Here, G^p = [iD^^iDy] denotes the gluon field strength tensor. For definition of the renormalization 



constants Zi_2(/u) we refer to [ |T7[ ] and references therein. For the sake of completeness we give the 



Feynman rules in the HQET in appendix |A.3.1| . The matrix elements of the above higher dimensional 
operators are given as 



B 



hiiDYh 



B 



2mB Ai , 
6mB A2 , 



(3.45) 



where B denotes the pseudoscalar i?-meson, see section |2.4.1| for a discussion of the values of the pa 



rameters Ai and A2. The second term on the r.h.s. of eq. ( |3.44[ ) vanishes by the lowest order equation of 
motion iv.Dh = 0. The on-shell condition of the heavy quark is ml = ml + 2mbV.k + k"^. Neglecting the 
last term, we have the simple condition v.k = 0. This is equivalent to the lowest order EOM of a heavy 
quark, thus we have the correspondence k ^ iD. 

'Note that here and only in this section we use this definition of the gluon field strength tensor G = G^^^, following the 
conventions of ref. |jl^. It is not consistent with the usual one appearing in QCD text books, eq. (A. 10 1, denoted here by G"^"^^. 
The two are related by G''^^ = igG^^'^. 
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Figure 3.3: The diagrams contributing to the operator product expansion. 

Suppressing the Lorentz indices for the time being, this operator product expansion (OPE) can be 
formally represented as: 



A'ye-'^-y {B\T{Tr{y)j2m\B) = - — 



+ 



{B\Oo\B) + -^{B\Oi\B) 
2 mi, 



{B\02\B) + 



4 nib' 



(3.46) 



The expressions for the operators Oq, Oi and O2 have been first derived in ref. [17], which we have 
checked and confirm. They are given as: 

1 



00 = -6ri(y-<f + m,)r26, 

X 

01 = -hTi^^T2iD^h-^{v-qThTi{i-i + ms)T2iD^h, 

X 



(3.47) 
(3.48) 



and 



Co 



(3.49) 



where 



16 4 

^{v- qTiv - qfhViii -i+ ms)T2iDaiD/3h - ^hTi{i - | + ms)T2{iDfh 
-^{v- qfhTi-t'"^2{iD^iDf3 + iDpiDa)h , (3.50) 
^m,/iriia„^r2G"^/i + ^ie^'^'^'iv - q)xhri-f^l5r2Gc.^h , (3.51) 
-/i(7^ri7°r2 + ri-f^r2^'')iDfsiD^h -^(v- qrh^r^i^ -i + ms)T2iDpiD^h 

X X'^ 

4 



X' 



{v - qThTi{i -i + rns)T2-i^iD^iDph 



(3.52) 



Here x = 1 + s — 2 {v ■ q) — ml + ie. The operator O3 responsible for \/m\ corrections can be seen in 



ref. [|65|]. 

The above Oj, i = 0, 1, 2 are obtained by expanding the Feynman diagrams shown in Fig. |3.3|, which 



contributes to the time-ordered product on the l.h.s. of eq. ( |3.46[ ). The diagram on the left is responsible 



for the operators Oo,0i,0^^ and To be definite, we write the intermediate s-quark propagator 

using 4-momentum conservation Ps = Pb — Q = fnbV — g + /c as 



1 



i - i+^/nib + rhs 



p1 — ml + ie mb x — 2q ■ k/mb + 2v ■ k/mb + k'^/m? 



(3.53) 
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and insert this into the diagrams Fig. 3.3. Expanding as well the propagator as the Dirac field b, which 
sandwiches the amplitudes, as the normalization of the states in powers of k/rrih, we obtain the desired 
OPE. Note that the leading operator Oq is defined in terms of the "full" four-component field b. The other 
two subleading operators Oi and O2 are, however, written in terms of the two-component effective fields 
h. In rewriting the operator Oi from b h fields by means of eq. (3.42), we obtain the operator 0^2^ 



Evaluation of the one-gluon diagram in Fig. 3.3 results in the operator O2 
The results of the power corrections to the structure functions Tj, i 



{9) 



1,2,3 can be decomposed 



into the sum of various terms, denoted by T- , which can be traced back to well defined pieces in the 
evaluation of the time-ordered product given above 



(3.54) 



j=0,l,2,s,g,5 



The expressions for t/-'^ {v.q, s) calculated up to 0{mB /mf) are given in appendix BJ_. They contain the 



parton model expressions Tf\v.q, s) and the power corrections in the HQE approach which depend on 



the two HQE-specific parameters Ai and A2 defined in eqs. ( |3.45| ). Note that the s-quark mass terms are 
explicitly kept in T^-'\v.q, s). The origin of the various terms in the expansion given in eq. ( 3.54 ) can be 
specified, as follows: 



T-^\v.q,s) 
Ti''\v.q,s) 
T^^'\v.q,s) 
T^^'\v.q,s) 



{B\Oo\B) for Ai 
{B\Oo\B 



A2 = 0, 



rPiv.qJ) 



{B\0^^'>\B) for j = 1,2,5 
{B\Oi\B) . 



(3.55) 



In the leading order in (1/771,5) the matrix element of Oi vanishes, but in the sub-leading order it receives 



a non-trivial contribution which can be calculated by using the equation of motion [17]. The contributions 



(s) ki 

arise from the matrix element of the scalar operator bb, i.e use of eq. ( [3.58D given below. We recall 
that the scalar current can be written in terms of the vector current plus higher dimensional operators as 



[15] 



1 



2mi 



bb = v^b-ff'b + {iDy - {v.iDY - (il2)a^''G^ 



fJ,U 



h + ... . 



With our normalization: 



{B\h^b\B) = 2{p_ 



it follows then 



{B\bb\B) = 2777^(1 + 



27772 



(Ai+3A2)) + 0(777b/7773) 



(3.56) 



(3.57) 



(3.58) 



Other possible Lorentz structures Uke 75, 757/^, 7/^7;^ sandwiched between b and b give zero after taking 
the i?-meson matrix element. 



.(i) 



From the expressions forT^-'^ given in appendix |B.1|, we see that T. 



.(0), 



1, 2, 3) are of order m,^ /777f, 
and the rest r/^\ Tl^\Tf'\T^^^ and t/^'' are all of order niBXi/mh^ or mB\2/'mb^ . Since the ratio 
niB/mb = 1 + 0{l/mb), we note that the Dalitz distribution in B ^ Xsi^l^ has linear corrections in 
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3.2 Power Corrections to the Dilepton Invariant Mass Distribution and FB 
Asymmetry 



The integration in the complex plane v ■ q can be done using the relation 

Im— oc ) \,, (^^"~^^ il-2vq + s-ml) . 
X" (n-1)! ^ ^ 



(3.59) 



Further, the integrand should be multiplied by the function 9{Av • — As"^ — v?) responsible for the 
correct integration boundary. ^ The resulting double differential branching ratio in i? ^ Xgi^i^ can be 
expressed as, 



dB 
ds du 



1 — rh\ 



-2 -2 

s — u 



i(2A,(- 



+3 A2 (- 1 + - 5mf - 8 s + 5 s^; 



-l + 2m2-m^-2s + s2) 



n&ii\2 I \ri |2 



+ 



4 1 



-2 - 4 - 6 

m„ — m„ + m„ 



|(2A,(- 



1 + rhi + fnt — m°„ + 2 s + lOm^s + + m^s 



o - 2 - -2 -2-2,-2,- 2 -2\ 
8m^s — s — nigS +u + nigU I 

^2-2\ 



+3 A2(3 + 5mi - 3m^ - 5m^ + 4 s + 28mp + 5 + Smp^' 



eff|2 



+ m^) - (1 - rhlf^ + + '^^1 - + s + mp) 



-A2(5m; 



Srrig + 2s + 5m^s] 



+2 
+4 



2 + Ai + 5 A2 1 n s Re{cf^) do 
2(1 + mfj + Ai(l + ml) + A2(3 + 5m^ u Re{do) cf^} 



Ei{s,u)6 u{s,ms)^ — u^ —E2{s,u)6' u{s,ms)'^ — ^ 



- n2 -2" 
u(s,ms) —u 

(3.60) 



where Ai = Ai/m^ and A2 = X2/1T1I. The auxiliary functions Ei{s,u) {i = 1,2), introduced here for 



ease of writing, are given explicitly in appendix B.2. The boundary of the Dalitz distribution is as usual 



determined by the argument of the 0-function and in the (n, s)-plane it has been specified in eq. (3.8) 



The analytic form of the result (3.6C) is very similar to the corresponding double differential distributions 



derived by Manohar and Wise in [ ]16[ ] for the semileptonic decays B {Xc, Xu)(.Vi- Further comparisons 
with this work in the V — A limit for the single differential and integrated rates are given a little later at 
the end of this section. 

Finally, after integrating over the variable u, we derive the differential branching ratio in the scaled 
dilepton invariant mass for B Xsl^i^ , 



2 1 

-n(s, ms){{l - ml)^ + s(l + ml) - 2s^) + -(1 - Ami + 6"i^ - 4m^ + 



+mls + rh'ts — m^„s — 38^ 



27hlf - Zmtf + 5s^ + bmls^ 



2s^ 



Ai 



u{s,rhs) 

+ {1 — 8ml + 18rfig — 16m^ + 5m^ — s — 3mp + Qm^s — 5m^s — ISs^ — ISm^s^ 



^This corresponds to q^ax = 'iE+E-, with lepton energies E± — v ■ q/2 ± u/{4mt). 
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{\c\ 



€ff|2 

9 I 



+ \C,o\' 



u{s,7hs) 

]u{s, ms)(2(l + m^)(l - m^)^ - (1 + 14m^ + - (1 + ^^s)'?^) 



+-(2 - 6ml + 4m^ + 4m^ - 6m^ + 2mJ° - 5s - 12m^s + 34m^s - 12m^s - 5mfs + 3s^ 



+20m^ - 12m^ - 14mf + Wmf + 3s + IGmls + 62m^s - 56m^s - 25mf s + 3s^ 



u{s, ifls 



eff 1 2 



+ 



8-^(5, ms)((l — m^)^ — (1 + + 4(1 — 2rTT,^ + — s — "fnls) u{s, rhs) \i 



+4 ( -5 + 30mf - 40m^ + 15mf - s + 21m^s + 25m^s - 45m^s + ISs^ + 22rhlf 



+45rfits^ — 7s'^ — Iftrhts 



2s3 



u{s,rhs) 

The leading power corrected expression for the FB-asymmetry A{s) is: 



(3.61) 



ds 



-2S, 



2(m(s, ms))^s + -(3 - 6rhl + 3m^ + 2s - 6m^s + 3s^)Ai 



+ s 
+ 



-9 - Qml + 15m^ - 14s - 30mp + ISs^) A2 
2 



iie(C9^ff) Cio 



4(n(s, rhs)f{l + m^) + - (1 + m^) (3 - 6m^ + 3rfif + 2s - 6m^s + 3s^)Ai 



+ 2(-7 - 3m2 - 5m^ + 15mf - 10s - 24m^s - 30m^s + 95^ + ISm^s^) A2I Re{Cio) C^^ ^ 



(3.62) 



The resuhs derived for the 0(qs) -improved and power-corrected Dalitz distribution, dilepton invariant 
mass, and FB-asymmetry in ^ Xsl'^l~ are the principal new results in this section. It is useful to write 
the corresponding expressions in the limit nis = 0. For the dilepton invariant mass distribution, we get 



2 Br, 



[I - s)2(l + 2s) (2 + Ai) + (1 - 15s2 + los^ A2 



{\C[ 



■eff|2 I \ri |2 
9 I + l^iol 



eff 1 2 



+ 



^(1 - s)2(2 + s) (2 + Ai) + 4 (-6 - 3s + Ss^) A2 
4(1 - s)2(2 + Ai) + 4 (-5 - 6s + if) A2] i?e(Cf^) C^^} . 



(3.63) 



The (unnormaUzed) FB asymmetry reads as. 



ds 



-2 Br, 



2(1 - sfs + -(3 + 2s + 3r)Ai + s (-9 - 14s + 15s^) A2 



Re{cf) Cio 



+ 



4(1 - sf + -(3 + 2s + 3s2)Ai + 2(-7 - 10s + 9f) A2 
3 



-eff 1 



(3.64) 



Our result [18] for the dilepton invariant mass distribution given in eq. (3.63) has been confirmed re- 



cently by [ETp in the = limit and is in disagreement with an earlier publication [ 17 1. (The differences 



between the previous result eq. (3.21) of the paper by [IM] have been discussed at length in [18].) 



Concerning the invariant dilepton mass spectrum derived by us and given in eq. (3.61), we would 
like to make the following observations: First, the leading order power corrections in the dilepton mass 
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distribution are found to be small over a good part of the dilepton mass s. However, we find that the 
power corrections become increasingly large and negative as one approaches s —>■ 5"^°^^, where s"^°'^ = 



(1 



Since the parton model spectrum falls steeply near the end-point s 



this leads to the 



uncomfortable result that the power corrected dilepton mass distribution becomes negative for the high 



dilepton masses. We show in Fig. 3.4 this distribution in the parton model and the HQE approach, using 



the central values of the parameters in Table A. 1 . Further, the power-corrected dilepton invariant mass 
distribution retains the characteristic 1/s behaviour following from the one-photon exchange in the parton 
model. We note that the correction proportional to the kinetic energy term Ai renormalizes the parton 
model invariant mass distribution multiplicatively by approximately the factor (1 + Ai/(2m^)), which is 
exact in the limit = and no new functional dependence in s is introduced (moreover, this factor 
is hardly different from 1). Hence, the negative probability near the end-point is largely driven by the 
magnetic moment term A2. 

The normalized FB asymmetry, dA{s)/ds, in the HQE-approach and the parton model are shown in 



Fig. 3.5. We find that this asymmetry is stable against leading order power corrections up to s < 0.6, 
but the corrections become increasingly large due to the unphysical behaviour of the HQE-based dilepton 
mass distribution as s approaches s™"^ (see Fig. 3.4). Based on these investigations, we must conclude 
that the HQE-based approach has a restrictive kinematical domain for its validity. In particular, it breaks 
down for the high dilepton invariant mass region in i? ^ Xs£~^l~ . This behaviour of the dilepton 



0.0000125 




Figure 3.4: Dilepton invariant mass spectrum dB{B Xse~^e ) /ds in the parton model (dashed curve) 
and with leading power corrections calculated in the HQE approach (solid curve). The parameters used 



are given in Table A. 1 



mass spectrum m B ^ Xs£~^£~ is not unexpected, as similar behaviours have been derived near the 
end-point of the lepton energy spectra in the decays B —>■ Xlv^ in the HQE approach [|^]. To stress 
these similarities, we show the power correction in the dilepton mass distribution as calculated in the HQE 
approach compared to the parton model through the ratio defined as: 



dB/ds{nQE) - de/(is(Parton Model) 



(3.65) 



dS/ds (Parton Model) 
The correction factor 

^HQE(-) foj. ^ ^ Xs£+e- shown in Fig. B]q is qualitatively similar to the 



corresponding factor in the lepton energy spectrum in the decay B — > Xc£iy£, given in Fig. 6 of [16]. We 
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Figure 3.5: FB asymmetry (normalized) dA{B Xge^e )/ds in the parton model (dashed curve) and 
with power corrections calculated in the HQE approach (solid curve). The parameters used are given in 



Table A.l 



note that we have been able to derive the power corrected rate for the semileptonic decays B —>■ Xciu^ 
obtained by Manohar and Wise in [p^. [| In doing this, we shall reduce the matrix element for the decay 
B Xs£^£~ to the one encountered in i? ^ Xcli'i, obtained by the replacements (V — A limit): 



a 



eff 



-C 



10 



a 



-eff 



(3.66) 
(3.67) 
(3.68) 



This amounts to keeping only the charged current V — A contribution in B ^ Xst^i" decays. 

Finally, since the HQE-improved expression for the decay rate cannot be given analytically due to the 
Wilson coefficient C|^(s) which is a complicated function of s, we give below the results in numerical 
form: 

pHQE ^ pfe^^ ^ ^^^^ ^ ^ (3_g9) 



where is the parton model decay width for b 



si^£ and the coefficients depend on the input 

0.501 



parameters. For the central values of the parameters given in Table A. 1 , they have the values ci 
and C2 = —7.425. This leads to a reduction in the decay width by —4.1%, using the values of Ai and A2 



given in Table A.l. Moreover, this reduction is mostly contributed by the A2-dependent term. We recall 
that the coefficient of the Ai term ci is (almost) the same as in the semileptonic width T{B Xui^e) 
obtained by Bigi et al [ p^ ] 



.HQE 



r 



si 



1 + ^Ai 



9 



Ac 



(3.70) 



where F^^ is the parton model decay width. This points towards the universality of this coefficient. The 
coefficient of the A2 term 02 for b — > si^i~ decay is larger than the corresponding one in the semileptonic 
decay width. Hence, the power corrections in T{B Xuii^e) and r{B Xgt^l^ ) are rather similar 
but not identical. 

2mh^Gb and 



-2mt^Kb, 3A2 



*The HQE matrix elements in our convention and the MW ones are related by Ai 

Al + 3A2 — —2mt^Eb — —2mt'^{Kb + Gt), likewise we have for the normalization of states \B) = ^/2mB\B)^''^ 
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Figure 3.6: The correction factor R^^^{s) (in percentage) as defined in eq. {3.65} for the dilepton mass 



spectrum dB{B — > Xsi~^i The parameters used are given in Table A.l 



3.3 Dilepton Invariant Mass and the FB Asymmetry in the Fermi Motion 
Model 



In this section, we present our estimates of the non-perturbative effects on the decay distributions in 
B —>■ Xsi^i^ ■ These effects aie connected with the bound state nature of the i?-hadron and the physical 
threshold in the B —>■ Xgi^i^ in the final state. In order to implement these effects on the decay distri- 
butions in B ^ Xsi^i^ , we resort to the Gaussian Fermi motion (FM) model [p2|] introduced in section 



2.4.2 



In the Fermi motion model, the problem of negative probabilities encountered in the HQE approach 
for the high dilepton masses near s —>■ Smax is not present, which motivates us to use this model as a 
reasonable approximation of the non-perturbative effects in the entire dilepton mass range. The success 
of this model in describing the inclusive lepton energy spectra in B —^ (X^, Xu)ii^e and B —>■ Xs"f 
strengthens this hope. 

In the decay B —>■ Xsl^i^ , the distribution dB/ds depends on the Lorentz-invariant variable s 
only. So, the Lorentz boost involved in the Fermi motion model (Doppler shift) leaves the dilepton mass 
distribution invariant. However, since the 6-quark mass niiy {p) is now a momentum-dependent quantity, 
this distribution is affected due to the difference {mi,{p) — m,},) (mass defect), which rescales the variable 
s and hence smears the dilepton distribution calculated in the parton model. For different choices of the 
model parameters {pp,mq) corresponding to the same effective 6-quark mass, m^^^ which is defined in 



eq. (2.47) the dilepton mass distributions should be very similar [23|, which indeed is the case as we have 
checked numerically but do not show the resulting distributions here. 

The situation with the FB asymmetry is, however, quite different. Being an angle-dependent quantity, 
it is not Lorentz-invariant and is sensitive to both the Doppler shift and the mass defect. We give in 
appendix B.3, the Dalitz distribution d^T{B Xsi^i~)/dsdu in the Fermi motion model. 

As we calculate the branching ratio for the inclusive decay B — > Xs£^£^ in terms of the semilep- 
tonic decay branching ratio B{B Xlv^), we have to correct the normalization due to the variable 



6-quark mass in both the decay rates. We recall that the decay widths for B ^ Xg 



and B 



Xiui 
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Figure 3.7: Differential branching ratio dB/ds for B Xgi^i^ (a) and normalized differential FB 
asymmetry dA{s)/ds (b) in the SM including the next-to-leading order QCD corrections. The dashed 
curve corresponds to the parton model with the parameters given in Table |4.i| and the solid curve results 
from the Fermi motion model with the model parameters {pF,nT-q) = (252, 300) MeV, yielding an effective 
b-quark mass rn'J^ = 4.85 GeV. 

in this model are proportional to {m^^)^ [^4[,|,53|. Hence the decay widths for both the decays in- 
dividually are rather sensitive to m^. This dependence largely (but not exactly) cancels out in the 
branching ratio B{B Xgi'^l'^ ). Thus, varying in the range = 4.8 ± 0.1 GeV results in 
Ar(i? Xgi'^i^ )/r = ±10.8%. However, the change in the branching ratio itself is rather modest, 
namely AB{B Xsl'^i^ )/B = ±2.3%. This is rather similar to what we have obtained in the HQE 
approach. 

The theoretical uncertainties in the branching ratios for B —>■ Xgl^i^ from the perturbative part, such 
as the ones from the indeterminacy in the top quark mass, the QCD scale Aqcd and the renormalization 
scale /i, have been investigated in the literature [pLPv]]. We have recalculated them for the indicated ranges 



of the parameters in Table A. 1 . The resulting (SD) branching ratios and their present uncertainties are 
found to be: 

B{B^Xse+e-) = (8.4 ± 1.9) x 10"^ , 
B{B ^ Xsp,~^fi~) = (5.7±0.9) X 10^^ , 

B{B X,r+r~) = (2.6 ± 0.4) x 10~^ , (3.71) 

where in calculating the branching ratio B{B —>■ XsT^t~), we have included the r-lepton mass terms in 
the matrix element [^]. These uncertainties, typically ±20%, are much larger than the wave function- 
dependent uncertainties, and so the theoretical accuracy of the SD-part in the SM in these decays is not 
compromised by the non-perturbative effects. 



We show the resulting dilepton invariant mass distribution in Fig. 3.7 (a) and the FB-asymmetry 
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(a) (b) 




Figure 3.8: Differential branching ratio dB/ds for B Xgi^i^ (a) and normalized differential FB 
asymmetry dA{s)/ds (b) using the Fermi motion model for three different pairs of the model parameters 
{pF,mq) = (450,0) MeV (solid curve), (310,300) MeV (long dashed curve), and {pF-,mq) = (310,0) 
MeV (short dashed curve) yielding the effective b-quark masses rn^^^ = 4.76 GeV, 4.80 GeV, and 4.92 
GeV, respectively. 



in Fig. 3/7 (b), where for the sake of illustration we have used the values {pp,mg) = (252,300) in 
(MeV,MeV), which correspond to an allowed set of parameters obtained from the analysis of the measured 
photon energy spectrum in — > Xs'y , using the same model [94|. We see that the dilepton mass distribu- 
tion is stable against Fermi motion effects over most part of this spectrum, as expected. We emphazise here 
that the end-point spectrum extends to the physical kinematic limit in i? ^ Xgi'^i" s™"^ = {mB — mxY 
with mx = max{mK, ms+mg) {niq is the spectator mass), which has to be imposed on the FM program. 
It corresponds to the invariant hadronic mass of the lowest physical state with total strangeness number 
s = 1, m{Xs) = rriK, as opposed to the parton model, in which s"^"-^ = (m;, — m^)^. The two thresholds 
can be made to coincide for only unrealistically values of rrif, and m^. The FB -asymmetry shows a more 
marked dependence on the model parameters, which becomes significant in the high dilepton mass region. 
As the parameters of the Fermi motion model are not presently very well-determined from the fits 



of the existing data [94,95|, one has to vary these parameters and estimate the resulting dispersion on 



the distributions m B ^ Xsi^l~ . We show in Figs. 3.8 the dilepton mass distribution (a) and the 
FB asymmetry (b), respectively, indicating also the ranges of the parameters {pp,mq). The resulting 
theoretical uncertainty in the distributions is found to be modest. 



3.4 LD Contributions in B ^ XJ+£- (I) 

Next, we implement the effects of LD contributions in the processes B Xgi'^i' . The issues involved 
here have been discussed in [96,28,?^]. The LD contributions due to the vector mesons J/t/j and i/j' 
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Figure 3.9: Differential branching ratio dB/ds for B Xgi^i^ (a) and normalized differential FB 
asymmetry dA{s) /ds (b) calculated in the SM using the next-to-leading order QCD corrections and Fermi 
motion effect (solid curve), and including the LD- contributions (dashed curve). The Fermi motion model 
parameters (pp, niq) in MeV are displayed in the figure. 



and higher resonances, as well as the (cc) continuum contribution, which we have already included in 
the coefficient Cf^^, appear in the (sL7^6L)(e7'^e) interaction term only, i.e., in the coefficient of the 
operator Og. This implies that such LD-contributions should change Cg effectively, but keep Cf^ and Cio 
unchanged. In principle, one has also a LD contribution in the effective coefficient Cf^^; this, however, has 
been discussed extensively in the context of the B — > Xg^ decay and estimated to be small [97,77,76.98]. 
The LD-contribution is negligible in Ciq. Hence, the three-coefficient fit of the data on B ^ Xgi^i" and 



B —>■ , proposed in ref. [|l2p on the basis of the SD-contributions, can be carried out also including 
the LD-effects. 

In accordance with this, to incorporate the LD-effects in i? — > Xgi^i^ , the function Y{s) introduced 
earlier is replaced by, 

Y{s) ^ Y'{s) = Y{s) + Yres{s) , (3.72) 
where Yres{s) accounts for the charmonium resonance contribution via B Xs{J/"^, . . .) ^ Xgi^i^ 



Its origin hes in the diagram displayed in Fig. 3.1, where the internal charm loop hadronizes before de- 



caying into a photon. We take the representation [13], 



E 



TTT(v,^e+e-)mv^ 



(3.73) 



l/i=V(ls),...,t/'{6s) 

where C(°) = 3Ci + C2 + 3C3 + C4 + 3C5 + Cq. We adopt k = 2.3 for the numerical calculations [||]. 
This is a fair representation of present data in the factorization approach [ ^9| ] ; also the phase of k, which is 
fixed in eq. (3.73), is now supported by data which finds it close to its perturbative value [ |10(]| ]. Of course, 
the data determines only the combination kC^^^ = 0.88. The relevant parameters of the charmonium 
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Figure 3.10: Dilepton invariant mass distribution in B Xgi^i^ (a) and normalized differential FB 
asymmetry dA{s)/ds (b) in the SM including next -to -leading order QCD correction and LD effects. The 
solid curve corresponds to the parton model and the short-dashed and long-dashed curves correspond to 
including the Fermi motion effects. The values of the FM parameters in MeV are indicated in the figure. 



resonances (15, . . . ,65) are given in the Particle Data Group [39], and we have averaged the leptonic 
widths for the decay modes V —>■ £^£^ for £ = e and £ = jj.. Note that in extrapolating the dilepton 
masses away from the resonance region, no extra g^-dependence is included in the 'y*{q'^)-Vi junction. 
(The (7^-dependence written explicitly in eq. (3.73) is due to the Breit-Wigner shape of the resonances.) 
This is an assumption and may lead to an underestimate of the LD-effects in the low-s region. However, 
as the present phenomenology is not equivocal on this issue, any other choice at this stage would have 
been on a similar footing. 



The resulting dilepton mass spectrum and the FB asymmetry are shown in Fig. |3^ (a) and Fig. ^ (b), 
respectively. We recall that the two curves labeled SD and SD+LD include explicit 0(as) -improvement, 
calculated in the parton model [PJ47[] and non-perturbative effects related with the bound state nature of the 
i?-hadrons and the physical threshold in the final state in B Xsi^i^ , using the Fermi motion model. 
In addition, the SD+LD case also includes the LD-effects due to the vector resonances, contributing 



to Cg^^ as discussed earlier. The parametric dependence due to the FM is shown in Figs. p.lO| for the 



dilepton mass spectrum (a) and the FB asymmetry (b), respectively, and compared with the case of the 
parton model in which case no wave function effects are included. These figures give a fair estimate of the 
kind of uncertainties present in these distributions from non-perturbative effects. In particular, we draw 
attention to the marked dependence of the FB asymmetry to both the LD-(resonances) and wave function 
effects, which is particularly noticeable in the region s > m^, . The dilepton invariant mass spectrum, on 
the other hand, is very stable except at the very end of the spectrum, which is clearly different in all three 
cases shown. This closes the first part of this chapter. 
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Bridge 

Concerning our approach to include resonant charm effects eq. ( 3.72 ), we compare it in section 3.10| 



with two other LD prescriptions, given in ref. [ ]28[ ] and | |29| ] and estimate the resulting uncertainties in the 
dilepton mass spectrum and the FB asymmetry. Further, we discuss a possible double counting, inherent in 
our procedure adding SD and LD amplitudes. The determination of these uncertainties is very important, 
as a measurement of the partly integrated spectra AB, AA in B —^ Xs£^i~ decay will be used to extract 
the SD coefficients, testing the SM. 

3.5 Introduction to Hadron Spectra and Spectral Moments in the Decay 

B XJ+i- 

In this second part of this chapter we present spectra in inclusive B — > Xsi~^£~ decay in kinematical vari- 
ables different from the dilepton invariant mass q'^, the hadronic energy and the hadronic invariant mass. 
Further, we calculate lowest moments in these hadronic variables. We include perturbative 0{as) correc- 
tions, 1 / mi) power corrections by means of the heavy quark expansion technique (HQET) and studies in 
the Fermi motion model (FM). This part is based on refs. [p5lj2^JP 



A similar program of investigations [93,23], [101-103] has been run for the charged current induced 



semileptonic B — > Xu^c^i^i decay. Here the main interest is focused on testing HQET and on the deter- 
mination of the CKM matrix elements Vcb and Vub- To be more specific, the HQE parameters Ai and A 
have been extracted from moments of the hadronic invariant mass spectrum in B ^ Xu^c^i^e decay [|6^. 
We recall that these non-perturbative parameters appear in the relation tub = mi, + A — (Ai + 3X2) /2mi, 



between the mass of the S-meson to the 5-quark mass (see section [2.4. 1[ ). Explicit calculation [ |25| , [26[ | 
shows that also in i? ^ Xs£'^i~ decay the hadronic invariant mass moments are sensitive to the HQET 
parameters Ai and A. This provides potentially an independent determination of these quantities. We 
think that the hadron spectra in i? — > Xgi^i^ and B Xuii'e can be related to each other over limited 
phase space and this could help in improving the present precision on Vub [ p^ ] and the parameters Ai and 
A [Hp. Of course, B Xu^c^i'e decays involve much less problems than FCNC B Xst^i decay, 
as the charged current mode has simpler short-distance (SD) couplings and no cc resonances present in the 
spectra. Besides these obvious differences, we will point out in the following sections similarities between 
rare B Xsi^i~ and the charged current B Xu^c^ve decays. 

What can we learn from the study of hadron spectra and moments in — > XsC-'^l~ ? Our motivation 
is manifold: 

• Hadron spectra have an interest on their own, they complete the profile of i? ^ Xsl'^£~ decay 
which has been given in the previous sections, i. e. the dilepton invariant mass distribution and the 
FB asymmetry. 



In their search for B Xgi £^ the CLEO collaboration [30] imposed a cut on the hadronic invari- 
ant mass Sh to suppress the BB background in measuring the dilepton invariant mass distribution. 
The hadronic invariant mass spectrum is absolutely necessary to acquire control over the signal after 
a cut in Sh- 
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• A possible determination of non-perturbative HQE parameters Ai, A from the first two moments of 
the hadronic invariant mass inB Xsi~^£~ decay, complementing the constraint from the charged 
current B Xclvt decay. The constraints from these decays can be used to reduce the present 
dispersion on Ai and A. 

• Test of the Fermi motion model in i? ^ Xsl^i~ decay. 

The power corrections presented here in the hadron spectrum and hadronic spectral moments in i? — > 
Xgl'^l'^ are the first results in this decay. 



3.5.1 Hadron kinematics 

Besides the parton level kinematics already introduced in section p.l.l[ , the corresponding kinematics at 
hadron level can be written as: 

B{pb) ^ Xs{pH)+t^{p+)+r{p^) . (3.74) 

The hadronic invariant mass is denoted by Sr = p\j and Eh denotes the hadron energy in the final 
state. The corresponding quantities at parton level are the invariant mass sq and the scaled parton energy 
xo = In parton model without gluon bremsstrahlung, this simplifies to sq = ml and xq becomes 
directly related to the dilepton invariant mass xq = 1/2(1 — s + m^). From momentum conservation the 
following equalities hold in the &-quark, equivalently i?-meson, rest frame {v = (1, 0, 0, 0)): 

xq = 1 — V ■ q , so = l — 2w-g + s, (3.75) 
Eh = niB — V ■ q , Sh = — 2mBV ■ q + s . (3.76) 

The relations between the kinematic variables of the parton model and the hadronic states , using the 
HQET mass relation, can be written as 

- A1 + 3A2 , f , Ai+3A2 ^ ^ 

Eh = A h - A H ] xq + . . . , 

ZniB V ZniB / 

Sh = m1 + A^ + (m| - 2Amij + A^ + Ai + 3A2) (so - rhi) 

+ (2AmB-2A2- Ai-3A2)xo + ... , (3.77) 

where the ellipses denote terms higher order in l/rrib. 



3.6 Perturbative 0{as) Corrected Hadron Spectra in B ^ Xst^l Decay 

In this section the 0{as) corrections to the hadron spectra are investigated. Following the argument given 
in section p.l.2[ , only O9 is subject to ag corrections and the corresponding Feynman diagrams can be 
seen in Fig. |3.2| . The effect of a finite s-quark mass on the 0{as) correction function is found to be very 
small. After showing this, we have neglected the s-quark mass in the numerical calculations of the 0{as) 
terms. 
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3.6.1 Hadron energy spectrum 



The explicit order Ug correction to Oq can be obtained by using the existing results in the literature as 
follows: The vector current Oq can be decomposed a.sV = {V — A) /2 + {V + A) /2. We recall that the 
{V — A) and (V + A) currents yield the same hadron energy spectrum [ |104| ] and there is no interference 
term present in this spectrum for massless leptons. So, the correction for the vector current case in B ^ 



Xs£~^l can be taken from the corresponding result for the charged {V — A) case 1 22,91], yielding 

Cgpixo) + Y{xo) (3.78) 



Cfixo) 



with 



p{x) 
a{x) 



1 + ^a{x) 
vr 



m'i 



(3.79) 
(3.80) 



(3x - 4x2 - 2m2 + Sm^x) 
where Y{xo) = Y{s) with s = 1 — 2xq + m^. The expression for Gi (x) with nis 7^ has been calculated 
in ^]\\. The effect of a finite rus is negligible in Gi(x), as can be seen in Fig. [3.11| , where this function is 
plotted both with a finite s-quark mass, mg = 0.2 GeV, and for the massless case, = 0. A numerical 
difference occurs at the lowest order end-point x™"^' = 1/2(1 + m^) (for mi = 0), where the function 
develops a singularity from above (xq > x™"^) and the position of which depends on the value of nis. 
The function Gi (x) for a massless s-quark is given and discussed below 



G^(x) 



Gi(x) 



+ 2(4x - 3) 



1 



TT 



84x^ + 585x2 



+ Li2(l-2x) 



- 1860x + 1215) + (8x - 9) ln(2x) 
} for < X < 1/2 , 



— (1 - x)(32x^ 
180^ 



136x^ + 1034x3 - 2946x2 + 1899x + 312) 



- — ln(2x - 
24 ^ 

+ x2(3-4x) 



l)(64x3 

V 
3 



48x2 - 24x - 5) 

4Li2(;^) +ln2(2x- 1) 
zx 



21n2(2x) 



for 1/2 < X < 1 . (3.81) 



The 0{as) correction has a double logarithmic (integrable) singularity for xq —>■ 1/2 from above 
(xo > 1/2). Further, the value of the order Og corrected Wilson coefficient Cf^^(xo) is reduced compared 
to its value with Ug = 0, therefore also the hadron energy spectrum is reduced after including the explicit 
order as QCD correction for < xq < 1/2. Note that the hadron energy spectrum for B —^ Xgi^i^ 
receives contributions for 1 > x > 1/2 only from the order Og bremsstrahlung corrections. 



3.6.2 Hadronic invariant mass spectrum 

We have calculated the order Ug perturbative QCD correction for the hadronic invariant mass in the range 
< sq < 1. Since the decay b —>■ s + £^ + i~ contributes in the parton model only at sq = ml, only 
the bremsstrahlung graphs h^s + g + l'^+l~ contribute in this range. This makes the calculation much 
simpler than in the full sq range including virtual gluon diagrams. We find 

^ = ^^o-^{^^^;^(93 - 41so - 955^ + hbsl) + \\n so(-3 - 5so + Osg - 2s^)}C| . (3.82) 
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Figure 3.11: The function Gi{x) is shown for rus = 0.2GeV (solid line) and for the massless case 



corresponding to eq. (3.81) (dashed line). 



Our result for the spectrum in B ^ Xs£~^i~ is in agreement with the corresponding result for the {V — A) 
current obtained for the decay B —>■ Xqiv^ in the mq = limit in [93] (their eq. (3.8)), once one takes 
into account the difference in the normalizations. We display the hadronic invariant mass distribution in 
Fig. 3.12 as a function of sq (with sq = m^SQ), where we also show the Sudakov improved spectrum, 
obtained from the 0{as) spectrum in which the double logarithms have been resummed. For the decay 
B — > Xuf-fi, this has been derived in [23], where all further details can be seen. We confirm eq. (17) 



of p3| ] for the Sudakov exponentiated double differential decay rate and use it after changing the 
normalization Fq Bq'^C^ for the decay B 
the kinematic variables (x, y) as 



Xsl I ■ The constant Bq is given in eq. (3.29). Defining 



V ■ q 



2 2 



the Sudakov-improved Dalitz distribution is given by 



dxdy 



[B^Xstr) 



4a^ ln(l - y) r 



x^f{l + 2x'^) expf-^ln 



3vr (1-2/) 



2a, 
37r 



{G{x) + H{y)) 



2as dH 
Svr dy 



(3.83) 



(3.84) 



(y) G 



where [ ]23[ ] 

G{x) = 

H{y) = 



[8x2(1 



2x^)lnx + 2(l 



(5 + 4x2) ln(l - x2) - (1 - x2)(5 + 9x^ 



'9 5 



6x 



2(l-x2)2(l + 2x2) 



+^2 + 2Li2(x2) -2Li2(l 



X 



(3.85) 



In 



2 - z{l - x) + K 



z 2 
x){3 + X + xz — z) [ 

(1 + X)2 



ln(l - z) - 2 In- 



z{l — x) + K 
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Figure 3.12: The differential branching ratio 



dB{B^Xse+e-) 

ds() 



in the parton model is shown in the 0{as) 



bremsstrahlung region. The dotted (solid) line corresponds to eq. ( p.82| ), (eq. f p.87| )j. The vertical line 
denotes the one particle pole from b si'^ 
as they tend to zero for larger values of sq. 



denotes the one particle pole from b s£^£ . We do not show the full spectra in the range < sq < 



2(l + x)2(l + 2x2) 



p7(l + x)(l + 2x2 
1-z 



+ (l-x)(3-2x2) 



(3.86) 



The quantity k in eq. ( p.86[ ) is defined as k = ^/z^JT^^x)'^~^^Axz. 

To get the hadronic invariant mass spectrum for a 6-quark decaying at rest we change variables from 
{x, y) to (g^, So) followed by an integration over q^, 



Arn^ 



1 



dxdy 2mfx{l — x)^ 



(3.87) 



The most significant effect of the bound state is the difference between and m^, which is domi- 
nated by A. Neglecting Ai, A2, i.e., using A = ms — rrib, the spectrum is obtained along the lines as 
given above for after changing variables from (x, y) to {q^ , Sh) and performing an integration over 
q^. It is valid in the region wtg "^^^"^^"^" < Sh < ni^B ''^B^ < Sh < m\, neglecting nis) which 
excludes the zeroth order and virtual gluon kinematics (sq = nil), yielding 



dB 

dSn 



4mf 



dq" 



d^B 



1 



dxdy 2m^mBx{\ — x)^ 



(3.88) 



The hadronic invariant mass spectrum thus found depends rather sensitively on m^, (or equivalently A), 
as can be seen from Fig. 3.13[ An analogous analysis for the charged current semileptonic B decays 
B Xuii^e has been performed in [101], with similar conclusions. 
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Figure 3.13: The differential branching ratio ^^^^'j^^ ^ ^^e hadronic invariant mass, Sh, shown 
for different values of nib in the range where only bremsstrahlung diagrams contribute. We do not show the 
result in the full kinematic range as the spectra tend monotonically to zero for larger values of Sh < m-'s- 

3.7 Power Corrected Hadron Spectra in B ^ Xgi'^i^ Decay 

We start directly with the structure functions Tj of the hadronic tensor, which are calculated up to O (m^ / ml ] 
in JTs]]. They have been decomposed into a sum of various terms, Ti{v.q, s) = J2j=o,i,2,s,g,s'^i''\''^-Q^ ^) 



(?) If-. 

where the individual parts can be seen in appendix |B.1| . After contracting the hadronic and leptonic 
tensors, one arrives at eq. (||27]). Now we are interested in a different set of kinematical variables. We 



transform {v • s) — > (xq, sq) with the help of the kinematic identities given in eq. ( [3.75[ ), and make the 
dependence on xq and sq explicit, 

1 .o 



--nib 



2(1 - 2x0 + so)Ti^/^ + 



so 



T2^/R^{l-2xo + so)uT3 



L/R 



(3.89) 



With this we are able to derive the double differential power corrected spectrum ^^^^^^ for the decay 
B —>■ Xsi^i^ ■ Integrating eq. (3.23) over u first, where the variable u is bounded by 



2Jxl- so<u< +2 J xl - So 



(3.90) 



we arrive at the following expression [ [26p 



dxo dso 
where 



-B^lm^xl - so I (1 - 2x0 + so)n{so, ocq) + ^^y^T2(so, xq) + 0(Aia,) , (3.91) 



i 8xo-4(| + A2) (|Cfff(s)p + |CioP 
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32(-2m^ - 2so - 4mgSo + xq + 5mfxo + sqXo + mpoXo) + ^^(^ + A2) 



X (-5 - ll?fi^ + 5so - "i^so + lOxo + 22m^xo - IOxq - lOm^Xo)^ 



eff 12 



(so - 2x0 + 1)^ 



+ 



-32 



-(m^ + so - XQ - mlxo) - 48(^ + A2) ) Re{Cl^{s)) cf^\ 



+ ^ I (^^(-2so - 3X0 + 5xg) + 8A2(-2so + xq + Sx^) j (\cf{s)\'' + jdop) 

— ^ — (6?Ti^ + 12so + ISm^so — 2sg — 2m3So — 3xo — 21m^xo — 13soXo — 19m^soXo 
— 3xo + 9m^Xo + SsqXo + Stti^soXq + 4xo + 4m^Xo) 

+ 32A2(— 2mg — 2mpo — 2sq — ^m^s^ + xq — m^xo — SsoXq — llmpoxo + Xg 

|(-<eff|2 

+ 13m^Xo + 5soXo + 5m^soXo) j ^ _ _^ ^^^2 

+ ( — %^{—2>m1 — 5so + 2mpo + 3xo + 6m^xo + 3soXo — Xq — Sm^Xg) 



\/2 2 2 2 22\\ -^^i^Q ('^)) ^-^7 

32A2(?7ig + So + 2m^so — xo + 2m^xo + 3soXo — Sxq — Sm^Xo) 



So -2x0 + 1 



1 \ ^2x [32X0 /|^eff/sM2 , 1^ |2 



-3Ai(so-x^)|^(|Cfi(s)|^ + |Cio 



"^^^/2 2 2 2\l I 

+ — — (— 2m^ — 2so — 4mgSo + xo + Sm^xo + sqXo + rngSQXo)-p- 



^ -128, .2^-^ .2 , i^e(Cgff(s))Cfff ] 

+ -^K + so-xo-m,xo) ^^_2^^ + i [ , 

2 



r2(so,xo) = i|(^16-40(^ + A2)^ (|C7eff(J)|2 + |Cio| 



^1 , \ \ Wi , ^2^ IC7 



eff 1 2 



+ -64 + 160(^ + A2) (1 + m 



3 so -2x0 + 1 



101^ 



(448 A \ |/^eff|2 

_ + 64A2(5xo - 1) (1 + m^) J ^ ' , - 64A2i?e(Cf"(s)) 
3 / So - 2x0 + 1 



eff/'s^^ /^eff I 



Here, x = so — rhi +ie, Xi = Ai/ ml and A2 = A2/ ?n.^ . As the structure function Ts does not contribute to 
the branching ratio, we did not consider it in the calculation of the hadron spectra. The Wilson coefficient 
Cg^^{s) depends both on the variables xo and so arising from the matrix element of the four-Fermi- 
operators. Here the normalization constant Bq, defined in eq. ( 3.29| ), expresses the branching ratio for B —>■ 



Xst^d- as usual in terms of the semileptonic decays B Xclu^. The double differential ratio given in 



eq. (3.91) agrees in the (V — A) limit given in eqs. (3.66) - (3.68) with the corresponding expression 
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derived for B —>■ Xclvp^ decay in [|93|] (their eq. (3.2)). 



The hadron energy spectrum can now be obtained by integrating over sq. Using eq. (3.59), the follow- 
ing replacements are equivalent to taking the imaginary part 

1 



X 

1 



The kinematic boundaries are given as: 



5{sq 
-5' (so - 



ml) 



ml) 



(3.93) 



maxirh , —1 + 2xo + 4m, ) < sq 



rhs < xq 



4mf 



(3.94) 



Here we keep rhi as a regulator wherever it is necessary and abbreviate Cf = Cq^^{s = 1 — 2xo + ml 



Including the leading power corrections, the hadron energy spectrum in the decay B 
below: 



dxQ 



+ 



+ 



„(9,10) 



[9o + Mgi + A252 



{\cf\' + \C^o\' 



eff |2 



(9) 



Xq 

eff 1 2 



i(l + m2^ 



+ 



„(7,9) , I (7,9) , I (7,9) 
% + Al9l + ^252 



dsQ 



, , j2|^eff|2 



(Ai/if'')+A2/i^' 
1 



(7,9). 



di?e(C7fff) ^eff , ^ J7,9) d^i?e(C|ff )^eff 



dsn 



dsl 



Xsl'^l is given 



eff N r^eff 



+ <5(xo - -(1 + - 4m/))/5(Ai, A2) + b\xQ - -(1 + - 4m/))/5,(Ai, A2) 



(3.95) 



The functions gf'"^^^ , g"^^ , gf^'^'* , /ij-^^ , ^i'^'^ , , k^'^^ in the above expression are the coefficients of the 
I /ml power expansion for different combinations of Wilson coefficients, with Qq'''^ being the lowest 
order (parton model) functions. They are functions of the variables xq and rhs and are given in appendix 
|C.1| . The singular functions 6, 6' have support only at the lowest order end-point of the spectrum, i.e., 
at x™"^ = ^(1 + rrig — 4mf). The auxiliary functions /5(Ai,A2) and /5'(Ai,A2) vanish in the limit 

Ai = A2 = 0. They are given in appendix |C.2[ 



The derivatives of Cg^ are defined as ^^-^ 



eff 



d"Cl 



eff 



-is = l 



-gjfi — i^c, — J. — 2xo + So; So = "T-s) ("^ = 1) 2). In the {V — A) limit our eq. ( 3.95D for the hadron energy 
spectrum in B ^ Xs£~^£~ agrees with the corresponding spectrum in B 
eq. (Al)). Integrating also over xq the resulting total width for B Xgi^ 



Xlvn given in ref. [93] (their 
" agrees again in the {V — A) 



limit with the well known result [15]. 



The power-corrected hadron energy spectrum '^^'^^~'^f ^ ^ (with Eq 



mhXQ) is displayed in 

Fig. 3.14 through the solid curve, however, without the singular 6,6' terms. Note that before reach- 
ing the kinematic lower end-point, the power-corrected spectrum becomes negative, as a result of the 
A2 term. This behavior is analogous to what has already been reported for the dilepton mass spectrum 
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Figure 3.14: Hadron energy spectrum "^^^^^^f ^ ^ in the parton model (dotted line) and including 
leading power corrections (solid line). For mb/2 < Eq < mi, the distributions coincide. The parameters 



used for this plot are the central values given in Table A.l and the default values of the HQET parameters 
specified in text. 

'^^''^^d^i^^^ ^ region []T8[], signaling a breakdown of the ^ expansion in this region. The 



terms with the derivatives of Cg^^ in eq. (3.95) give rise to a singularity in the hadron energy spectrum at 



the charm threshold due to the cusp in the function Y{s), when approached from either side. The hadron 
energy spectrum for the parton model is also shown in Fig. p. 14 , which is finite for all ranges of Eq. 



What is the region of validity of the hadron energy spectrum derived in HQET? It is known that in 
B —^ Xgi'^l" decay there are resonances present, from which the known six [39] populate the xq (or 



Eq) range between the lower end-point and the charm threshold. Taking this into account and what has 
been remarked earlier, one concludes that the HQET spectrum cannot be used near the resonances, near 
the charm threshold and around the lower endpoint. Excluding these regions, the spectrum calculated in 
HQET is close to the partonic perturbative spectrum as the power corrections are shown to be small. The 
authors of [^], who have performed an l/rric expansion for the dilepton mass spectrum '^^^'^"^^"^ ^ ^ 
and who also found a charm-threshold singularity, expect a reliable prediction of the spectrum for < 
3m^ corresponding to £^0 ^ -^(1 + "^1 — 3m,^) w 1.8 GeV. In this region, the effect of the l/m^ power 
corrections on the energy spectrum is small and various spectra in i? — > Xsi'^i^ calculated here and in 



[18] can be compared with data. 



The leading power corrections to the invariant mass spectrum is found by integrating eq. (3.91) with 
respect to xq. We have already discussed the non-trivial hadronic invariant mass spectrum which results 
from the 0{as) bremsstrahlung and its Sudakov-improved version. Since we have consistently dropped 



everywhere terms of 0{\ias) (see eq. (3.91)), this is the only contribution to the invariant mass spectrum 
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also in HQET away from sq = ml, as the result of integrating the terms involving power corrections 
in eq. (3.91) over xq is a singular function with support only at sq = ml. Of course, these corrections 
contribute to the normalization (i.e., branching ratio) but leave the perturbative spectrum intact for sq ^ 



mt 



3.8 Hadronic Moments in B ^ XJ+£- in HQET 

We start with the derivation of the lowest spectral moments in the decay B —^ Xgl^t^ at the parton level. 
These moments are worked out by taking into account the leading power l/mt and the perturbative 0{as) 
corrections. To that end, we define: 



M 



(n,m) 

e+e- 



2\n'm 







dsodxo 



oux-Q , 



(3.96) 



for integers n and m. These moments are related to the corresponding moments (x^(so — "^i)") obtained 
at the parton level by a scaling factor which yields the corrected branching ratio B = BqM!^^^ . Thus, 

(3.97) 



{x^{so-mlr) = ^M^;^T^ 



The correction factor Bq/B i?> given below in eq. ( 3.103 ). We remind that one has to Taylor expand it in 
terms of the 0{as) and power corrections. The moments can be expressed as double expansion in 0{as) 
and l/mft and to the accuracy of our calculations can be represented in the following form: 



M 



{n,m) 

e+e- 



vr 



with a further decomposition into pieces from different Wilson coefficients for i = 0, 1, 2: 



D. 



{n,m) 



a. 



(n,m) ^eff 



7 



+ pt'"'^C!o + 7r"'''C^'' + 5, 



(n,m) ^eff I An,m) 



(3.98) 



(3.99) 



The terms 7^ '""^ and 5^ '""^ in eq. ( ^3% result from the terms proportional to Re{Cl^)C^^ and \Cl^\^ 
in eq. ( |3.91| ), respectively. The results for q,("'™)^^|"'™)^-)/|"'"^)^(^("'™) are presented in appendix |C.3[ 
Out of these, the functions a^-"'™^ and are given analytically, but the other two and S^"^'^^ 

are given in terms of a one-dimensional integral over xq, as these latter functions involve the coefficient 
C|^, which is a complicated function of xq. 

The leading perturbative contributions for the hadronic invariant mass and hadron energy moments 



can be obtained analytically by integrating eq. ( [3.82| ) and eq. ( |3.81[ ), respectively, yielding 

^(0,0) ^ ^5 ~ ^(1,0) ^ ^(2,0) ^ _±_ 

9 ' 675 ' 486 ' 

^(0,1) ^ 1381 - 21071^ ^(0 2) ^ 2257 - 320^^ 



(3.100) 



The zeroth moment n = 
derived some time ago 



1350 ' 5400 

m = is needed for the normalization and we recall that the result for A^^'^^ was 
Likewise, the first mixed moment A^^'^^ can be extracted from the results 



(3.101) 



given in ||93|] for the decay B X£i>£ after changing the normalization, 

^(1,1) = A. 
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For the lowest order parton model contribution Dq ' , we find, in agreement with []93D, that the first two 



hadronic invariant mass moments (sq 
vanish: 



mz 



{{so 



m 



2\2\ 



and the first mixed moment {xq{so — m^)) 



^(n,o) ^ for n 



1,2 and D, 



(1,1) 





0. 



(3.102) 



We remark that we have included the s-quark mass dependence in the leading term and in the power cor- 
rections, but omitted it throughout our work in the calculation of the explicit Ug term. All the expressions 
derived here for the moments agree in the V — A limit (and with rfis = in the perturbative correction 
term) with the corresponding expressions given in ref. [^]. From here the full 0{asms) expressions can 
be inferred after adjusting the normalization Fq —>■ Bo^Cg. We have checked that a finite s-quark mass 



effects the values of the given in eqs. ( |3lOOH3?ToTl ) by less than 8% for = 0.2 GeV. 

We can eliminate the hidden dependence on the non-perturbative parameters resulting from the b- 
quark mass in the moments M.^+^J with the help of the HQET mass relation. As rus is of order Aqcd, 
to be consistent we keep only terms up to order m^/m^ [66]. An additional mb-dependence is in the mass 
ratios rhi = Substituting by the i?-meson mass using the HQET relation introduces additional 
©(l/m^, l/m"^) terms in the Taylor expansion of eq. (3.97). We get for the following normalization 
factor for e/eo = A^l°'°^ ■ 



+ 



+ 



+ 



+ 



+ 
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+ / ^ dxo(64(-l - 4x0 + 7xi)Re{cf^)C^^ + 16(-1 + 15x^ - 20xg)|C| 



-Yeff|2i 



(3.103) 



ml 



Here, the ^ and ^ terms proportional to result from the expansion of 



ln( 



4mf , 
ml ' 



ln( 



4mr 



+ 2- 



m 



B 



K_ 

rriB 



+ 



A2 Ai + 3A2 



+ ... 



m 



B 



m. 



(3.104) 



The first two moments and the first mixed moment, {xq)B/Bq, {xq)B/Bo, {sq 
ml)'^)B/BQ and (xo(so — m^))B/Bo are presented in appendix |C4 . 



7nl)B/Bo, {{so 



With this we obtain the moments for the physical quantities valid up to 0{as/m'^, 1/m^), where 
the second equation corresponds to a further use of m^ = 0{Aqcd)- We get for the first two hadronic 
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invariant mass moments ^ 

{Sh) = m2 + A2 + (m|-2AmB)(so-m^) + (2Amij-2A2-Ai-3A2)(xo), 

= + 2A2m^ + 2m^(m|-2AmB)(so-m^) + 2m^(2AmB-2A2-Ai-3A2)(xo) 

+ {m% - 4Am|)((so - m^sf) + 4A2m|(xg) + 4Am|(a;o(so - m'i)) , (3.105) 
= {m% - 4Am|)((so - m^sf) + AA^mlixl) + 4Am|(xo(so - m^)) , 

and for the hadron energy moments: 

= A2 + (2AmB-2A2-Ai-3A2)(xo) (3.106) 
+(m| - 2Ami3 + A^ + Ai + 3A2)(a;g) . 



3.8.1 Numerical estimates of the hadronic moments in HQET 

Using the expressions for the HQET moments given in appendix |C.4| , we present the numerical results for 
the hadronic moments in 5 — > Xs£^£~ , valid up to 0{as/m'^, 1/m^). We find: 
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ml 

+ 0.641^ + 0.589 



A9 
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B 



m 



2 ' 
B 



{{so 



m 



2\2\ 



0.0071 — (1 + 0.083—) - 0.196^ 

vr mB mi 



(3.107) 



As already discussed earlier, the normalizing factor B/Bq is also expanded in a Taylor series. Thus, in 
deriving the above results, we have used 

J- = 25.277 (1 - 1.108— -0.083— - 0.041^ + 0.546^ - 3.439^). 

Bq TT mB m^^ m^ 



The parameters used in arriving at the numerical coefficients are given in Table A.l and Table 3.1 



Inserting the expressions for the moments calculated at the partonic level into eq. ( |3. 105[ ) and eq. ( |3.106[ ), 
we find the following expressions for the short-distance hadronic moments, valid up to ©(a^/m^, l/m%): 

-.2 A ^, A ^2 



{Sh) 



mz 



m 



+ 0.093- 



m 



B 



vr 



0.069—— + 0.735— + 0.243 

mB vr mB 



+0.273^ 



0.513- 



m 



B 



Our first expression for (Su), eq. (3.105), does not agree in the coefficient of {sq — rh^) with the one given in ref. [ |93| ] (their 
eq. (4.1)). We point out that should have been replaced by in this expression. This has been confirmed by Adam Falk 
(private communication). Dropping the higher order terms given in their expressions, the hadronic moments in HQET derived 
here and in \9M agree. 
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(S^) = m|(0.0071^ + 0.138 

vr mB TT 



+ COST- 



AL 



0.196- 



ml 



= 0.367mij(l + 0.148 



A A 
— -0.352 + 1.691 — 



{E'jj) = 0.147m|(l + 0.324— - 0.128 



vr 
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niB TT TUB mj^ rriB rn 



(3.108) 

A2 . 



A as A 
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niB TT niB 



2.740- 
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A2 . 

2- + 0.162^ 



m 



B 



One sees that there are linear power corrections, 0{K/mB), present in all these hadronic quantities except 
(Sir) which starts in 

Setting mg = changes the numerical value of the coefficients in the expansion given above (in 
which we already neglected a^ms) by at most 1%. With the help of the expressions given above, we 
have calculated numerically the hadronic moments in HQET for the decay B —^ Xs£^i~, i = ijL,e and 



have estimated the errors by varying the parameters within their ±lo" ranges given in Table A.l. They 



are presented in Table |3.2| where we have used A = 0.39 GeV, Ai = -0.2 GeV^ and A2 = 0.12 GeV+ 
Further, using ^^(m;,) = 0.21, the explicit dependence of the hadronic moments given in eq. (3.108) on 



the HQET parameters Ai and A can be worked out [25| 
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{El) 
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TUB 
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niB 
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m. 



(3.109) 



A A^ A 

0.372771^(1 + 1.64 + 0.01^ + 0.02^) 

rrifi rrin rrir. 



TUB 

A 



0.150m|j(l + 2.88 + 2.68 



nil 
A2 



B 



0.29- 



Ai 



rrir 



rriB mg no^ 

While interpreting these numbers, one should bear in mind that there are two comparable expansion pa- 
rameters K/niB and Os/vr and we have fixed the latter in showing the numbers. As expected, the depen- 
dence of the energy moments {E^) on A and Ai is very weak. The correlations on the HQET parameters 
Ai and A which follow from (assumed) fixed values of the hadronic invariant mass moments {Sr) and 
{SD are shown in Fig. 3.15. We have taken the values for the decay B X^/i^^" from Table for 



the sake of illustration and have also shown the presently irreducible theoretical errors on these moments 



following from the input parameters mt, ag and the scale /x, given in Table A.l. The errors were calcu- 
lated by varying these parameters in the indicated range, one at a time, and adding the individual errors in 
quadrature. Further the correlation following from the analysis of data on semileptonic B Xivi decays 
[ pTl ] is shown in Fig. 3.15 (ellipse). As can be seen, it gives a complementary constraint to the one from 
B Xgi^i^ decay [25 1, which allows in principle a precise determination of A, Ai from data on the 
latter. 

The theoretical stability of the moments has to be checked against higher order corrections and the 
error estimates presented here will have to be improved. The "BLM-enhanced" two-loop corrections 
[105] proportional to a^Po, where Pq = 11 — 2nj/3 is the first coefficient in the QCD beta function. 



can be included at the parton level as has been done in other decays [ p3|j68[ ], but not being crucial to our 
point we have not done this. More importantly, higher order corrections in Us and 1 / are not included 
here. While we do not think that the higher orders in will have a significant influence, the second 
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moment {Sfj) is susceptible to the presence of 1/m^ corrections as shown for the decay B Xlui 
[106]. This will considerably enlarge the theoretical error represented by the dashed band for {Sjj) in 
Fig. 3.15| . Fortunately, the coefficient of the A/ms term in (Sh) is large. Hence, a measurement of this 
moment alone constrains A effectively. Of course, the utility of the hadronic moments calculated above 
is only in conjunction with the experimental cuts. Since the optimal experimental cuts in B —^ Xsl^l~ 
remain to be defined, we hope to return to this and related issue of doing an improved theoretical error 
estimate in a future publication. We remark here that care has to be taken in a general HQE calculation 
with cuts. For an extraction of meaningful observables the calculation must be smeared by integration. If 
the remaining phase space gets too restricted the OPE, which is based on parton-hadron duality, breaks 
down. This happens for example near the high-(7^ end-point of the invariant dilepton mass spectrum in 



B Xsi+i- decay [jT8|]. 

Related issues in other decays have been studied in literature. The classification of the operators 
contributing in 0{l/ml), estimates of their matrix elements, and effects on the decay rates and spectra in 
the decays B Xiui and B ^ {D, D*)£i'£ have been studied in [ ]107|jl08|JT09| ]. Spectral moments of the 
photon energy in the decay B Xg^ have been studied in [110]. For studies of 0{l/ml) contributions 
in this decay and the effects of the experimental cut (on the photon energy) on the photon energy moments. 



see [65]. An HQE analysis of the first two hadronic invariant mass moments with a lepton energy cut has 
been worked out in ref. [11061]. 




Figure 3.15: {Sh) (solid bands) and (Sff) (dashed bands) correlation in fAi-Aj space for fixed values 
(Sh) = 1-64 GeV^ and (Sfj) = 4.48 GeV^, corresponding to the central values in Table [jT^ The curves 



are forced to meet at the point Ai = —0.2 GeV^ and A = 0.39 GeV. The correlation from B X£v£ [3]_] 
is also shown here (ellipse). 

Finally, concerning the power corrections related to the cc loop mB ^ Xst^(-~ , it has been suggested 
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HQET 


(Sh) 




(Eh) 


(Eh) 




(GeV^) 




(GeV) 


(GeV^) 


e+e" 


1.64 ±0.06 
1.79 ±0.07 


4.48 ±0.29 
4.98 ±0.29 


2.21 ±0.04 
2.41 ± 0.06 


5.14 ±0.16 
6.09 ±0.29 



Table 3.2: Hadronic spectral moments for B — > X^/U^/x and B Xge^e in HQET with A = 
0.39 GeV, Ai = —0.2 GeV'^, and A2 = 0.12 GeV'^. The quoted errors result from varying /i, Ug and the 



top mass within the ranges given in Table 4. 1 



in [|^] that an 0{K^(j]j/'m^) expansion in the context of HQET can be carried out to take into account 
such effects in the invariant mass spectrum away from the resonances. Using the expressions (obtained 
with nis = 0) for the 1/m^ amplitude, we have calculated the partonic energy moments A(a;Q), which 
correct the short-distance result at order A2 / ml : 
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(3.110) 



(3.111) 



2Vr-(r- 1) \ 1 ± VI - 1/^ 

The invariant mass and mixed moments give zero contribution in the order we are working, with = 0. 
Thus, the correction to the hadronic mass moments are vanishing, if we further neglect terms proportional 
to and ^Aj, with i = 1,2. For the hadron energy moments we obtain numerically 



HEh)i/. 
A {El 



mijA(xo) = -0.007 GeV , 
m|A(x^) = -0.013 GeV^ , 



(3.112) 



leading to a correction of order —0.3% to the short-distance values presented in Table 3.4 



3.9 Hadron Spectra in the Fermi Motion Model 



In this section, we study the non-perturbative effects associated with the bound state nature of the B- 
hadron on the hadronic invariant mass and hadron energy distributions in the decay B X^i'^t' . These 



effects are studied in the Fermi motion model (FM) [22] introduced in section 2.4.2. In the context of rare 
B decays, this model has been employed to calculate the energy spectra in the decay i? — > ± 7 in 



,24], which was used subsequently by the CLEO collaboration in their successful search of this decay 



[11]. It has also been used in calculating the dilepton invariant mass spectrum and FB asymmetry in 
B Xs£+i- lITSl], see section |33. 
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Comparison with HQET 



The FM has received a lot of phenomenological attention in B decays, partly boosted by studies in the 
context of HQET showing that this model can be made to mimic the effects associated with the HQET 



parameters A and Ai [ |111| , |16| ]. We can further quantify this correspondence. The HQET parameters are 
calculable in terms of the FM parameters pp and rUg with 



A 
Ai 



dp (j){p) m"^ 
dpp'^cl){p) = - 



-P 
3 



(3.113) 



10 r 

In addition, for niq = 0, one can show that A = 2ppl There is, however, no parameter in the FM 
model analogous to A2 in HQET. Curiously, much of the HQET malaise in describing the spectra in the 



end-point regions is related to A2, as also shown in 1 16,18]. 



The relation between m^, ra^. A, Ai and A2 in HQET has already been stated (eq. (2.41)). With the 
quantity rr^ defined in eq. (2.47) and the relations in eqs. (3.113) for Ai and A, the relation 



+ A 



Ai/(2mf) 



TUB = mj;" + Ai/[-zmj;") , (3.114) 
is found to be satisfied in the FM model to a high accuracy (better than 0.7%), which is shown in Table 



3.3| for some representative values of the HQET parameters and their FM model equivalents. We shall use 
the HQET parameters A and Ai to characterize also the FM model parameters, with the relations given in 



eqs. ( I3.113D and ( |2.47| ) and in Table ^ 



PF,rnq (MeV,MeV) 


mf^ (GeV) 


Ai (GeV2) 


A (GeV) 


(450, 0) 


4.76 


-0.304 


0.507 


(252, 300) 


4.85 


-0.095 


0.422 


(310,0) 


4.92 


-0.144 


0.350 


(450, 150) 


4.73 


-0.304 


0.534 


(500, 150) 


4.68 


-0.375 


0.588 


(570, 150) 


4.60 


-0.487 


0.664 



eff 

Table 3.3: Values of non perturbative parameters m, , Ai and for different sets of the FM model 



parameters (pp, ruq) taken from various fits of the data on B ^ Xg + 7) decays discussed in [14 ]. 



Calculation of the hadron spectra 

We turn to discuss the hadron energy spectrum in the decay B —>■ Xgi^i^ in the FM model including the 
0{as) QCD corrections. The spectrum -^^[B Xsl'^i^) is composed of a Sudakov improved piece 
from Cq and the remaining lowest order contribution. The latter is based on the parton model distribution, 
which is well known and given below for the sake of completeness: 

^ = Bo^!^^{mt-2mlml + mt + mls + m^,s-2s^){\cf{s)\'' + \Cio\^) 
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16(m^ - 2m2m^ + - mfs - mls)Re{cf^{s)) C7^^} 
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ml , 



s 

(3.115) 



(3.116) 



where u is given in eq. (^). Note that in the lowest order expression just given, we have |Cf^^(s)p = 
+ 2CgRe{Y{s)) with the rest of Cf^(s) now included in the Sudakov-improved piece as can be 
seen in eq. ( [3.84| ). To be consistent, the total semileptonic width T^/, which enters via the normalization 
constant Bq, has also to be calculated in the FM model with the same set of the model parameters. We 
implement the correction in the decay width by replacing the 6-quark mass in Tgi given in eq. ( 3.116| ) by 
[18]. The hadronic invariant mass spectrum in the decay B Xgi^i^ in this model is calculated 



m 



eff 



very much along the same lines. The kinematically allowed ranges for the distributions are mx < Eh < 
niB and < Sh < rri^, and we recall here that the physical threshold has been implemented by 
demanding that the lowest hadronic invariant mass produced in the decay B —>■ Xgi^i^ satisfies mx = 
max{mK-,rnq + ms). The results for the hadron energy and the hadronic invariant mass spectra are 
presented in Figs. 3.16| and 3.18 , respectively. We do not show the Sh distribution in the entire range, as 
it tends monotonically to zero for larger values of Sh- 




Eh [GeV] 



Figure 3.16: Hadron energy spectrum in B Xgi'^i in the Fermi motion model based on the per- 
turbative contribution only. The solid, dotted, dashed curve corresponds to the parameters (Ai, A) = 
(-0.3,0.5), (-0.1,0.4), (-0.15,0.35) in (GeV^, GeV), respectively. 

A number of remarks is in order: 
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Figure 3.17: Hadron energy spectrum in B Xgi'^i based on the perturbative contribution only, in 

eff 



the Fermi motion model (dotted curve) for (ppjiriq) = (252, 300) {MeV, MeV), yielding m^^ 
GeV, and in the parton model (long-short dashed curve) for nifj = 4.85 GeV. 



4.85 



The hadron energy spectrum in B Xsi^£ is rather insensitive to the model parameters. Also, 
the difference between the spectra in the FM and the parton model is rather small as can be seen 



in Fig. 3.17. Since, away from the lower end-point and the cc threshold, the parton model and 
HQET have very similar spectra (see Fig. 3.14), the estimates presented in Fig. 3.16| provide a 
good phenomenological profile of this spectrum for the short-distance contribution. Very similar 
conclusions were drawn in [23] for the corresponding spectrum in the decay B Xu£v£, where, 
of course, the added complication of the cc threshold is not present. 



In contrast to the hadron energy spectrum, the hadronic invariant mass spectrum m B ^ X, 



is sensitive to the model parameters, as can be seen in Fig. 3.18. Again, one sees a close parallel 
in the hadronic invariant mass spectra m B ^ Xgl'^l" and B X^ii^i, with the latter worked 



out in [101]. We think that the present theoretical dispersion on the hadron spectra in the decay 



B Xgl'^i can be considerably reduced by the analysis of data in B ^ Xy 



• The hadronic invariant mass distribution obtained by the 0(0^) -corrected partonic spectrum and 
the HQET mass relation can only be calculated over a limited range of Sh, Sh > rnB^, as shown 



in Fig. 3.12. The larger is the value of A, the smaller is this region. Also, in the range where it 
can be calculated, it depends on the non-perturbative parameter nib (or A). A comparison of this 
distribution and the one in the FM model may be made for the same values of nih and mf^^. This is 
shown for mi, = 4.85 GeV in Fig. 3.18 for HQET (long-short dashed curve) to be compared with 



the dotted curve in the FM model, which corresponds to mf^^ = 4.85 GeV. We see that the two 
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S„ [GeV^] 



Figure 3.18: Hadronic invariant mass spectrum in the Fermi motion model and parton model, based 
on the perturbative contribution only. The solid, dotted, dashed curve corresponds to the parameters 
(Ai, A) = (-0.3,0.5), (-0.1,0.4), (-0.15,0.35) in (GeV"^, GeV), respectively. The parton model (long- 
short dashed) curve is drawn for nib = 4.85 GeV. 

distributions differ though they are qualitatively similar for larger values of Sh > 3GeV^. 



3.10 LD Contributions in B ^ XJ+£- (II) 

This section is devoted to various aspects of the cc resonance effects. Following the procedure adopted in 



[18], we include the long-distance (LD) resonance effects in the decay B Xsi^i and simply add the 



cc resonant contribution with the perturbative cc contribution expressed through the function g{mc, s) in 



section 3.1.2 (see, eq. (3.18)). Thus, in our method, 



Cf{s) = Cgvis) + Y{S) + Yresis) ■ (3.117) 

The function Yresis) accounts for the cc resonance contribution via B —>■ Xs{J/"^, . . .) ^ Xgl'^t^ 



and can be seen in eq. (3.73). Note that in this approach, the effective coefficient Cg^(s) has a s- 
dependence, which is not entirely due to the propagators in the function Yres{s) as also the perturbative 
cc contribution g{rhc, s) is a function of s. In the resonant region, the perturbative part is not noticeable 
due to the fact that the resonant part in C^^{s) completely dominates. However, when the cc pair is suffi- 
ciently off-shell, the s-dependence of the function Cg^^(s) is not (and should not be) entirely determined 
by the cc resonant contribution. This is the motivation of the representation in eq. ( p.ll7| ). 

We start with an analysis of the constraints from existing data on the FM model parameters. Especially 
the question if the FM reproduces the measured J/'0 momentum distribution in B XsJ/ip will be 
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0.5 1 1.5 2 



Figure 3.19: Momentum distribution of J/^j in the decay B XsJ/ip in the FM model. The solid, 
dotted, dashed curve corresponds to the parameters (Ai,A) = (—0.3, 0.5), (—0.3, 0.53), (—0.38, 0.59) 
in (GeV^, GeV), respectively. The data points are from the CLEO measurements 

investigated. Then we turn to the effect of the Lorentz boost on the hadron spectra including LD effects 
according to eq. ( [3.1 17| ) and present Eh, Sh distributions for B Xgi'^l" decay in the FM. Further, we 
study the uncertainties in the B Xsi^£^ spectra resulting from the ambiguities in the parametrization 
of the LD effects. Differences in the distributions from different approaches to treat the cc resonances 
are shown as well for hadron spectra as for the spectra discussed before, namely, the dilepton invariant 
mass distribution and the FB asymmetry. Finally, the hadronic moments are calculated in the FM and 
compared with the ones in the HQE approach for identical values of equivalent parameters. 

3.10.1 Constraints on the FM model parameters from existing data 

The FM model parameters pp and rUg (equivalently Ai and A) can, in principle, be determined from 
an analysis of the energy spectra in the decays B —^ Xuii^e and S — > + 7, as all of them involve 
the decay of a 6-quark into (an almost) massless (u or s) quark. Assuming that the parameters of the 
FM models are universal, these parameters can also be constrained from the lepton energy spectrum in 
the decay B Xcii'e and from the shape of the J/ip- and ip'- momentum distributions in the decays 
B — > Xs{J/ip, ip')- We review the presently available analyses of the photon- and lepton-energy spectra 
in B decays in the FM model (and HQET, as the two are very similar) and also present an analysis of the 
J/V'-momentum spectrum in i? ^ XsJ/4>. 

• Analysis of the photon energy spectrum mB ^ Xs + ^ 
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Figure 3.20: Hadron energy spectrum in B Xgt'^i^ including the resonance and perturbative contri- 
butions in the Fermi motion model (dotted curve) for (Ai, A) = (—0.1 GeV^,OA GeV), and in the parton 
model (long-short dashed curve) for m?, = 4.85 GeV. 



The photon energy- and invariant hadronic mass distributions in i? — > Xs'y were calculated in the FM 
model using the leading order (in q^) corrections in [24,|]. These spectra were used in the analysis of the 
CLEO data on ^ + 7 [ 11], in which the values pp = 270 it 40 MeV suggested by the analysis of 
the CLEO data on i? ^ Xii'i were used, together with the effective 6-quark mass m^^^ = 4.87 ± 0.10 
GeV, which gave reasonable fits of the data. We translate these parameters in terms of Ai and A using the 



relations given in eqs. (3.1 13) and (3.114), yielding 



Ai 



-0.11 



-0.035 
hO.030 



GeV^ A = 0.40 ±0.1 GeV. 



(3.118) 



The same data was fitted in [Q] in the FM model, yielding {pF,nT'q) = (0.45 GeV, GeV) as the best-fit 
solution, with {pp, mq) = (0.310 GeV, 0.3 GeV) differing from the best-fit solution by one unit in x^- 
The quality of the CLEO data [ pT| ] is not good enough to draw very quantitative conclusions. The best-fit 
values translate into 

Ai = -0.3GeV2, A = 0.5 GeV. (3.119) 
• Analysis of the lepton energy spectrum m B ^ Xlv^ 
A fit of the lepton energy spectrum in the semileptonic decay B 



Xlvi in the context of HQET has been 



performed in [|31[]. Using the CLEO data [ |1 13[ ], the authors of [ |31[ ] find: 



Ai 



-0.19 ± 0.10 GeV% A = 0.39 ±0.11 GeV 



(3.120) 



Since the FM model and HQET yield very similar lepton energy spectra apart from the end-point, one can 
take the analysis of also holding approximately for the FM model. 
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Figure 3.21: Hadronic invariant mass spectrum in B Xgi'^l including the perturbative and res- 
onance contributions in the Fermi motion model. The solid, dotted, dashed curve corresponds to the 
parameters (Ai, A) = (-0.3,0.5), (-0.1,0.4), (-0.15,0.35) in (GeV"^, GeV), respectively. 



• Analysis of the J/V'-momentum spectrum in B ^ Xg J /ip 
An analysis of the J/V'-momentum spectrum in B ^ Xs( J/V', ip') measured by the CLEO collaboration 



Ql 120 in the FM model has been reported in [ ]114| ]. The authors of [ |114| ] addressed both the shape and 
normaUzation of the J/V'-data, using the NRQCD formalism for the inclusive color singlet and color 
octet charmonium production in i? —^ XsJ/ip and the FM model. The preferred FM parameters from 
this analysis are: {pp,mg) = (0.57 GeV, 0.15 GeV), where nig only plays a role in determining the 
position of the peak but otherwise does not influence the small momentum tail of the J/ip momentum 



distribution. This yields values of the parameter pi? which are consistent with the ones obtained in |115| 
0.541q |5, GeV based on an analysis of the CLEO data on S — > Xiu^ [113]. The central values of 



PF 



Pf in Oll5| ] as well as in Jl 14| ] correspond to ~ 4.6 GeV, which is on the lower side of the present 



theoretical estimate of m;, pole mass, namely mb = 4.8 ± 0.2 GeV 

We have redone an analysis of the J / ip-momentum distribution which is shown in Fig. 3.19[ As 



shown in this figure, and also discussed in [114], the low-momentum J/ip, in particular in the region 
< 0.6 GeV, are problematic for inclusive decay models, including also the FM model. The \kji^\- 
spectrum appears to have a secondary bump; an inclusive spectrum behaving as a Gaussian tail or having 
a power-like behavior oc \kj/^\''^ in this region is hard put to explain this data. There are also suggestions 
in literature [ ]116| ] that the spectrum in this region is dominated by the three-body decay B — > J /ipKp 
and hence the bump reflects the underlying dynamics of this exclusive decay. In view of this, we have 
taken out the first six points in the low-|/cj/^| spectrum and fitted the FM model parameters in the rest 
of the I A; j/^ I -spectrum. The three curves shown correspond to the FM model parameters {pp^niq) = 
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Figure 3.22: The real part (a) and the absolute value (b) ofCg {s) are shown as a function of s, where 
Cg^^{s) = Cgr](s) + Y(s) + Yresis)- The solid line corresponds to Y{s) calculated using the perturba- 



tive cc contribution g{ific, s) given in eq. (3.18), and the dotted curve corresponds to using g^ific, s) in 



eq. (3.121 ). The dashed one corresponds to the approach by [ 28 /. 



(0.45 GeV,0 GeV) (solid curve), {pF,mq) = (0.45 GeV,0.15 GeV) (dotted curve) and {pF,mq) = 
(0.50 GeV, 0.15 GeV) (dashed curve). They all have reasonable x^> with x^/dof = 1.6,1.6 and 1.1, 
respectively. Excluding also the seventh lowest point, the improves marginally, with the resulting 
being x^/dof = 1.4, 1.4 and 0.94. Including the sixth point, the fits become slightly worse. However, they 
are all acceptable fits. It is interesting that the best-fit solution of the photon energy spectrum in B ^ Xs + 
7, {pf, THq) = (0.45 GeV, GeV), is also an acceptable fit of the |fcj/^|-data. The corresponding Ai, A 
and rrife values from these two analyses are compatible within ±la with the HQET-based constraints from 
the semileptonic B decays ^]\\, quoted above. Thus, the values in eq. ( 3.1 19 ) appear to be a reasonable 



guess of the FM model parameters. But, more importantly for the present study, the phenomenological 
profile of the LD contribution B — > Xs{J /il),^'' , ■■■) — > Xst'^i^ presented here is certainly consistent 
with present data and theoretical constraints. 



3.10.2 Effects of the Lorentz boost on the hadron spectra in S — X^f+f 

We now discuss the B-meson wave function effects in the FM model on the hadron spectra m B ^ 
Xs£^i~ . Since the resonances in i? ^ Xsi^i^ are in the dilepton invariant mass variable s and not in 
Sh, and noting that neither £"0 (partonic energy) nor Eh are Lorentz-invariant quantities, it is expected 
on general grounds that the effect of the Lorentz boost in the FM model on £'//-and S"// -distributions will 
be more marked than what was found on the invariant dilepton mass spectrum in [|l^]. We recall that for 
the dilepton invariant mass, the Lorentz boost involved in the FM model leaves the spectrum invariant 
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Figure 3.23: Hadron energy spectrum in B Xgi^i^ including the resonance and perturbative 
contributions in the Fermi motion model. In (a), the FM model parameters are fixed at (Ai,A) = 
(—0.1 GeV^,OA GeV). The almost overlapping curves differ in the perturbative cc contribution with the 



solid curve obtained using eq. ( p.iq )/or g{rhc, s), the dotted curve using g{rhc, s) given in eq. ^3.121]) . 
The dashed curve corresponds to the approach by /|2^. In (b), the solid, dotted, dashed curve corresponds 
to the parameters (Ai, A) = (-0.3, 0.5), (-0.1, 0.4), (-0.15, 0.35) in (GeV^, GeV), respectively. 



and there is only a subleading effect due to the momentum dependent 6-quark mass. Not so in the hadron 

spectra. In the hadron energy spectrum, the cc-resonances, which are narrowly peaked in the parton model, 

are broadened by the Lorentz boost of the FM model. To show this, the hadron energy spectrum in the FM 

,eff 



model is compared with the spectrum in the parton model in Fig. 3.2C for identical values of and , 
taken as 4.85 GeV. In terms of the hadronic invariant mass, the resonance structure is greatly smeared. The 
reason is that each g^-bin contributes to a range of Eh and Sh- The different-g^ regions overlap in Sh 



resulting in a smearing of the resonances over a wide range. This can be seen in Fig. 3.21 for the hadronic 
invariant mass. Various curves illustrate the sensitivity of this spectrum on the FM model parameters. 



3.10.3 Ambiguities in adding LD and SD contributions in 5 ^ X^^"*"^ 

Since we are simply adding the short-distance (SD) and resonant charmonium amplitudes, it can not be 
ruled out that possibly some double counting has crept in in the coefficient C^^{s) , once as a continuum 
cc contribution and then again as J/?/', V'') ••• resonances. In the absence of a clear separation of the LD 
and SD physics in the spectra, we can not plead the case one way or the other. In the meanwhile, the 
question is whether the addition of the cc-continuum and resonating pieces as being done here and in 



[18] compromises the resulting theoretical precision significantly. This can only be studied by comparing 



the theoretical scenario in question with other trial constructions which have no cc double counting. For 
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Figure 3.24: Dilepton invariant spectrum (a) and the (normalized) Forward-Backward asymmetry (b) in 
B Xsi^l^ including the resonance and perturbative contributions in the Fermi motion model. The 
FM model parameters are fixed at (Ai, A) = (—0.1 GeV^, 0.4 GeV). The curves differ in the perturbative 



CO 



contribution with the solid curve obtained using eq. for g{mc,s), and the dotted curve using 



g{rhc, s) given in eq. (3.121). The dashed curve corresponds to the approach given in ref. [28]. 



example, one could retain in the perturbative function g{rhc, s) just the constant part in s by replacing 

g{mc, s) by g{mc, s), where 



8 ITlb 

9 // 



— In rhr H . 

9 27 



(3.121) 



This function (with ^ = mi,) has been proposed in [^,28] as an alternative representation of the cc 
perturbative contribution and represents the (minimal) short-distance contribution. We denote this ansatz 
for Cf^ defined as C^^is) = C9r]{s) + Y g{mc, s) g{mc,s)") + Yresis) by LSW. Another 
approach is based on a dispersion relation, as proposed by [p8|], here and in the following denoted by KS. 
The KS parametrization of the cc resonant part differs from ours (eq. ( 3.72| )), and the non-resonant part has 



been extracted from data, see [28] for details. The advantage of the KS procedure is that there is certainly 
no double counting. To study the difference numerically, we plot both the real part i?eCg^^(s) and the 
absolute value 1 6*9^^(5)! as functions of s in Fig. 3.22| by using the complete perturbative expression for 
g{rhc, s) in eq. ( 3.16D (our approach) and g{rhc, s) given in eq. (3.121) (LSW) and the KS parametrization. 
Both figures (a) and (b) show that the KS curve is always below our, which is again always below the 



LSW one. One sees from Figs. 3.22 (a) and 3.22 (b) that the difference in these functions in the variable 
s is visible. However, the three parametrizations of the perturbative cc part give almost identical hadron 
spectra, with the resulting uncertainty in the hadron energy and the hadronic invariant mass spectra being at 
most 12.1(4.5)% and 4.1(2.5)%, respectively. The difference between our approach and the KS one (first 
numbers) is larger than the one between ours and the LSW one, given in parentheses. These differences are 
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already difficult to see in tlie hadron energy spectrum in Fig. 3.23 (a); the effect on tiie hadronic invariant 



mass is even less noticeable and hence is not shown. Since, other uncertainties on the hadronic distribution 



are much larger, see, for example. Fig. 3.23 (b) showing the sensitivity of the hadron energy spectra on the 



i?-meson wave function parameters, the much talked about cc-continuum related ambiguity in literature 



is numerically small. Fig. 3.23 shows that it is not the dominant uncertainty in predicting the theoretical 



profiles of hadron spectra in B ^ Xgl^t 



LD uncertainties in the dilepton invariant mass and the FB asymmetry 

We further analyse the uncertainties resulting from different parametrizations of the short and long- 
distance amplitudes in g^-spectra which we have investigated in the previous section |3.4|. We show in 



Fig. 3.24 the dilepton invariant mass (a) and the Forward-Backward (FB) asymmetry (b) with all three 



"SD+LD" approaches discussed before. Unfortunately, the difference between our and the KS one [28] 
in the dilepton spectrum is maximal in the low-q^ region below the J/ip-peak, < 9GeV^. It amounts 
up to 15%. In the other distribution, the FB asymmetry, the difference is found to be moderate over the 
full -range . It does not exceed 5% in the range 5GeV2 < g2 < gCeV^. However, we remark here that 
the position of the first zero of the FB asymmetry is affected by the parametrization of the cc states. 

3.10.4 Numerical estimates of the hadronic moments in FM model and HQET 

To underline the similarity of the HQET and FM descriptions in S — > Xsi~^i~ , and also to make com- 
parison with data, we have calculated the hadronic moments in the FM model using the spectra which we 
have presented in the previous sections. The moments based on the SD-contribution are defined as: 

(X^) = ( J Xl^dXH)/B fov X = S,E. (3.122) 

The moments {X^)cc are defined by taking into account in addition to the SD-contribution also the con- 
tributions from the cc resonances. The values of the moments in both the HQET approach and the FM for 



n = 1, 2 are shown in Table 3.4, with the numbers in the parentheses corresponding to the former. They 



are again based on using the central values of the parameters given in Table |A.1|, and are calculated for 



the same values of the HQET parameters A and Ai, using the transcriptions given in eqs. (3.113). Both 
the HQET and the FM model lead to strikingly similar results for the hadronic moments shown in this 
table. However, the moments {Xjj)cc with X = S, E are significantly lower than their SD-counterparts 
(Xj^) calculated for the same values of the FM model parameters. This shows, at least in this model 
study, that the cc resonances are important also in moments. The hadronic invariant mass spectra in 
B — > Xs£'^i~ for both the SD and inclusive contributions are expected to be dominated by multi-body 
states, with (Sh) ^ (1-5 - 2.1) GeV^ and {Sh)cc ^ (1-2 - 1.5) GeV^ . Note that the difference in the nu- 



merical values of the hadronic mass moments {Sh)cc and {Sjj)cc shown in Table 3.4 caused by different 



LD parametrizations is less than 0.22%, 0.42%, respectively, as can be see in Table 3.6 
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(Sh) 


( ) cc 




i. S TJ ^ r'r- 


(Ai, A) in (GeV^, GeV) 


(GeV 


^) 


(GeV 


^) 
I 


(-0.3,0.5) 


2.03 (2.09) 


1.51 


6.43 (6.93) 


3.10 


(-0.1,0.4) 


1.75 (1.80) 


1.36 


4.04 (4.38) 


2.17 


(-0.14,0.35) 


1.54(1.49) 


1.19 


3.65 (3.64) 


1.92 




{Eh) 


{Eh)cc 


\ H 1 


{ EI-T ") nn 

\ H 1 cc 


(Ai, A) in (GeV^, GeV) 


(GeV) 


(GeV 




(-0.3, 0.5) 


2.23 (2.28) 


1.87 


5.27 (5.46) 


3.52 


(-0.1,0.4) 


2.21 (2.22) 


1.85 


5.19(5.23) 


3.43 


(-0.14,0.35) 


2.15 (2.18) 


1.84 


4.94 (5.04) 


3.39 



Table 3.4: Hadronic spectral moments for B 
indicated values of the parameters (Ai, A). 



Xgfi^H in the Fermi motion model (HQET) for the 



3.11 Branching Ratios and Hadron Spectra in B 
Invariant Masses 



XJ+i- with Cuts on 



In experimental searches for the decay B —^ Xgi'^i^ , the short-distance contribution is expected to be 
visible away from the resonances. So, cuts on the invariant dilepton mass are imposed to get rid of the 
dilepton mass range where the charmonium resonances J/ip and ip' are dominant. For example, the cuts 



imposed in the recent CLEO analysis |3C] given below are typical: 



cut^ : g2 < (^^^^ _ 0.1 GeV)2 = 8.98 GeV^, 
cut B : < {mj/^ - 0.3 GeV)^ = 7.82 GeV^ , 
cutC : g2 > + 0.1GeV)2 = 14.33 GeV^ . 



(3.123) 



The cuts A and B have been chosen to take into account the QED radiative corrections as these effects 
are different in the e~^e~ and modes. In the following, we compare the hadron spectra with and 

without the resonances after imposing these experimental cuts. For the low-g^ cut for muons (cut A), the 



hadron energy spectra and the hadronic invariant mass spectra are shown in Fig. 3.25 (a), (b) and Fig. 3.26 
(a), (b), respectively. The results for the low-q^ cut for electrons (cut B), are shown in Fig. 3.25 (c), (d) 
and Fig. 3.26| (c), (d), respectively. Finally, the hadronic spectra for the high-g^ cut (cut C) for e+e" and 
can be seen in Fig. 3.25 (e), (f) for the hadronic energy and in Fig. 3.26| (e), (f) for the hadronic 
invariant mass. We see that the above cuts in greatly reduce the resonance contributions. Hence, the 
resulting distributions essentially test the non-resonant cc and short-distance contributions. These figures 
will be used later to quantify the model dependence of the integrated branching ratios in B — > Xs£^i~ ■ 
As mentioned in [^], the dominant BB background to the decay B —>■ Xgl'^l^ comes from two 
semileptonic decays of or I? mesons, which produce the lepton pair with two undetected neutrinos. To 
suppress this BB background, it is required that the invariant mass of the final hadronic state is less than 
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t = 1.8 GeV, which approximately equals ttid. We define the survival probability of the B Xgi^i 
signal after the hadronic invariant mass cut: 

(IB 



Sit) 



.| dSH 



dSH)/B , 



(3.124) 



and present S{t = 1.8 GeV)) as the fraction of the branching ratio for B —>■ Xsi~^£~ surviving these cuts 
in Table |3.5| . To estimate the model dependence of this probability, we vary the FM model parameters. 
Concentrating on the SD piece, we note that the effect of this cut alone is that between 83% to 92% 
of the signal for B —>■ Xg^'^ ji" and between 79% to 90% of the signal m B ^ Xge^e" survives, 
depending on the FM model parameters. The corresponding numbers for the inclusive spectrum including 
the SD and LD contribution, here and in the following abbreviated as tot =SD+LD, is 96% to 99.7% for 
both the dimuon and dielectron case. This shows that while this cut removes a good fraction of the BB 
background, it allows a very large fraction of the B — > Xst^£~ signal to survive. However, this cut does 
not discriminate between the SD and (SD+LD) contributions, for which the cuts ^ - C are effective. The 
numbers for the survival probability S{t = 1.8 GeV) reflect that the hadronic invariant mass distribution 
of the LD-contribution is more steep than the one from the SD contribution. 



(Ai,A) 


B ■ 10-6 


B ■ 10-6 


No s-cut 


No s-cut 


cut A 


cutB 


cut C 


cutC 


GeV2, GeV 












e^e- 




e^e- 


(-0.3,0.5) 


5.8 


8.6 
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79 % 


57% 


57% 


6.4% 


4.5% 


(-0.1,0.4) 


5.7 


8.4 


93% 
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63% 


68% 


8.3% 


5.8% 


(-0.14,0.35) 


5.6 


8.3 


92% 


90% 


65% 


67% 


7.9% 


5.5% 


(-0.3,0.5)tot 


562.5 


563.9 


96% 


96 % 


0.8% 


1.0% 


0.06% 


0.06% 


(-0.1,0.4)tot 


564.0 


565.6 


99.7% 


99.7% 


0.8% 


1.1% 


0.08% 


0.08% 


(-0.14,0.35)tot 


566.5 


568.2 


99% 


99 % 


0.9% 


1.2% 


0.08% 


0.08% 



Table 3.5: Branching ratios for B Xgl^l" , i = fi, efor different FM model parameters are given in 
the second and third columns. The values given in percentage in the fourth to ninth columns represent the 
survival probability S{t = 1.8 GeV) defined in eq. (3.124) for different FM model parameters and cuts 



on the dilepton invariant mass as defined in eq. (3.123). The subscript tot = SD + LD denotes that both 
the short and the long-distance contributions are included in the branching ratios and S{t). 

With the additional cut A (B) imposed on the dimuon (dielectron) invariant mass, between 57% to 
65% (57% to 68%) of the B —>■ Xsl'^i~ signal survives the additional cut on the hadronic invariant mass 
for the SD contribution. However, as expected, the cuts A and B result in drastic reduction of the inclusive 
branching ratio for the decay B — > , as they effectively remove the dominant cc-resonant part. 

In this case only 0.8% to 0.9% (1.0% to 1.2% of the inclusive signal survives for the cut A (B). The 
theoretical branching ratios for both the dielectron and dimuon cases, calculated using the central values 



in Table A.l are also given in Table 3.5. As estimated in eq. (3.71), the uncertainty on the branching 



ratios resulting from the errors on the parameters in Table A.l is about ±23% (for the dielectron mode) 
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and ±16% (for the dimuon case). The wave function-related uncertainty in the branching ratios is smaller, 



as can be seen in Table 3.5. With the help of the theoretical branching ratio and the survival probability 
S{t = 1.8) GeV, calculated for three sets of the FM parameters, the cross section can be calculated for all 
six cases: 

(i) no cut on the dimuon invariant mass [(SD) and (SD + LD)], (ii) no cut on the dielectron invariant mass 
[(SD) and (SD + LD)], (iii) cut A on the dimuon invariant mass [(SD) and (SD + LD)], (iv) cut B on the 
dielectron invariant mass [(SD) and (SD + LD)], (v) cut C on the dimuon invariant mass [(SD) and (SD + 
LD)], (vi) cut C on the dielectron invariant mass [(SD) and (SD + LD)]. This gives a fair estimate of the 
theoretical uncertainties on the partially integrated branching ratios from the S-meson wave function and 
cc resonances. This table shows that with 10^ BB events, C(70) dimuon and (0(100) dielectron) signal 
events from B Xs£'^e~ should survive the CLEO cuts A (B) with m{Xs) < 1.8 GeV. With cut C, 
one expects an order of magnitude less events, making this region interesting for the LHC experiments. 
Given enough data, one can compare the experimental distributions in B ^ Xs£^i~ directly with the 
ones presented here. 



3.11.1 Hadronic spectral moments with cuts in the FM 



We have calculated the first two moments of the hadronic invariant mass in the FM model by imposing a 
cut Sh < with t = 1.8 GeV and an optional cut on q^. 



{S^h) = ( 



,1 " dSudq^ 



dSndq^ 



dSudq^ 



for n = 1, 2 



(3.125) 



Here the subscript cutX indicates whether we evaluated {Sh) and {S'jj) with the cuts on the invariant 
dilepton mass as defined in eq. ( 3.123 ), or without any cut on the dilepton mass. The results are collected 



in Table |3.6| . The moments given in Table |3^ can be compared directly with the data to extract the FM 
model parameters. The entries in this table give a fairly good idea of what the effects of the experimental 
cuts on the corresponding moments in HQET will be, as the FM and HQET yield very similar moments for 
equivalent values of the parameters. The functional dependence of the hadronic moments on the HQET 
parameters taking into account the experimental cuts still remains to be worked out. 

Further, we have calculated {Sh) and {S'jj) with a cut Sh < 3.24 GeV^ and optional ones on q"^ (cut 
A-C according to eq. ( pl3| )) with the approaches KS [|28|] and LSW ^ for (Ai, A) = (-0.1, 0.4) in 
GeV^, GeV. They differ from ours (eq. (3.117)) in the parametrization of the resonant and non-resonant 



cc contributions, as discussed in section |3.10.3| . We compare the values of the moments for the same set 
of FM parameters. Denoting our approach by y, we define by A the maximal deviation in % between 
ours and KS and LSW, generically written as: A = max{\y — LSW|/|y|, \y — KS|/|y|) and present it 



in the last row of Table [3.6| . We see that the uncertainties in the hadronic mass moments from different 
"SD+LD" parametrizations are small, namely below 1.6% in the worst case. 



3.12 Summary and Concluding Remarks on the Decay B Xgi'^i 

In this chapter we have investigated distributions, decay rates and moments in rare inclusive B Xsi^i^ 
decay in the standard model. In the first part, we have concentrated mainly on the study of distributions in 
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lA 




(-0.14,0.35) 


J. J. 


2 34 


1 38 


2 55 


1.47 


2 83 


1 52 


2 97 


0.66 


0.54 


(—0 3 5)/„/ 

V vj.o, yj-^ Jtot 


1.41 


2.61 


1.41 


2.62 


1.61 


3.32 


1.66 


3.41 


0.74 


0.68 


(-0.1,0.4)/n/ 


1.35 


2.14 


1.36 


2.15 


1.77 


3.60 


1.87 


3.94 


0.74 


0.62 


(-0.14,0.35)tot 


1.17 


1.84 


1.18 


1.85 


1.45 


2.16 


1.51 


2.95 


0.66 


0.54 


A(%) 


0.15 


0.19 


0.22 


0.42 


0.90 


1.56 


0.32 


0.58 


0.01 


0.32 



Table 3.6: {Sh) and {Sfj) for B Xsi'^i , £ = fi, efor different FM model parameters and a hadronic 
invariant mass cut Sh < S.2AGeV^ are given with and without additional cuts on the dilepton invariant 



mass as defined in eq. (' p.i2j| ). The Sn-tnoments with cuts are defined in eq. (' p.i25| ). The subscript 
tot = SD + LD denotes that both the short and the long-distance contributions are included in these 
moments. The value of A estimates the uncertainty from different approaches to take into account the 
effect of the cc continuum and resonances, see text. 

the dilepton invariant mass q^, the differential branching ratio in this variable and the FB asymmetry. Our 
findings can be summarized as follows [|l^] : 

• We have calculated the leading 1 / mf, power corrections with HQE techniques in the dilepton in- 
variant mass distribution in B ^ Xsi~^£~ decay and have explicitly kept the s-quark dependence. 
Our calculation is at variance with an earlier one [17] in the limit = and has been confirmed 
recently by [^] for the massless s-quark case. 

• We find that the 1 / corrections are stable over a good part of the dilepton mass spectrum. How- 
ever, near the high end-point the distribution becomes unphysical due to the HQE parameter A2, 
signaling a breakdown of the heavy quark expansion. 

• The B Xgi^i^ decay rate in the HQE method decreases by about 4% and the branching ratio 
by about 1.5% from their corresponding parton model values. 

• Alternatively, we have implemented a Gaussian Fermi motion model in the decay B — > Xsi~^£~ to 
model i?-meson bound state effects. We have analysed the dilepton invariant mass distribution and 
the FB asymmetry within this framework, showing the dependence on the parameters of this model. 
Non-perturbative effects are found to be perceptible in both distributions in the high q^ region. 

• The theoretical uncertainties in the short-distance branching ratios in i? ^ Xsi^i^ decay are found 
to be ±23% (±16%) for the electron (muon) case in the FM. 
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• We have modeled the long-distance contributions from intermediate charmonium resonances with 
a VMD ansatz and presented the dilepton invariant mass distribution and the FB asymmetry in the 
FM, including next-to-leading order perturbative QCD corrections in figures. 

We have completed the description of final states in B ^ Xgi'^i^ decay in the second part of this chapter, 
which is devoted to the study of hadron spectra and hadronic spectral moments. We summarize [ p5|j26| , p4| ] : 

• We have calculated the 0{as) perturbative QCD and leading 0{l/mi,) corrections to the hadron 
spectra in the decay B Xgi^i" , including the Sudakov-improvements in the perturbative part. 

• We find that the hadronic invariant mass spectrum is calculable in HQET over a limited range 
Sh > msA. and it depends sensitively on the parameter A (equivalently mf,). These features are 
qualitatively very similar to the ones found for the hadronic invariant mass spectrum in the decay 
B ^ XJue rtToTI]. 



The l/?7ifc -corrections to the parton model hadron energy spectrum in i? — > Xsi~^£~ are small over 
most part of this spectrum. However, heavy quark expansion breaks down near the lower end-point 
of this spectrum and close to the cc threshold. The behavior in the former case has a similar origin 
as the breakdown of HQET near the high end-point in the dilepton invariant mass spectrum, which 



we have presented here and in ref. [18] 



We have calculated the hadronic spectral moments {Sjj) and (Ej^) for n = 1,2 using HQET. The 
dependence of these moments on the HQET parameters is worked out numerically. In particular, the 
moments {Sjj) are sensitive to the parameters A and Ai and they provide complementary constraints 
on them than the ones following from the analysis of the decay B — > X£v£. The simultaneous fit 
of the data in i? ^ Xsi'^i' and B — > Xluf^ could then be used to determine these parameters very 



precisely. This is illustrated in Fig. 3.15 



The corrections to the hadron energy moments A{EH)i/m^ A{Ejj)i/^2 from the leading 
0(AQ(^£,/m^) power corrections have been worked out, using the results of [27]. We find that 
these corrections are very small. The corresponding corrections in A(5^)i/„2 vanish in the theo- 
retical accuracy we are working. 

We think that the quantitative knowledge of A and Ai from the moments can be used to remove 
much of the theoretical uncertainties in the partially integrated decay rates in B ^ X^ii^e and 
B —>■ Xsi^i^ ■ Relating the two decay rates would enable a precise determination of the CKM 
matrix element Vub- 

As a phenomenological alternative to HQET, we have worked out the hadron spectra and moments 
in S — > Xs£~^i~ in the Fermi motion model. We find that the hadron energy spectrum is stable 
against the variation of the FM model parameters. However, the hadronic invariant mass is sensitive 
to the input parameters. Present theoretical dispersion on this spectrum can be reduced with the 
experimental measurements of the corresponding spectrum in the decay B — > Xuii^e, which will 
determine these parameters. Conversely, with good measurements of the decay B Xgi^i^ , one 
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could fix the input parameters in tlie decay B Xuii'e and determine the CKM matrix element 
Vub rather precisely. 

For equivalent values of the FM and HQET parameters, the hadronic spectral moments in B ^ 
Xsi^l^ decay are found to be remarkably close to each other. 

We have worked out the effect of the cc resonances in the effective coefficient ^^(s) on the hadron 
spectra in B Xsi^i^ , by parametrizing the present data on the resonant part from the decays 
B Xs{J /ijj^i^' , •••)■ The resonances are incorporated at the parton level and the broadening 
of these resonances from the wave function effects in the FM model are then worked out. These 
spectra will provide an important test of the FM model in B ^ Xsi~^£~ ■ 

We find that the cc long-distance effects in i? ^ Xgi^i^ decay are also important in the hadronic 
moments. 

We have quantitatively studied possible double counting effects which may have entered in sim- 
ply adding the resonant contribution via Breit-Wigner functions and the complete perturbative cc 
contribution in the coefficient Cf^^(s) in B ^ Xsi'^i~ . The numerical difference between this 



approach, followed here [18,14|, and alternative ones |28] and [p^], are found to be small in the 
dilepton invariant mass spectrum and negligible in the hadron spectra and moments. Theoretical 
spectra are found to be more sensitive to the wave function effects, which dominate the uncertainty 
in the shape. 

• We have worked out the hadron spectra by imposing the experimental cuts designed to suppress the 
resonant cc contributions, as well as the dominant BB background leading to the final state BB 
Xsi^£~ {+ missing energy). The parametric dependences of the resulting spectra are studied in the 
FM model. In particular, the survival probability of the B —^ Xst^l~ signal by imposing a cut on 
the hadronic invariant mass Sh < 3.24 GeV^, as used in the CLEO analysis, is estimated and its 
model dependence studied. This quantifies the statement that with the indicated cuts, these spectra 
essentially test the physics of the short-distance (and non-resonant cc) contribution. 

The CLEO collaboration has already been searched for inclusive b — > s£~^i~ decay with i = e, ^. Their 



results are [30]: 



B{b^ se+e-f^^^ < 5.7 -10-^ 

B{b sfi+fi-f'^^^ < 5.8 • IQ-^ . (3.126) 
Comparing this with our estimates for the SD branching ratios in the decay B —>■ Xs£~^i~ given in 



eq. (3.71), we see that the above CLEO upper bounds are approximately one order of magnitude away 



from the theoretical standard model prediction. 
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(a) (b) 




Figure 3.25: Hadron energy spectrum in B X.sl'^i in the Fermi motion model with the cuts on the 



dilepton mass defined in eq. (3.123); (a),(c),(e) without and (b),(d),(f) with the cc-resonance contribution 
corresponding to cut A,B,C, respectively. The solid, dotted, dashed curves correspond to the parameters 
(Ai, A) = (-0.3,0.5), (-0.1,0.4), (-0.15,0.35) in (GeV^, GeV), respectively 
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Figure 3.26: Hadronic invariant mass spectrum in B Xgi'^i in the Fermi motion model with the 



cuts on the dilepton mass defined in eq. (3.123); (a),(c),(e) without and (b),(d),(f) with the cc-resonance 
contribution corresponding to cut A,B,C, respectively. The solid, dotted, dashed curves correspond to the 
parameters (Ai, A) = (-0.3,0.5), (-0.1,0.4), (-0.15,0.35) in (GeV^, GeV), respectively. 



Chapter 4 



The Decay Bg 77 



Besides the rare decays B —>■ Xgj and B —>■ Xsi~^£~ , i?s ^ 77 is another potential candidate to explore 
perturbative and non-perturbative aspects of QCD and test the standard model (SM). The L3 collaboration 
has already been searched for Bs —>■ 77 decay. Their upper bound is the best present limit of this channel 



[117] 



J3(Bs 77) < 1.48 • 10"^ . (4.1) 
The first theoretical analysis of a rare decay into 2 photons is contained in the pioneering work by Gaillard 



and Lee [118], who considered Ks^l 77- Exclusive Bs — > 77 decay has been investigated in the 



lowest order in refs. [^2[-35]. The branching ratio found is 4.5 • 10^"^ in the SM context for = 0.5 
GeV and other parameters given in Table |A. 1[ The large value of the s-quark mass here results from using 
the constituent quark mass mg ~ mx- As learned from studies of i? — > Xgj decay, the flavour changing 
neutral current (FCNC) b —>■ vertex receives in leading logarithmic (LLog) QCD approximation a 
large enhancement about a factor of 2-3. We suggest similar effects for the b — > 577 transition. To see 
whether Bs —>■ 77 decay is worth more effort to be searched for at future experiments, a more advanced 
analysis of its branching ratio is required. A branching ratio of order 10^^ is a benchmark of a decay 
to be measured at present B experiments like CLEO or Hera-B with reasonable statistics. Upgrades and 
planned i?-factories will be sensitive to branching ratios of order 10^^. However, the BsBs pair is too 
heavy to be produced at the T(4s) resonance. 

The Bs 77 final state consists of a CP-odd and a CP-even T+ amplitude. It offers, besides 
the branching ratio, another observable, the so called CP ratio |T+p/|r~p [ ^4p5| ]. With it CP violating 
effects can be studied. 

The work reported in this thesis, which has already been published [^], [ p8[ ] differs from the previous 



ones [35,33,33] with respect to three points: a) We calculate and use the QCD-improved LLog amplitudes. 



b) In contrast to previous works using the constituent quark model we model the bound state effects of 



the Bs meson through an heavy quark expansion technique (HQET) inspired approach following [16]. 
This introduces an additional dispersion on the Bs 77 branching ratio and CP ratio, which serves as 
an estimate of the hadronic uncertainties, c) We include long distance effects due to decay chains via 
intermediate vector mesons in our analysis. To be definite, we estimate the additional contribution in the 



decay Bs 77 through Bs 07 followed by — > 7 using Vector Meson Dominance (VMD) [ 74 ] 
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Further, the Bg — > </>'0 decay is modeled by inclusive b ^ sip decay. Using the VMD model, the amplitude 
for the chain process Bg — > (pi/j —>■ (pj —>■ 'yj is presented. 

Leading logarithmic QCD corrections for the short-distance part of the decay Bs 77 have also 



been calculated by Chang et al. [119|. They fix A = — irib, which corresponds to rus in the naive 
constituent quark model, and the renormalization scale /x = mi,. We emphasize here that the decay rate 
(and the CP ratio) is sensitive to both of these parameters and requires further theoretical investigation. 
Soni et al. [ ]120| ] calculated also LLog QCD calculations in Bg —>■ 77 and in addition in the decay 
B Xg^j. Analyses of Bg 77 decay in non standard models have been done in ref. [RSl] in the 



lowest order and in ref. [121] including LLog QCD corrections in the two Higgs doublet model (2HDM), 



in ref. [122] in the minimal supersymmetric model (MSSM) and in ref. [123] in the 2HDM with flavour 



changing neutral currents allowed at tree level. 

4.1 Leading Logarithmic Improved Short-Distance Contributions in Bg 
77 Decay 

In this section we present the leading logarithmic QCD-improved rates for exclusive Bg 77 decay. We 
use the free quark model and make the connection between quark and mesons states by means of the Bg 
meson decay constant Jb^ ■ Dot products of kinematical variables, which are not fixed by this are estimated 
in an approach inspired by HQET. 

QCD-improved rates in 6-quark decays can in general be obtained through the following procedure: 
Matching of the full theory with an effective theory at a scale /x = mw, using an effective Hamiltonian 
and performing an evolution of the Wilson coefficients from mw down to /i ~ 0{mb), thus resumming all 
large logarithms of the form a" (m;,) ln™(^^), where m < n (n = 0, 1, 2, . . .). In the leading logarithmic 
approximation, which we use here, m = n. In our case, which is Bg —>■ 77 now an enormous short cut is 



possible from observing that the effective Hamiltonian in eq. ( |2. 18| ) for 6 — > 57 is identical for b — > 577 
to this order of — ^ : 

Heffib ^ S7) = neffib ^ S77) + O(^) . (4.2) 

The proof goes as follows: One can ask, if there are more operators needed for b — > 577 than included 
in Tleff{b —>■ 57) and try to find an operator sXb with dim{X) < 3 and containing two photons! Here X 
must be a gauge and Lorentz invariant structure made out of quark and photon fields, masses and covariant 
derivatives = d^ + ieQqA^. (For the moment we shall work in zeroth order of the strong interactions.) 
When constructing the full set of physical operators, the equations of motion (EOM) can be used to reduce 
the operator basis: 

ipq = mgq, D ^F'^'^ = eQ^qy" q , d^F^" = . (4.3) 
Here F^^'^ , F^^^ = ^e^^"^^ Fap denote the photon field strength tensor and its dual, respectively. For chiral 



fermions we have the following EOM, which can be obtained with the identities given in eq. (A.41) 



ipbL = mbbji , SLtp = rUgsL = rUgSR . (4.4) 
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Other useful identities in this context are: 

D'^ = -]^eQqaF , D,, = ^{p-f^ + j^p) , [D^, D,] = ieQgF^, . (4.5) 



As a result, using the EOM gives either a mass, a current or remains in the operator basis O1...8 eq. (2.18) 
or gives contributions to the FCNC self energy. The latter will be absorbed in the on-shell renormalization 
and does not give any contribution to 6 ^ 577 

= ml)b = sY:{p'^ = m^) = . (4.6) 

Here p, p' is the 4-momentum of the incoming 6-quark, outgoing s-quark, respectively. To display the 
foregoing we give some examples of EOM operator identities: 

• s{ip)^b = mlsb 

• nibsD'^b = -mlsb - ^eQgSaFb 

. s^,{D^F^''')b = eQgSj,bqrq 
Since after applying the equations of motion there exists no gauge-invariant FCNC-2-photon operator 



with field dimension < 6, the set of operators given in eq. (2.18) is a complete basis for both b —>■ and 



b —>■ S77 decay JTO|,^. Hence all the results obtained for the former, a collection of which can be seen 



in section 2.2, can be used for the latter, like the LLog evolution of the Wilson coefficients. 



4.1.1 ^ 77 decay in the effective Hamiltonian theory 

Having convinced us that the set-up eq. ( |4.2| ) is correct, we can now turn to an explicit calculation of LLog 
improved FCNC 2 photon amplitudes. The amplitude for the decay Bs —>■ 77 can be decomposed as 

[HiB 

A{Bs ^ 77) = 4{h)e'2(.k2){A+g^, + i^-e^,„^A;f fc^) , (4.7) 

where the ki and e\{ki) denote the four-momenta and the polarization vectors of the outgoing photons, 
respectively Alternatively, we can write the amplitude in terms of photon field strength tensors: 

A{Bs ^ 77) = T+F^,Ff + iT-Fl^F^" . (4.8) 

We have for real photons kf = 0, e^./cj = with i = 1,2, ki.k2 = m'^j2 and further for a Bs decaying 
at rest the additional conditions ei.k2 = e2-ki = 0. Then the following equations result: 



-m 



T+ , A- = 2T- . (4.9) 



Also, in the rest frame of the Bs meson, the CP = — 1 amplitude A^ is proportional to the perpendicular 
spin polarization el x ei, and the CP = 1 amplitude is proportional to the parallel spin polarization 
el. €2. The ratio 



4A- 



|2 



rcp = |r+|7|r-|^ = , .' , (4.10) 



'We adopt the convention Tr{'y>'Y'y°''y'^'y5) = 4ie^'''"'^ , with e°^^^ = +1. 
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can be used to study CP violating effects in Bg 77 decays. We will discuss it together with the 



Bs —>■ 77 branching ratio in section 4.4 



Using the effective Hamiltonian in eq. (2.18), the CP-even (^1^) and CP-odd (A ) parts in the SM can 



be written as (for diagrams see Fig. 4.1 and Fig. 4.2) in a HQET inspired approach: 



4 171^2 \ 

y m^^ ■' + m 



V2tt \3mBs-^s[mB, - As) 

- ^g^/(mg)Cg(/x)+ / ^imbA{mb) + msA{ms))Difi)] , (4.11) 

q 9(mj,-'-' +ms-'-') J 



where we have used the unitarity of the CKM-matrix J2i=u c t ^■ts^ib = and have neglected the contri- 
bution due to V*sVub < V-tlVtb = At. In eq. ( [4.1l[ ) Nc is the colour factor {Nc = 3 for QCD) and Qg = | 
for g = ti, c and Qq = — ^ for g = d,s,b. The QCD-corrected Wilson coefficients in leading logarithmic 



approximation which are discussed in section 2.2, Ci...6(/^) and CY^^{fi), enter the amplitudes in 

the combinations 

Cui^^) = Caifi) = (C3(/u) - C5ifi))N, + C^ifi) - Ceifi) , 

Cifl) = + C-sifi) - C5(/i))7Vc + C2i^l) + C^ip) - C6(/U) , 

Csifi) = Cbip) = (C3(/i) + C^imNc + 1) - N.C^iij) - Ce{^i) , 

Dii^) = C^ifx) + Ce{fi)N, . (4.12) 

The Feynman rules used are given in appendix |A.3| . Note that the chromomagnetic operator Og does not 
contribute here in this order of a^. The functions I{mq), J{mq) and A{mq) come from the irreducible 
diagrams with an internal q-type quark propagating, see Fig. [4.1| , and are defined as 

2 

TTl 

I{mq) = 1 + -^A{mq) , 
m%, — Amj. 

AK) = |ln( \ / ^ )-^^ f^^i^^^l' 





4m2 




4m2 


y4m2- 





2 

A(m„) = - I 2arctan(-^^ ^) - vr 1 for < 1. (4.13) 

Detailed analysis shows that the diagram Fig. ^ requires the calculation of three different types of inser- 
tions: The current-current operators Oi^2 only give a contribution to a charm loop. They have the structure 
7^L 7^L, which leads after integration over the internal quark momentum to the function I{mq) given 
above. In contrast, each of the penguin operators O3...6 has two possible insertions, a "direct" one and one 
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after Fierz ordering of the fields. In the former just internal b- and s-quarks appear and the operators have 
to be "turned" by ±90° to generate a diagram consisting of one continuous b s fermion line. In the 
latter the four-Fermi operators, which contribute to all 5 active flavours q = u,d, s, c, b, are rearranged 
with the help of the Fierz transformation given in appendix |A.3| . This simplifies the calculation as it cir- 
cumvents a trace over 7-matrices. The procedure is legitimate since the resulting amplitude is (IR and 
UV) finite. The operators O3 4 have the same Dirac structure as Oi^2 ~ 7^iL "ff_iL which is reproduced 



after Fierz transformation, see eq. (A.21). Therefore, here no new integrals appear. The operators O5 e are 
of j^L^j^R type. "Direct" insertion leads to the functions A(mg), J{mq) given above. They contribute 
to the Bg —>■ 77 amplitude only via an internal s- and 6-quark. Here care must be taken of the left-right 
structure, which is different for the s- and the 6-quark and results in the sign difference in the correspond- 
ing term proportional to D{fi) in the CP-even amplitude A'^ given in eq. ( |4.1 1[ ). The Fierz transformation 



of O5 6 results in a scalar/pseudoscalar coupling ^ R® L, see eq. { ^.22^ . The analytical expression for 



such an insertion is minus the one for 7^L 7^L, which can be checked after explicit calculation. From 



here the minus signs in the functions eq. (4.12) can be understood. 

The parameter enters eq. ( [4. 1 1[ ) through the bound state kinematics. At the quark and meson level, 
the decay kinematics are given 

Kp) ^ S(P )7(A:i, £1)7(^2, £2) , (4.14) 
B,{P) ^ l{ki,eMk2,e2) , (4.15) 



respectively. A problem lies now in the intermediate propagators of the reducible diagrams, see Fig. 4.2, 
where we need to evaluate p.ki and p' .ki, i = 1,2. The answer cannot be given just by using kinematics, 
energy/momentum conservation in a chosen frame and a model is necessary here. For definiteness, we 
consider the decay Bs = (bs) — > 77. We write the momentum of the 6-quark inside the meson as 
p = rubV + k, where /c is a small residual momentum, v is the 4-velocity, which connects the quark with 
the meson kinematics through P = ruB^v and P is the momentum of the meson. In the Bs rest frame, 
V = (1, 0, 0, 0). Now following [p^, we average the residual momentum of the 6-quark through 

<ka> = -77— (^1 + 3A2)Va , 

znib 

<kakf3> = Y(9al3 - VaV/3) , (4.16) 

where Ai, A2 are matrix elements from the heavy quark expansion. Using P = p — p' , P.ki = 
v.ki = ^ and the HQET relation 

ms, =mfe + A, - -^(Ai +3A2) (4.17) 
zrrib 



one gets: 



Jk- - -^A 
h"?^p' = "^g-3A2, 



{ml^ff=p'^ = {rrfjff -ml +2mBAs ■ (4-18) 
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The non-perturbative parameter A^ can be related to A, which has been extracted (together with Ai) from 
data on semileptonic , B^ decays in ref. [31], and the measured mass difference Am = uib^ — ms = 



90 MeV [p9|], defining A^ = A + Am. The matrix element A2 is well determined from the Bf^ — B(^g-^ 



mass splitting, A2 = 0.12 GeV . With the help of eq. (4.17), the correlated values of A and Ai can 
be transcribed into a correlation between A(g) and m^. We select 3 representative values |^ {mi,, As) = 
(5.03, 370), (4.91, 480), (4.79, 590) in (GeV, MeV) to study the hadronic uncertainties of our approach. 
Note that we assume here that Ai, A2 are flavour independent. Furthermore, we have used the definition 

< 0| 57^756! B.(P) > = , (4.19) 

which leads together with the off-shellness of the quarks inside the meson to the matrix element of the 
pseudoscalar current 



m 



<0\sj5b\Bs{P)> = -ifB.^fT^^^ ■ (4.20) 

+ ms 

The auxiliary function = g^{m^J^ ,Ks) is defined as 

g-=mBAml^^ + mlfff + Ks{ml^-{rnJ^ + mlfff) . (4.21) 

Note that in the limit A^ m^, m^J^ — > m\,^s and using m^^ = m\, + m^ we recover the result obtained 
by the constituent quark model [3^33 35], ignoring QCD corrections. Using the above expressions, the 



partial decay width is then given by : 

- — (4|A+|2 + i 
327rmBs 2 



T{Bs ^ 77) = 00-^ (41-4+1' + -m%JA~\^) . (4.22) 



Now, there are 2 new observations to be made: 
First, the Wilson coefficients in eq. ( [4.1 1[ ) depend on the scale fi. Therefore, since the behaviour of these 



short-distance (SD) coefficients under renormalization is known from the studies of 5 — > Xsj I l43| , [44[ 
one can give an improved width for Bs 77 by including the leading logarithmic QCD corrections by 
renormalizing the coefficients C1...6 and Cj^^ from n = mw down to the relevant scale fi ^ 0{mh). 
The explicit 0{as) improvement in the decay width r{Bs —>■ 77) requires the calculation of a large 
number of virtual corrections, which we have not taken into account. Varying the scale in the range 
^ < /i < 2mii, one introduces an uncertainty, which can be reduced only when the complete next-to- 
leading order (NLO)-analysis is available, similar to the recently completed calculation for the B Xsj 
decay [|3],||]. 

The second point concerns the strong dependence of the decay width T{Bs 77) on A^, F ~ C(]p-) 



in eq. ( 4.22). It originates in the s-quark propagator in the diagram with an intermediate s-quark in Fig. nJi 



In the earUer work the authors of e.g. [ ]32| ] evaluated the decay width with mg ~ mx, assuming that the 
constituent quarks are to be treated as static quarks in the meson. This is a questionable assumption. In the 
HQET inspired approach, this gets replaced by A^, which is well-defined experimentally. This formalism 
implies, that the decay width T{Bci — > 77) will involve the parameter A, which avoids the unwanted 
uncertainty on m^. 

^We choose (Ai,A) = (-0.09,280), (-0.19,390), (-0.29,500) in (GeV2,MeV) from Fig. 1 in [0]. 
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Figure 4. 1: The generic diagram contributing toh ^ 577 in the effective theory due to the (Fieri ordered) 
four-quark operators. The diagram with interchanged photons is not shown. 




Figure 4.2: The reducible diagrams contributing to b ^ S77. Tlie blob denotes the FCNC operator O7. 
The diagrams with interchanged photons are not shown. 

The branching ratio and the CP ratio as a function of the scale fi for different values of [rrih, Ag) are 



discussed in section 14 including the Oy-type long-distance (LD) estimate. 

In the lowest order, which can be recovered at fi = myy, Cj and C2 are the only remaining non- 
zero Wilson coefficients. The reducible diagram (IPR), which is proportional to Cy, contributes to A^, 
however, the irreducible one (IPI) (~ C2{mw) = 1), represented by the charm loop, enters only the 
CP-odd amplitude . In contrast, in neutral pion decays vr'^ 77 the electromagnetic vector coupling 
results in an CP-odd amplitude and hence = and the CP ratio ~ l^+p/l^d" p vanishes. The authors 



of [32], analysed the h — > 577 transition in the lowest order in the full theory (SM). This amounts in the 
calculation of in total 2 x 34 diagrams (the factor 2 is due to the diagrams with interchanged photons). 
They further interpreted the total IPI amplitude as a (local) FCNC 2 photon operator with canonical field 
dimension 8: 

O,^ = {Q,efs [Fl^d^Ff + F^.^d^Ft) IpLh . (4.23) 

After applying the EOM to the photonic part of O^^ the following Lorentz structure is obtained in the Bs 
rest frame 

X = [ei.e2(^i - h) - iiUh - M L , (4.24) 
which leads to the relation Oyy = {QceYki.k2sXb and the IPI amplitude: 



A{Bs ii)ipi{n = mw) = j=-XtC^^{mw){0^^) . (4.25) 



Note that we used on-shell conditions for the photons. The initial value of the "Wilson" coefficient is 
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given as 

Note that O^j is the only operator with dimension < 8 after applying the EOM containing 2 photons and 



2 fermions [33]. One can try to renormalize this dim 8 operator as a point. As a result, the leading order 



anomalous dimension of O^^ vanishes. 



4.2 QCD Sum Rule for the B, 07 Form Factor 

In the description of exclusive B decays hadronic matrix elements < X\Oi\B > are involved. Here X 
is any hadron (with mass mx) and the amplitudes just mentioned are purely non-perturbative objects. 
At present there are two methods to calculate them in a model independent way and, depending on the 
mass value of X, one can be chosen. For so called "heavy-to-heavy" transitions, where X contains one 
heavy quark, the heavy quark expansion technique (HQET) is most appropriate. Heavy quark symmetry 
implies that the form factors, which appear in the Lorentz decomposition of a decay into a hadron with 
certain spin, are all related to one single function, the Isgur-Wise function ^{v.v') [ |124| ]. Here v, v' denote 
the velocities of the B, X, respectively. Moreover, it can be shown that at the zero recoil point, that is 
where the final hadron is at rest in the rest frame of the decaying B, we have the normalization ^{l) = 1. 
The behaviour of ^{v.v') for general values of the argument v.v' / 1 cannot be calculated. It is of non- 
perturbative nature. The advantage is that there is just one universal function describing all transitions in 
the heavy quark limit. From kinematical considerations one can get the possible range of the dot product 
as v.v' = {m\ + m\ — q^) /2mBmx, where = {tubv — mxv'Y- The maximal momentum transfer 
Qmax — {f^B — m-x)'^ Corresponds to the minimal value of v.v' = 1. The decay under consideration 
Bs — > (t)j requires = for an on-shell photon, and we have v.v' = 2.73. For such "heavy-to-light" 
decay an extrapolation too far from the zero recoil point is needed, and the QCD sum rule method is more 
useful. 



The decay Bs </>7 is CKM allowed and, in the language of the operator basis given in eq. (2.18), 



involves the operator Oj, (see eq. ( 4.29| ) below). It has been studied in the literature in the framework of 



Light-cone QCD sum rules [ |125[ ], which is based on the approximate conformal invariance of QCD. Here 
the sum rule is evaluated in terms of meson wave functions on the light-cone. These universal functions 
with increasing twist replace the expansion in "classical" QCD sum rules into many vacuum expectation 
values of operators with increasing dimension. We show the calculation of the form factor Fi in the decay 



Bs —>■ 07 in the framework of the ordinary QCD sum rules [126], including the contribution from the 
gluon condensate [p7|]. 

4.2.1 Calculation of the sum rule 

The amphtude for the Bg — > ^7 transition A{Bs 4>^) =< (l)"f\Heff\Bs > reduces to 



A{Bs ^ = e^Crrn, < (i){p')\sa^^Rq^b\Bs{p) > 



(4.27) 
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with the constant C 



^ = %^^tsVtbCf{li) , (4.28) 
where we just take the contribution due to the electromagnetic penguin operator O7 

O7 = Y^s^a^^imbR + msL)baF^''' , (4.29) 

into account and put = 0, justified by <C mi). Here e and q are the photon polarization and the 
(outgoing) photon momentum, respectively. Lorentz decomposition gives further: 

< (l){p')\sa^,Rq''b\Bs{p) > = ie^upae'^'v'v"' Fi{q^) 

+ {e'lp.q-p^q.e'^)G{q^) , (4.30) 

where p, p' denote the four-momenta of the initial i?s-meson and the outgoing (j), respectively and is 
the polarization vector of the 0-meson. At = both form factors coincide [127] and it is sufficient 



to calculate -Fi(O). Note, that the form factors introduced above are in general functions of two variables 
and p'^. Since the (f) is on-shell, we abbreviate here and in the following unless otherwise stated 

Fi{q^) = Fi{q\p'^ = ml). 



The starting point for the sum rule is the three-point function [128] 



To,^. = - J d'^xe'P^-'P'y < 0\T[J^{x)Tf,{0)My)]\0 > , (4.31) 

where Jq = sjaS, J5 = sij^b and = s^a^i^q'^b correspond to the electromagnetic, pseudoscalar 
currents and the penguin operator, respectively. Performing now an operator product expansion (OPE) 
of Ta^, we obtain a perturbative term, the so-called bare loop, and non-perturbative power corrections. 



diagrammatically shown in Fig. 1.3. The bare loop diagram can be obtained using a double dispersion 
relation in p"^ and p'^, 



' bare 



1 roc rcX) p(g gl\ 

' ds / ds- „ , ,' , — -r- + subtractions . (4.32) 



J ml Jo (s - _ 



Technically, the spectral density p{s, s') can be calculated by using the Cutkosky rule, namely, by replac- 
ing the usual propagator denominator by a delta function: 

i.-2}j^-2 —2m5{k'^ — m?)9{kQ). As a result we get 

8 (s — s'Y 

OPE enables us further to parametrize the non-perturbative effects in terms of vacuum expectation 
values of gauge-invariant operators up to a certain dimension, the so-called condensates. We consider 
up to dimension-5 operators; i.e. the quark condensate, gluon condensate and the quark-gluon (mixed) 



condensate contributions (Fig. |4.3| ). This calculation is carried out in the fixed point gauge, i.e. A^x'' = 0. 
We get 
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Figure 4.3: Contributions of perturbation theory and of vacuum condensates to the Bg 
dashed lines denote soft gluons. 



(pj decay. The 



a. 



' dim— 4 



<G^> 



dx 



1447r JO JO 



1—x nco 

dy I daa^ 
Jo 



im—5 



+ 



1 



rriT 



3(p2 



2(p2 — m'lYp''^ ' 



where 



ci 

C3 



4 4 
mijX 



mlx^{l -x-y) , 
mba;^(3 + y)(l-x-y) , 



mix 



x{l-x-y) , 
y{l-x-y) . 



Here we used the exponential representation for the gluon condensate contribution 

1 _ 1 r 

~~ TyJo 



da a^-^e""^ 



(4.34) 



(4.35) 



(4.36) 



D"" (n 

The momenta P, P' in eq. ( [4.34 ) are Euclidean. 

For the calculation of the physical part of the sum rules we insert a complete set of on-shell states with 



the same quantum numbers as Bs and (j) in eq. ( 4.31 ) and get a double dispersion relation 



'■phys 



mb 



t,m^ 



(p2 _ m|J(p'2 - m^) 



Fi(0) + continuum , 



(4.37) 
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where and fs^ are the leptonic decay constants of the (p and Bs mesons respectively, defined as usual 
by 



< 0\Ja\4> > 
< 0\J5\Bs{p) > 



(4.38) 



We have absorbed all higher order states and resonances in the continuum. 

Now, we equate the hadron- world with the quark-world by Tphys = T^are + + T4 + T5. Using 
quark-hadron duality, we model the continuum contribution by purely perturbative QCD. To be definite, it 
is the part in eq. ( |4.32 ) above the so-called continuum thresholds sq and Sg. To get rid of subtractions and 



to suppress the contribution of higher order states, we apply a double Borel transformation B [126] with 
respect to and p'^. We make use of the following properties of the Borel transform: 



1 



{P 



m 



2\n ■ 



Bier'^P 



(n-1)! [W^Y 



(4.39) 
(4.40) 



Finally, this yields the sum rule 

Fi(0) 



exp( 

m, 

2 

a. 



M2 M'2 fBsfmrnj^^ tt^ J ml Jo 



< ss > e 



m. 



°^4M4 3M2M'2 2M2 



mt 



1 



)] 



+ — < > 

TT JO 



N{x)dx} 



where s = min{s — mf, Sq) and Xmax = m^+m''^ ■ Here we used the parametrization 



g < saGs > 



rriQ < ss > 



(4.41) 



(4.42) 



The last term in eq. ( 4.41 ) is due to the gluon condensate contribution and the function N{x) is defined 
by: 



N{x) ^gexp( M^^l - X -xAiyM'^) 



mlM'^x{mlM'^ - iM'^M'^ + SM^M'^ 



X 



5M^M'^x - M^M'^x^ - 2M^M''^x^ - M^x^) / {M^{-M''^ + M'^x + M'^xf'l^A?>) 



4.2.2 Analysis of the sum rule 



First we list the values of the input parameters entering the sum rules (eq. ( 4.41 )), which are not included 

-0.011 GeV^ |fT30ll, ^ <G^ >= 0.012 GeV^ [ITIel], 



ml 



0.8 GeV^ 1129], < ss > 



in Table |AJ 

= 1.019 GeV and = 0.23 GeV [ |I3T| 1 



We do the calculations for two different continuum threshold values sq = 33 GeV^ and sq = 35 GeV^ 



and take Sq = 1.8 GeV . In Fig. 4.4 we present the dependence of Fi(0) on M and M for sq 



33 GeV^ According to the QCD sum rules method, it is necessary to find a range of K'P and M'^, where 
the dependence of -Fi(O) on these parameters is very weak and, at the same time, the power corrections 
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Figure 4.4: The dependence of the decay constant Fi (0) on the Borel parameters and M'^ for 
So = 33 GeV'^. 



and the continuum contribution remain under control. From Fig. i.4 and Fig. follows that the best 
stabiUty region for Fi(0) is 7 GeV^ < < 9 GeV^, 2 GeV^ < M'^ < 3 GeV^ for sq = 33, 35 GeV^. 
We get: 



Fi(0) = 0.24 ±0.02 



(4.44) 



This agrees for our value of within errors with the result given in the literature, based on Light-cone 



QCD sum rule calculations [125] 



Numerical analysis shows, as also mentioned in [ |128p , that the natural hierarchy of the bare loop, 
the power corrections and continuum contributions does not hold due to the smallness of the integration 
region, and the power corrections exceed the bare loop contribution. The gluon condensate contribution is 
< 1% of the dim-3 + dim-5 condensate contributions and can therefore be safely neglected in numerical 
calculations. 



4.3 The Bg ^ 07 ^ 77 Amplitude using VMD 

Starting with the amplitude for the decay Bg —>■ 07 as input, we calculate the CP-odd and CP-even 
amplitudes in Bg —>■ 77 by using a (/> — > 7 conversion factor supplied by the VMD model. Here an 
extrapolation of the Bs c/ry decay amplitude from p''^ = (needed for Bg 4>^) to p'"^ = 
(required for Bg —>■ 77) is necessary, such that the (p meson propagates as a massless virtual particle 
before converting into a photon. Note that we suppressed in our notation the dependence of the form 



factor = Fi{q^,p' 



/2 



ra. 



on the second argument p . We define here Fi (Q ) = Fi {q = 0,Q ) 
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Figure 4.5: The dependence of the decay constant -Fi(O) on the Borel parameter M"^ for fixed M'^ at 
So = 33 GeV^ (solid) and sq = 35 GeV^ (dashed). 

for virtual momenta = —p'"^. We assume, that the form factor Fi(Q^) is dominated by a single pole, 
which is a good approximation for light mesons and write: 

Fi(0) 



1 



(4.45) 



pole 



Using an nipoie of order 1.7 — 1.9 GeV, which corresponds to the mass of the higher resonances of the 4> 
meson, we estimate -Fi(O) = 0.16 ± 0.02. 

With the help of VMD [^,77,76] and factorization we can now present the amplitude for Bg 77. 



Using the intermediate propagator 
mechanism 



at Q2 



0, the 



7 conversion vertex from the VMD 



< 0|J^ em|</'(p',e) >= eQsf^iO)m^e^ , 
and the A{Bs (pj) amplitude, see eq. ( |4.27| ), we get: 

A(^Bs ^ 07 ^ 77) = e^(A:i)e^(A;2)(A+^^^5(^^ + ^^^£,07 
with the CP-even (Aj^j^^ ) and CP-odd (Aj^^^ ) parts: 



rbi Ail 



A 



LDo^ 



2xCmb- 



m 



Bs 



rrij 



V2^A(0)/,(0)A,!i'''"'« 



TT 



3m^ 



A 



LD, 



■07 



2xCmfeFi(0) 



TT 



' 3mw 



(4.46) 



(4.47) 



(4.48) 



where /^(O) = 0.18 GeV [75], Qs = —1/3 and C is defined in eq. ( ^.28| ). The factor 2 stems from the 
addition of the diagrams with interchanged photons. While for the analysis of the sum rule for Bg (pj 



we have used 



(m?), here we take into account the suppression in going from Q"^ 



m. 



to = 0. We treated the polarization vector e"^ as transversal and replaced e ei, e*^ €2, q 



ki, p' ^ k2. The conversion factor x is defined as x = —eQ, 
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2.2 



1.6 



1.4 



(5.03,370) 

(4.91,460) 

(4.79,590) 




Figure 4.6: Scale dependence of the ratio R{fJ,) defined in eq. (\^.5G ). The solid, short-dashed and long- 
dashed lines correspond to the values {rrih, Ag) in {GeV,MeV) as indicated in the figure. The dotted line 
depicts the suggested choice of the scale fifrom B Xsj studies in NLO [43 4^^^. The parameters 



used are given in Table 4. 1 



4.4 Numerical Estimates of the /i, {rrib, A^) Uncertainties and the O7 Medi- 
ated LD Effects 

We combine in this section the results obtained in the previous sections p~Tl - p3| , i. e., the LLog QCD 
corrections, the HQET inspired bound state model and the Oy-type LD effects in Bg 77 decay. We 
give numbers for the Bs — > 77 branching ratio and CP ratio and discuss the dependences on and the 
uncertainties due to the renormalization scale fi and the bound state parameters {mi,, As). 

Adding the 07-type LD amplitudes (eq. ( [4.48 )) to the short-distance ones (eq. ( [4. 11[ )), we obtain the 
Bs 77 width including the Bs —>■ (j)^ —>■ 77 contribution: 



r{Bs ll)sD+LDo, 



1 



.{A\A+ + Atj,^^]-" + lmUA- + A 



(4.49) 



327rmB^ 

Here a comment about double counting is in order. The "LD" amplitudes considered here, which involve 
the Wilson coefficient C7, also contain a piece from perturbation theory. It originates in the bare loop 
diagram in the calculation of the sum rule and enters the value of the form factor Fi . Without this pertur- 
bative part, it is not possible to perform an operator product expansion; it corresponds to the leading term 
in the sum rule and hence there is no way to avoid it. The "SD" amplitudes on the other hand include also 
contributions from small momenta. As a consequence, by adding the perturbative and the non-perturbative 



parts in eq. (4^.49) there is certainly some double counting present. However, as usual, it is assumed that 
the SD parts are small in regions where the LD ones are large and hence the effect of this is small. 
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First we study the leading logarithmic ^u-dependence of the ratio 

T{Bs ll){mw)sD+LDoj 
In the numerical analysis we neglect the masses of the light quarks, i. e. we use I{mg) 



(4.50) 



1 for q 



u,d,s and msA{ms) = msJ{ms) = in eq. ( 113 ). From Fig. 4^ we find an enhancement factor 
of 1.3 — 2.3 relative to the lowest order result obtained by setting /i = mw, depending on the model 
parameter (mfe, A^). Varying fj, in the range 2.5 GeV < /x < 10.0 GeV, gives an uncertainty AR/R{^ = 
5 GeV) ±(17, 19,22)% for = (590,480,370) MeV, respectively. Here one can argue, that the 
choice /i = ^ takes into account effectively the bulk of the NLO correction as suggested by the NLO 



calculation for B Xg^ [43,44|. 



Table 4. 1 shows the combined /i and model parameter dependence of the branching ratio 

V{Bs -f7)sD+LDo, 



B{Bs 7j)sD+LDo^ 



(4.51) 



The dependence of the form factor Film"^) on the 6-quark mass has been extrapolated from Fig. 3 [125]. 
Here Fi(0) = 0.14, 0.15, 0.16 has been used for ruh = (5.03, 4.91, 4.79) GeV, respectively. Quahtatively, 



(GeV) 


A, = 370 MeV 
rub = 5.03 GeV 


As = 480 MeV 
rrib = 4.91 GeV 


A^ = 590 MeV 
rub = 4.79 GeV 


2.5 
5.0 
10.0 


1.43 • 10"^ 
1.18 • 10"^ 
0.99 • IQ-^ 


8.1 • lO"'^ 

6.8- lO"'^ 

5.9- 10-"^ 


5.0 • 10-^ 
4.3 • lO-'^ 
3.8 • 10-^ 



Table 4.1: Branching ratio B{Bs ^^)sd+lDoj fa^ selected values {nib, ^s) and the renormalization 
scale fi. 

the influence of the LD contribution through Bg ^ (f)^ ^ 77 reduces the width because of the destructive 
interference of the LD + SD contributions. To quantify this, we define 



B{Bs 77)50 



with r{Bs n)sD given in eq. ( 4.22 ). We find, that k lies in the range: 

- 15% < K < -27% , 



(4.52) 



(4.53) 



depending mainly on {mb, A^). To summarize, lowering the scale /i and A^ enhances the branching ratio 

B{Bs ^ 77). 



The dependence of the CP ratio 033], [121], here including our LD Oy-type estimate 
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A* 


As = 370 MeV 


As = 480 MeV 


As = 590 MeV 


(GeV) 


rub = 5.03 GeV 


nib = 4.91 GeV 


nib = 4.79 GeV 


2.5 


79 (0 80) 


88 fO 881 


89 fO 90) 


5.0 


0.69 (0.71) 


0.73 (0.75) 


0.70 (0.73) 


10.0 


0.61(0.63) 


0.60(0.63) 


0.55 (0.60) 


mw 


0.38 (0.41) 


0.33(0.36) 


0.26 (0.33) 



Table 4.2: The CP ratio rcp sd+lDq? given in eq. ( 4.54) for selected values {nib, ^s) <^nd the renormal- 



ization scale /x. The values in parentheses correspond to rcp as defined in eq. I j^.lOj) without taking into 
account the 0-j-type LD effects. 



on the renormalization scale and on the bound state parameters can be inferred from Table The values 
of rc p without taking into account the LD contribution from the decay chain Bg —>■ (j)^ —>■ 77 are also 
shown in parentheses. As can be seen, the Oy-type LD effects reduce the ratio [ ]121| ]. Further including 
the LLog QCD corrections enhance rcp sd+lDq? by a factor of 1.6 — 3.4 compared to the lowest order 
result (fi = ni\Y), depending on the {nib, ^s) parameter set. As a rule, both lowering ^ and increasing 
enlarge the value of the CP ratio. 



4.5 Estimate of the Long-Distance Contribution through h ^ sil) in Bg ^ 
77 Decay 

In this section we estimate the additional LD effect due to the dominant four-quark operators Oi and O2 
(see eq. ( 4.56| )) through the Bg 4>'ip — > (^7 ^ 77 chain decay. We use at quark level 6 — > s-^ followed 
by 6 ^ 57 decay [ f?^ and we pass to the hadronic level using the transition form factor Fi (0) from the 
amplitude A{Bs (/)) [37], [127] given in eq. ( 4.44 ). For both the conversions tp ^ ^ and ^ 7 
we employ the Vector Meson Dominance (VMD) model [^]. The conversion i/; — > 7 needs further 
manipulation because of the strong contribution from the longitudinal part of the Tp meson. We extract the 
transverse part using the Golowich-Pakvasa procedure [^], [77]. Further, we calculate the Oi 2-type LD 
effect to the Bs — > 07 decay using the method given in ref. |[3^], namely, by taking into account the virtual 
c-quark loop instead of the hadronization of the cc pair. This procedure was originally applied to estimate 
the LD effect in i? — > K*^ decay and uses operator product expansion and QCD sum rule techniques. 
Finally we present amplitudes for the decay chain Bs — > iptlj ^ (/>7 — > 77. 



4.5.1 The chain process Bs ^ 4>ip ^ 4"^ - 

We first consider the additional contribution to b 



77 



57 from b — > stpi — > 57, where ipi are all cc J = 1 



bound states, see Fig. 4.7. The relevant part of the effective Hamiltonian describing this process is given 



93 



as 



with the dominant four-Fermi operators 

01 = 

02 = 



(4.55) 



1 



1 



2 

-75 



75-1 
— C/3 



75 



1 -75 



b . 



(4.56) 



Here a, (3 are 5C/(3) colour indices and V^^*^ are the relevant elements of the quark mixing matrix. The 



initial values of the corresponding Wilson coefficients are Ci{my/) = and C2{my/) = 1. To include 
leading logarithmic QCD corrections we evaluate Ci,2(/u) at the relevant scale, jjL nn, for 5-decays, and 
this takes into account short-distance effects from single gluon exchange. The analytical expressions can 
be found in [|6P]. Further we have used the unitarity of the CKM matrix V^, Kb = -V^^^tb - and 



Figure 4.7: The diagram contributing to Bg (pi^ — > (/>7. 
have neglected the contribution due to an internal u-quark, since V*gVub <^ Ws'^tb = 



Using factorization, we obtain the inclusive decay amplitude for the process h ^ sip [ J760 as 

A{b^ siP{ki,e'^')) = -iCf^{ml)m^-s^^'{l-^^)hef, . (4.57) 

Here 



with, assuming naive factorization. 



C = --|Ata2(/i) (4.58) 



a2{^i) = Cl{^i) + ^^ , (4.59) 



where A'^c = 3 in colour SU{'i) and ki,e^ are the momentum and the polarization vector of the ip, 
respectively. In eq. ( 4.57D we used the matrix element 



< 0|c7^c|^(/ci,e'^) >= f^{m^^)m^ef . (4.60) 
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At this stage there is a critical remark about factorization in order, concerning the value of 02 (/i) used. 
The decay under consideration is a class II decay following the classification of [ 132]. In general eq. ( [^.59 ) 
is written as 

af = (Ci(^) + ^) [1 + ei(/.)] + C2{^i)e2{^A , (4.61) 
where ei(/u) and e2(ju) parametrize the non-factorizable contributions to the hadronic matrix elements. 

|-|- 

of takes into account all contributions of the matrix elements in contrast to 02 (/^). which assumes naive 
factorization ei(/i) = e2(/i) = 0. Especially e2(/u), which is the colour octet piece, has sizable contribu- 
tions to naive factorization in class II decays [^9|]. Furthermore, the additional problem is not to know the 
correct factorization scale. In order to include the non-factorizable corrections we use the effective coeffi- 
cient o|^^, which is determined experimentally from the world average branching ratio of B ^ K^*'^'ip as 



[991 



if^ = 0.21 . (4.62) 



This choice restores the correct scale and is fj, independent. Writing the ansatz 



4^ = Climb) + (Ciimb) , (4.63) 



it follows that ^ w 0.41 with Ci {rrib) = —0.25 and C2{nib) = 1.11 for the input values given in Table |A.l 
For comparison, naive factorization would give 02 (m;,) = 0.12. 

Our aim is to replace the '0 meson with the photon 7 and to construct a gauge invariant amplitude. 
We remove the longitudinal component of the meson tj) and then can be converted into the polarization 



vector of the photon 7. We utiUze the Golowich-Pakvasa [|7g],[|77|] procedure making use of the Gordon 
identity, namely 7/^7^ = g^j^a — i(^fj.a- We start with the vertex 57/^(1 — 75)6 and using the equation of 
motion |i5 = mbb and momentum conservation p = p' + ki, we get 

S7/.(1 - I5)b = — {s7^/(l + 75)& + s7m^i(1 + 75)&} , (4.64) 
nib 



where p,p' are the momenta of the h- and s-quark, respectively. We neglect the first term in eq. (164) 
since ^ ^ 1 and p'^ej^ = 0, which follows from e^p^ = in the rest frame of the 5-quark and the 
transversality condition e^k'^ = 0, where is the transversal polarization vector of the if) meson [76]. 
The second term can be written as 

— S7^^i(l + 75)6 = — {s(l + l^)ki^h - iscj^„/cf (1 + 75)^ • (4.65) 
mb nib 



Only the (T^q, term in eq. ( 4.65| ) couples to the transversal component of the and we obtain the corre- 
sponding amplitude as 

A{h ^ sV^) = -2Cf4nil)'^^sa^akf Rbel , (4.66) 

where R = denotes the chiral right projection. Note that the coupling structure is the same as due to 
a direct use of O7 = j^sa ^j^ynibRhF^^'^ with the photon field strength tensor F^^'^ and nis = 0. For 
the ijj^ — > 7 conversion following the VMD mechanism we have 

< 0|J/,,edV'^(A;i,e^) >= eQcU{0)ni^el , (4.67) 
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where Qc = 2/3 and /^(O) is the couphng at kl = 0, see eq. (4.71). Using the intermediate propagator 
of the ip meson at /cf = 0, we get 



^ rrib 



(4.68) 



The expression for the amplitude eq. ( |4.68[ ) can be completed by summing over all cc resonant states 

ip{lS),%b{2S),ip{3770),%b{A040), ^^(4160) and ^(4415) 



^ nib 

The various decay couplings = ) are calculated using 



fl=T{^P,^e+e-) 



3m. 



'ipi 



Q247ra2 ' 



and the measured widths from \ p% and given in Table ^v3| 



Ipi 


/^jGeV] 


U(1S) 


0.405 


f^(2S) 


0.282 


fil:(3770) 


0.099 


/^(4040) 


0.175 


/V'(4160) 


0.180 


/i/)(4415) 


0.145 



(4.69) 



(4.70) 



Table 4.3: Vector meson coupling constants used in the numerical calculations. 
Now we have to extrapolate the couplings /^^(fcf ) from k\ = m^. to kl = 0. We use the suppression 



factor [76] 



^ = /'(15)(0)//|(15)K)=0.12 



(4.71) 



obtained from data on photo production of the ijj and assume k to be universal for the other (higher) 
resonances. We now use eq. ( |4.69| ) to find the matrix element of Bs — > </>7 through the b — > sV'^ 



S7 



transition at quark level. The matrix element [ ]127| ] is given as 

<(t>{p')\-sa^o.Rk'^b\Bs{p)> = ie^upae^''p''p"'Fi{kl) 

+ {e'^^p.ki - P^.kl.e'^)G{kl) , 

and we get the amplitude 

#2 /-Q^ 

A{Bs^^j) = 2Ce>te'^''Y.'-^^eQSe^.upaKp"' 

^ rub 



(4.72) 



+ 9pu 



- ml 



^This is consistent with k = 0.11 |^] based on a dispersion relation calculation. 



(4.73) 
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where ei^, are the polarization vectors and ki, p' are the momenta of the photon and (j) meson, respec- 



tively. We used = 0) = i^i(A;f = 0) | |127 ]. Note, that the form factors introduced above are in 



general functions of two variables kf and p . We abbreviated here Fi{ki) = Fi 



nil) and use 



the value Fi(0) = 0.24 ± 0.02 [37] obtained in section 4.2.2 



VMD vs soft gluon interaction 

Now we want to compare our result for A{Bs 



eq. ( |4.73| ) with the same amplitude calculated by 



the method worked out in [p6|]. This method is based on the new effective quark-gluon operator obtained 
by the interaction of the virtual charm quark loop with soft gluons, in contrast to a phenomenological 
description in terms of t/j resonances converting into a photon, as we used. In this approach, the operator 
Oi does not give any contribution to the matrix element of Bg —>■ 07 for an on-shell photon. The Fierz 



transformation of the operator O2 reads (using eq. (A.13) and eq. (A.21)) 



O2 = 1/N, Oi + 1/2 Ooctet , 



(4.74) 



where 



octet 



4(c7M^r'^c)(.-7^- 



75 



h) 



(4.75) 



and T°- = A"/2 are the SU (3) colour generators. Then the only contribution comes from the colour octet 
part Ooctet- Using the operator Ooctet as a vertex of the virtual charm quark loop, which emits a real 
photon, and taking into account the c-quark-soft gluon interaction, a new effective operator is obtained. 
The matrix element of this operator between Bs and cj) meson states gives the long distance amplitude of 
Bs 07 decay due to the Oi, 2 operators and it is written as (see [Q for details; there the amplitude for 
the decay B ^ K*"fis, given) 



A:[Bs ^ 07) = 2C'4e't"'{ie^,pX^p"' L + 



m 



(4.76) 



where C = 'gm'^'' ■ form factors L and L are calculated using QCD sum rules [36], 



L 



L 



nib 



,nii 



exp( 



m4>ni%jBj<i, 



+ 



M2 M'2- 



,nit 



ds e-' \^ 



^niQ < ss > 



m2/M2 S 

AttOs < SS >2 nib 



M2 + M'2 



nit 



n^_ 

sM"^ 



)(1 + 



M'2^ 
7W' 



12 

nib 



27 



exp(- 



,m 



Bs 



m^rnj^jBsU 
mb .as ^ ^2 . 

i-TTT-o < SS > ml 
^ 12 



+ 



M2 M'2 ' 
ds e~' 



,nir 



+ 



IGvTQs < ss > nib 
27 



M2 + M'2 



mt 



sM"^ 



)(1 + 



M' 



/2 



5M2' 



(i+iik)i-p(-5i)} 



(4.77) 
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The Borel parameters M and M' are varied to find tlie stability region for L and L. We use in the 



evaluation of the sum rules the input parameters given at the beginning of section \.1.2 and Table A.l 
The stability region is reached for eCeV^ < < OGeV^ and 2GeV2 < M'^ < AGdW'^ and we get 

L = (0.30 ± 0.05) GeV^ , 

L = (0.35 ± 0.05) GeV^ . (4.78) 
Writing the amplitude for Bg —>■ 07 as 

A('){Bs - 07) = e>^e'f"'{ie^,p^k1p"^ + g^, A+O) , (4.79) 
and using eq. ( P773| ), ( F!76| ) and ( ^J^ , we can compare the coefficients obtained by the two different 



methods and get 

lA- - A'- 



A~ 
\A+-A'+\ 



< 10% , 

< 5% . (4.80) 



A+ 

This means, that the amplitudes agree within 10%. 

In our approach, the structure of the transition b — > sip'^ 57 is proportional to a^a ^-^^kf (see 



eq. 4^.66), since we removed the longitudinal part of the tp meson from the amplitude. Further, the form 



factors Fi {kf) and G{kj) in eq. ( ^.72[ ) are related for a real photon {kf = 0), Fi (0) = G(0). Therefore, in 



the amplitude A{Bs — > (pj) only one form factor appears, which is Fi(0) in eq. ( [4.73| ). However, the form 



factors L and L in A'{Bs — > ^7) given in eq. (4.76) are not related. They are calculated separately using 



QCD sum rules and this causes the difference between the ratios in eq. ( 180| ). In spite of the fact that the 
amplitudes and A'^ are different from each other, they coincide within the given approximation and 
theoretical uncertainties lying in both methods. 

Oi 2-type LD amplitudes in Bs 77 

We can now present the amplitude for B^ —>■ 77 due to the chain reaction Bs — > (pip — > ^7 — > 77. We use 
the intermediate propagator at zero momentum transfer and the 0^7 conversion vertex from the VMD 
model, 

< 0| eiW, e^) >= eQsUiO)m^ef, , (4.81) 

where the polarization vector is treated as transversal. To apply the VMD mechanism to the amplitude 
eq. ( |4.73 ), we have to know the form factor at Fi{k'l = 0,p'^ = 0). We employ the extrapolated value 



Fi(0) = Fi(0, 0) = 0.16 ± 0.02 [||] from section gj. Then the amplitude can be written with p' k2, 

62 as 

A{Bs -^(P^^cp^^ 77) = e^(/ci)e^(A:2)(5^.^+-^^^ + ie^,,^A:f ^2^^"-^^^) , (4.82) 
with the CP-even and CP-odd A~ parts 



2 



^l-Do = 4X— A(0)^4(0)eQ, 



^l-Do, = 4x^A(0)E4(0)eQc, (4.83) 
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where C is defined in eq. ( [4.58[ ) and tlie conversion factor x is given as x = ~^Qs^^^- Here /^(O) 



0.18 GeV [75] and = —1/3. The extra factor 2 comes from the addition of the diagram with inter- 



changed photons. 

4.6 Final Numbers and Conclusion on the Decay Bs 77 

In conclusion, we have reanalysed the decay rate Bg — > 77 in the SM and we included leading logarithmic 
QCD corrections. Our model to incorporate the bound state effects in the Bg meson is inspired by HQET, 
resulting in the parameters (mfc, A^). The strong parametric dependence of the decay rate T{Bs 77) 
and the CP ratio rcp on (m^, A^) and on the renormalization scale ^ has been studied by us. Further we 
investigated the influence of the LD contributions due to the chain i?s — > ^7 — > 77. Depending on A^j, 
the LD-contributions induced by the operator O7 become sizeable. 

For typical values (mb,As) = (5 GeV, 370 Me V) and /i=5 GeV, we get (including long-distance 
effects through 0-j) the branching ratio B{Bs — > '^^)sd+lDo^ — l-l^ " 10~^, which is a factor 1.9 larger 
compared to the lowest order estimate for the same values of the parameters. However, varying (m;,, A^) 
and i_L in the allowed range results in significant variation on the branching ratio (see Table pTl| ), yielding 

0.38 • 10-6 < B{Bs ^ ll)sD+LDo, < 1-43 • 10"^ . (4.84) 
Improving this requires NLO calculation in the decay rate Bg —>■ 77 and further study of the bound state 



effects. The present best limit on the branching ratio in Bs 77 decay [ |117[ ] given in eq. ( [4. 1[ ) is still a 
factor 100 — 400 away from the estimates given here. 

Likewise the CP ratio rcp sd+lDq? rather uncertain. Varying mh/2 < 11 < 2mb and (mf,, A^) in 
the allowed range, we get in the SM 

0.55 < rcPSD+LDo7 < 0-89 . (4.85) 

Further we presented a VMD model based calculation of the LD contribution to CP-even A'^ and 
CP-odd A~ decay amplitudes for Bg — > 77 decay due to the inclusive process b ^ sip via Bs — > — > 
4'7 ~^ 77 decay. The conversions to photons from both the tpi resonances and the (j) meson lead to two 



suppressions and make the amplitudes in eq. (4.83) smaller than the ones from the LD effect from the 



Bs —>■ (j)^ —>■ 77 chain decay [p7|]. To quantify this we estimated the ratio 



and found 



^l^t^Bs - 07 - 77) \cf{t^)\ r < 



% < /o < 4% , (4.87) 



while varying ^ < /i < 2mb and allowing af^^ to have a theoretical error of 25% as stated in [^]. We 
compared the LD-contribution to Bs — > 77 decay resulting from intermediate tpi production with the one 



obtained by the interaction of the virtual charm loop with soft gluons [36|. We see that both amplitudes 
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are in good agreement within the accuracy of the calculation. The new LD contribution resulting from the 



four-quark operators Oi and O2 is smaller compared to the one of the operator O7 [ 37 ] and affects our old 



estimate given in eq (4-.84) for the branching ratio B{Bs Jj)sd+lDo^ [37] by less than 1%. 

Another LD effect in Bs 77 decay is the one due to intermediate Dg, D* mesons [ 133]. At quark 
level this involves the four-Fermi operator transition b — > ccs. The calculation cannot be done by first 
principles and hence is not straight forward, as the diagrammatic structure is the one of charmed mesons 
in a loop, which are no fundamental particles. The BsDg Dl vertex can be treated with a factorization 
approach, which is an approximation. The next task is to give a prescription of the electromagnetic cou- 
pling of the (charged) Ds mesons, which can be solved by minimal substitution as a first approximation. 
The authors of [133] found a contribution to the branching ratio Bg 77, which is even larger than the 
SM short-distant one. The channel Bg " — > ° 77 surely needs further investigation. 

"New physics"-effects in Bg 77 decay are found to be small as in 6 ^ 577 they are mainly driven 
by the Wilson coefficient Cj^^, for which a strong constraint from data on B ^ Xg'^ decay exists. This 



has accordingly been studied in ref. 0121| ] in the 2HDM and in ref. 0122| ] in the MSSM. However, a small 
enhancement of the branching ratio B{Bg — > 77) compared to the SM one is still possible in some regions 
of the parameter space. In models with an extended operator basis the branching ratio and the CP ratio 



can be much larger than the SM estimates [ ]123| ]. 

Once the necessary machines are running, Bg —>■ 77 will certainly get the same attention as the single 
photon decay b ^ s^y has at present. In particular, the branching ratio B{Bs — > 77) ~ 0(10"^) is large 
enough to be observed at the LHC. 



Chapter 5 



Summary & Future 



Rare B decays are one of the most active fields in recent particle physics. Its main theoretical principles 
and developments of the last 40 years, which are still used, are roughly given as: Of course, the stan- 
dard model (SM) (1961) [|l|] and the quark mixing (CKM) matrix (1963) [9]; further phenomenological 



approaches like vector meson dominance (VMD)(1969) [74] and the Fermi motion model (FM) (1979) 



[22|. The description of low energy weak processes (1981) [134] made progress with the inclusion of 
QCD improved perturbative corrections by using renormalization group equation methods yielding the 
effective Hamiltonian theory (1991) [§-§], [p^, together with the onset of the heavy quark expansion 
(HQE)(1990)[1T5|]. 

We outlined these methods and applied them to the decays B Xgl'^l" with i = e, ^ and (partly) 
to exclusive Bg —>■ 77 decay. We presented quantitative SM based results in terms of distributions, decay 
rates and moments which can be compared with experimental results. 

Concerning B —>■ Xs£^£~ decays, the invariant dilepton mass spectrum and the Forward-Backward 
(FB) asymmetry can be used to extract the short-distance coefficients from data in conjunction with the 
branching ratio in B Xs^ decay. In this work we have analysed these spectra and their present 
uncertainties. 

Further, apart from being a test of the SM, the decay B Xgt^l^ can help to improve our knowledge 
in certain aspects of long-distance effects: 

a) HQE enables a description of 5-meson bound state effects in terms of higher dimensional operator 
matrix elements. We proposed in this thesis the determination of the non-perturbative HQE parameters 
A and Ai from moments of the hadronic invariant mass in B ^ Xst^£~ decays, as it has been done 
for the charged current induced decays B Xlug [|6^. Given enough data, these parameters can be 
extracted from B —^ Xgl'^l^ decays and assuming universality, can be used in the analysis of, e.g., the 
decay B — > Xulvn. Likewise, B Xs£'^£^ decay can be used to test the FM, which can be seen as a 



model dependent resummation of the theory into a so-called shape function [ ]19[ , |135p , and/or to determine 
its parameters. We remark here that some of the HQE and FM parameters are related. 

b) Long-distance (LD) effects occur inB —^ Xs£'^i~ decays via the decay chain B Xs( J/'f, , . . .) 
— > Xs£~^£'' , which we have taken into account with a VMD ansatz. Since in the literature there is no agree- 
ment about the implementation of this LD contribution together with the short-distance one, we compared 
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our approach [|1^] with alternative ones [ |28| , [29| ] and estimated the resulting uncertainties in the observ- 
ables. We find that these uncertainties are not the dominant ones in B Xsi~^£~ decay. Further we have 
shown, that one can reduce the influence of the cc resonances by kinematical cuts. At present, only an 
experimental analysis can identify the correct procedure to implement the charmonium resonances into an 
a parton model based calculation in the decays b —>■ s£~^£~ . 

Finally, by means of building appropriate ratios of (partly) integrated spectra of rare B Xgi^i^ and 
semileptonic B —>■ Xu,c^i^i decays, the uncertainties resulting from bound state effects in the individual 
decays are expected to cancel out to a large extent. 

The essential points reported in this thesis are: 

• The calculation of leading power corrections in spectra and hadronic spectral moments in the decay 



B Xgi^i , including next-to-leading order perturbative 0{as) corrections [[18, 25, 26, T 



The presentation of leading logarithmic QCD corrections to Bg — > 77 decay and an estimate of the 



long-distance effects due to intermediate neutral vector mesons in Bg —>■ 77 decay [p7,38]. 



Besides Bg — > 77, other exclusive decay modes relevant for future ^-experiments are Bg - 
Bg vv and, of course, B {K, K*)t^l^ , B — > {K, K*)vu. The transitions h — > svv are the cleanest 
theoretically among other b ^ s decays. The expected branching ratio is larger (~ 4 • 10^^ [ ]136| ]; note, 
that one has summed over all neutrino flavours) than the one from b s£~^£~. However, the decay 
b — > suu is difficult to observe. Moreover, the inclusive B —>■ XgT^T~ channel is interesting. In the SM, 
its branching ratio is smaller than the one involving light lepton species (e, /x), however, in non-standard 
multi-Higgs models it can be enhanced through large Higgs coupling of the r-lepton. 

Upgrades of present experiments and planned ^-facilities like Hera-B, CLEO, BaBar and Belle are 
about to start soon and they will be sensitive to branching ratios of order 10~^ and below. For an overview 



see Table 5.1 



This work will help the search for flavour changing neutral current B Xgt^i" and Bg — > 77 
decays and in particular, will contribute to precise determinations of the HQET parameters and Vub using 
the inclusive decays B Xg£^£^ and B Xud-f^ in forthcoming B-facilities. 
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CHAPTER 5. SUMMARY & FUTURE 



Expt. 


Collider 


Beams 


1 — 


Year 


£ (1033 


a[bb) 


1, L 

00 pairs 


O 

P^CT 


a[bb) 








(GeV) 


online 




(nh) 


(10^/yr) 


(fim) 




CLEO III 


CESR 


e~^e~ 


10 


1999 


1.2 




1.2 


30 


3 • 10-1 




CESR-IV 




10 


? 


30 




30 


30 


3 • 10-1 


BaBar 


PEP-II 


e+e-t 


10 


1999 


3-10 




3-10 


270 


3 • 10-1 


Belle 


KEK-B 


e+e-t 


10 


1999 


3-10 




3-10 


200 


3 • 10-1 


HERA-B 


HERA 


pN 


40 


1998 




6-12 


50-100 


9000 


1 • 10-*^ 


CDF II 


Tevatron 


pp 


1800 


2000 


0.2-1.0 




20000 


500 


1 • 10-3 


DO 




















BTeV^ 








2004 


0.2 






5000 




LHC-B ^ 


LHC 


pp 


14000 


2005 


0.15 


5 • 10^ 


75000 


7000 


5 • 10-3 


Atlas 
















500 




CMS 





















Table 5.1: Future B experiments. Parameters which do not change between different experiments at the 
same collider are entered only once. ^ asymmetric beam energies, forward detector. ^ 
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Appendix A 

Generalities 



A.l Input Parameters 



Parameter 


Value 


Parameter 


Value 


mw 
mz 

ms 
mc 
mb 
mt 

^(5) 


80.26 GeV 
91.19 GeV 

0.2325 

0.2 GeV 

1.4 GeV 

4.8 GeV 
175 ± 5 GeV 

0.214^°:°^^ GeV 


OLsijnz) 
Ai 

A2 


129 

0.117 ±0.005 
1 

(10.4 ±0.4) % 
-0.20 GeV2 
±0.12 GeV^ 
4.09 • IQ-^^ GeV 
0.2 GeV 
5.369 GeV 



Table A.l: Values of the input parameters used in the numerical calculations, unless otherwise specified. 



A.2 QCD 

The QCD Lagrangian reads in covariant gauge []137|] {A";^: gluon field) 



Cqcd= E q{ip-mq)q-h31,G^^'^ + Csr,+Cghosts, (A.l) 

where the gauge fixing term and the one for ghosts c", c° are given as 

Cfi^ = -^^{d-A'^f, (A.2) 

Cghosts = cfd-D'^'c". (A.3) 
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The chromomagnetic field strength tensor and covariant derivative are written as 

= d,Al-d,Al + gryAlAl, (A.4) 

= d^-igT^Al, (A.5) 

D-^ = 5-'d^-gr'^A-^, (A.6) 

where /"^^ are the structure constants of SU{3), defined by 

[T^ T*] = i/"''^T^ . (A.7) 

We have the identities 

[D^,,D,] = -igT^G^,, (A.8) 

D--Df-Dl-Df = -gr'^G^. (A.9) 

Often the abbreviation is used 

G^, = G%T^ . (A. 10) 



T", a = 1, ... ,8 are the generators of QCD. They are related to the Gell-Mann (3 x 3) matrices A*^ 

2 



through = ^. The T°- obey the following relations k = 1,2, 3) 



Sp{Ta) = 0, (A. 11) 

SpiTaTb) = 6at/2, (A.12) 
"""" -Sij^ki + i:SikSji , (A.13) 



2N, 



TiT,l = 6ikCF, (A.14) 



with the invariant C^? is in an arbitrary SU{Nc) given as 



iV^ _ 1 4 



The coefficients of the QCD beta function (see eq. ( |2.30| )) are written as: 



3 

Here, Nc denotes the number of colours (Nc = 3 for QCD) and Nf denotes the number of active flavours 
(Nf = 5 for the effective Hamiltonian theory relevant for b decays). 



A.3 Feynman Rules 

The covariant derivative consistent with our definition of the operator basis and the corresponding Wilson 
coefficients given in section 2^ is 1 51 ] 

D^ = d^ + igT'^A'l + ieQA^ , (A. 18) 
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where A"" , A denote the polarization four- vectors of the gluon, photon respectively. Note that the sign 



convention of the strong coupling here is opposite to the usual one appearing in QCD text books [137, 



138] given in eq. ([\.5[), but can be made coincident with the substitution g —g. The Feynman rules 



consistent with eq. (A. 18) are given here with boson propagators in Feynman gauge. In a general gauge 
with gauge parameter ^ they are written as: 

. g^, + (e-l)fc^V(F + ze) 
with (^ = 1,0 corresponding to Feynman, Landau gauge, respectively. The rules should be complemented 



I , 1 ^ 



1:9 

■AN\N\N\l\f\r —'i-^^^ab 




by 

• evaluate fermion lines against the momentum flow 

• add a (—1) for a closed fermion loop and perform the trace over the string of 7 matrices 

The rule for an O7 operator insertion is, using 9^ = iq^ for an out going photon and further e • g = for 
a real photon, and F^"' = d^'A'' - d" A^" 

aF = a^.F'^'^ = f 4] = 27^ ^e^ . (A.20) 

The Fierz transformation in d = 4 dimensions is defined as: 

{qi7fMLq2){q3^f,Lq4) = +{qi-f^,Lqi){q^-i^Lq2) , (A.21) 
(^i7^L(i?)g2)fe7M^(^)94) = {-2){qiR{L)qi){q:iL{R)q2) . (A.22) 

A.3.1 Feynman Rules in the Heavy Quark Limit of QCD 

The effective Lagrangian in the limit of an infinitely heavy quark h with mass mg ^ 00 is given by 
C-HQET = hiv.Dh. Here v denotes the velocity (f^ = 1) of a heavy quark h with momentum p = 
rnqv + k and small residual momentum k of order Aqcd- These rules are consistent with the definition 
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h 



J 



k 



2 v-k+ic -i^ 



j k 



of the covariant QCD derivative in eq. (A^), which causes a sign difference in the quark gluon/photon 
coupling compared to the weak effective Hamiltonian rules above, based on the convention eq. ( A.18| ). A 
heavy quark h is represented by a double line. 



A.4 Utilities 



A variety of tools for 1-loop calculations is collected in the appendix of ref. [ ]13 
Distributions 



6ix) 
9{x) 



1 

2^JR 



lim — dq 

e^o 27r JR q- it 



dq 



dB_ 

dx 



(x) = S{x) . 



Geometrical series 



1 

l±x 



n=3 



Fit quality x (see also [39| statistics section), dof: degrees of freedom 



X 
dof 



E 

data points 



, data-curve , 



error 



^data points ^Fit parameters 



(A.23) 
(A.24) 
(A.25) 

(A.26) 



(A.27) 
(A.28) 



^ should be around 1, if it is much smaller, the errors are underestimated, if it is large, the model fails. 

Special Functions useful for loops 
Poly logarithms: 



oo 



k=l 



Li2{z) 



/ -ln(l-t). 
Jo t 



(A.29) 
(A.30) 
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Spence function 

/■I dt 

Sp{z) = Li2{z) = - — In(l-zt), (A.31) 

Jo t 

2 2 

Sp{0) = , Sp{l) = y , Sp{-1) = ^ , (A.32) 

Sp(z) = -Sp(l -z) + ln(z) ln(l - z) , (A.33) 

6 

Sp{z) = -Spi^)-^-^ln\-z). (A.34) 

Useful identities for loops: 

1 1 + iz 1 1 + z 

arctan(z) = — In , arctanh(z) = - In . (A.35) 

Phase space element, S'p= |p|^d|p|(icos ,cos0 G [— 1,+1] ,0 G [0,27r[. 

^= f (fp5{p^ - m^)e{E) ; E = Jp^ + m? . (A.36) 
2E J * 

Dirac algebra identities, for more see ref. []139|], especially the appendix. 



{ipi^lu} = '^Qfiu , (y^y = -[7m>7i^] , l^llu = 9f,u - icr^y . (A.37) 



A useful tool within this context is the TRACER routine [ |140| ] running under the symbolic algebra program 
mathematica. 

Chiral projectors L{R) = (1 =F 75)/2: 

ll = 1, 7^ = 75, (A.38) 

{L{R)f = L{R) , LR = RL = 0, = L{R) . (A.39) 

Further we have 

7o = 1, 70 = 70- (A.40) 

Fermion fields 

= -0^70 = (V'*)'^7o , ipLiR) = L{R)iIj , 

= L{R^={L{R)i;yjo = tpHL{R)yjo = i^hoR{L)=i^R{L) . (A.41) 



Appendix B 



B Xsi^i Dilepton Invariant Mass 



Distributions and FB Asymmetry 



B.l The Functions t/^'^ {v.q, s] 



In this appendix we list the functions Tp\v.q, s), {i = 1,2,3) with j = 0, 1, 2, s, 5, (5, (defined in 



eq. ( 3.55 )), representing the power corrections in 6 ^ si'^l up to and including terms of order m^/ 



ml. 



The parton model contributions T^^^ are given in eqs. ( |B.1| ) - ( |B.3D - 
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(B.l) 
(B.2) 
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(B.6) 



1 rne 
3 mfe'^ 
4 



101 



— 4s — 12^55 + Ssx + QnigSX — 4s v ■ q — ArUgS v ■ q + 7sxv ■ q 



+7m5SX'u • q + 12sv ■ (f + 2Qrh1sv ■ if — Qxv ■ q^ — Gm^xi; • q^ + 4^^ • ^ 

+Amlv ■ - 4xv ■ f - Arhlxv ■ f - 8v ■ q^ - 8ih1v ■ q^ ICf^l^ 

4 



+- 



SX"^ 



4s^ — 5s,T — Asv ■ q + Aifilsv ■ q + 3xv ■ q — Sfhlxv ■ q — Asv ■ (f + 2xv ■ 



+Av ■ q^ - Amjv ■ q' 



Re 



2 niB 

3 nib^ 

\cf 

1 ruB 
3 rrift^ 
4 



Ai 



eff|2' 



V • g 



(B.7) 



(B.8) 



Ai 



{s-{v q)') 



x^ 



+- 



s^x^ 



-As^ + hsx + 8sv ■ q — Qxv ■ q + Asv ■ q^ — Axv ■ q^ — Sv ■ q"^ (1 — m^) |Cy 



^eff|2 



+- 



sx-' 



{As -3x-Av + ml) - 2xv 



Re 



2 ms 
sx rrib^ 

0, 

2 me 



(A1 + 3A2) {s-vq)Re (Cf^^ ^ Cio)* Cf" +2m^|(7^ 



-.effl 



,2|--ieff|2 



s X nif, 
1 ms 



3 (Ai + 3A2)Re \icf ^ Cio)* cf 



(B.9) 

(B.IO) 
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(B.18) 



Here the variable x is defined as x = I + s — 2 {v ■ q) — ml + ie. 



B.2 Auxiliary Functions Ei{s, u) and £"2(5, u) 



In this appendix we give the auxiliary functions Ei{s,u) and E2{s,u), multiplying the delta-function 
6[u{s, rhg) — -u?] and its first derivative 6'{u{s, rhg) — u'^], respectively, appearing in the power corrected 



Dalitz distribution (PB/dsdu{h — > sl'^l ) in the HQE approach given in eq. ( 3.60 ) 



1 - + Qnvt - 4m^ + - 2m^J + 4m!s - 2rh\s + 2m^s^ - 
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E2{s,u) = -Xi(l- TJil + sj u{s,ms)'^ 



l-2mi + mt-s'-u'][ \C^''\' + \a 
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B.3 Dalitz Distribution d'^r{B X^tr) /dsdu and FB Asymmetry in the 
Fermi Motion Model 




Figure B.l: Phase space boundaries for the u' and p integrations with fixed values of s and u drawn for 
s = 15 GeV^ and u = 8.9 GeV^. The integration region (solid curve) is given by the intersection ofu'_^ 
(short dashed) and ±u(^s,p) (long dashed curve). The Fermi motion parameters used are {pF^^q) = 
(450, 0) in MeV. 

We start with the differential decay rate d^rs/ds du dp, describing the decay b — > s£~^£' of a moving 
6-quark having a mass mi, = nn^ip) and three momentum \p\ = p with a distribution (^{p), which will be 
taken as a Gaussian [22], 

dsdudp Ju'^,^ ^ ^u'^ + 4mb^s L^*^^' 

Here, d^r^/ds dn' is the double differential decay rate of a b-quark at rest and can be written in the case 
of 6 ^ si+e- as. 



(B.l) 



d^n 

ds du' 



= \^t*s ^^»\' ^if-s ^ Tr-2 l^^^'^P) + u' + F,is,p) u"] , (B.2) 

192 VT'^ rrife'^ lovr^ L J 



and the three functions have the following expressions. 
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F2(s,p) = 4sRe{cfCio)+8{mb' + ms')CioC^^\ 



^eff 



F3is,p) 



"^b lr'eff|2 
— 1^7 I ' 



(B.4) 
(B.5) 



which can be read off directly from eq. ( 3.60| ) in the limit A, = 0; i = 1,2. Note that the Wilson 
coefficient Cg^^^ also has an implicit m,;, dependence, as can be seen in the text. The integration limit for 
u' is determined through the equations 
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A typical situation in the phase space is displayed in Fig. |B.1| . Integration over p gives the double 
differential decay rate including the Fermi motion. The result is, 
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(B.ll) 



Note that the upper limit in p integration, Pmax is determined such that p satisfies, 

^max(Pmax) -S, "u) = li j (pmax ) ■S, m) . 

Lastly, the normalized differential FB asymmetry including the Fermi motion becomes, 

_ J-Upi, dsdu Jo AsAu 



where 



and 



mx = Max [niK , Jti^ + mq] , 



(B.12) 

(B.13) 

(B.14) 
(B.15) 
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with rriq the spectator quark mass and niK the kaon mass. Since the calculations are being done for an 
inclusive decay B X^i^i^ , we should have put this threshold higher, say starting from tuk + ^1^^, 
but as this effects the very end of a steeply falling dilepton mass spectrum, we have kept the threshold in 
B Xs£+e- at m{Xs) = m^. 



Appendix C 

B Xsi^i~ Hadron Spectra and 
Moments 



C.1 Coefficient Functions '^^^ , gP , '^^ , , h!^^^ , /cf ^ , k^'^ 

These functions enter in the derivation of the leading (1/m^) corrections to the hadron energy spectrum 
in S ^ Xsi^l^ , given in eq. ( |3^ . 



9o'^^^ = - m2--(-2m^ + 3x0 + Sm^xo - 4xo) , (C.l) 

^(9,10) ^ ^ -^(9"^^ + ^3?Tif - ISmgxQ - 18x1 - 52m^xg + 36xg + 203;^) , (C.2) 
Xq — rhg 

^(9,10) ^ ; ^ ^(3mg + 23mf-3xo-21mgxo-6xg-52mgxg + 36x|] + 20x^) , (C.3) 
Xq — 

gj^^ = Jxq — m^ — (lOm^ + lOm^ — 3xo — ISm^xo — 3m^xo + 2xo + 2m^Xo) , (C.4) 

= ; / 7 1-^-— ^^(97n2+34m^ + 104m6 + 110mf + 31mi° 



- 132?7ifxo - 312mf xo - ISOmfxo - 18xg - ITOm^xg - 58mfxg + 74rfif xg - 20mfxg 

+ 72x11 + fyUmlxl + 576m^xg + 228m^x|j - 116x^ - 676771^3;^ - 436m^x^ - 20m^x^ 

+ 72x^ + 2Amlxl + 24mfx^) , (C.5) 

02'^ = — 1^=^= Y', ^ — (27m^ + 93m^ + 97rfig + 31mf — 3xo — 63mgXo 



- 189mf xo - 129mf xo - 18xo - 108m^Xo - 62m^Xo - 20mf Xq + 72xo + 324^^X0 

+ 180m^x|] - mxl - lb2mlxl - 20mfx^) , (C.6) 



5o = Vx2-m2i28(-2m2 + xo + m2xo) , (C.7) 

^(7,9) ^ ^ ^ Ujrhl + 3m^ + 2m2xo - 2xg - ^mlxl) , (C.8) 
Xq — 
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C.2 Auxiliary Functions /^(Ai, A2), fs'i^i: ^2) 

The auxiliary functions given below are the coefficients of the singular terms in the derivation of the 
leading (l/m2) corrections to the hadron energy spectrum in B ^ X. 

"2 



, given in eq. ( [3.95[ ) 
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/5,(Ai,A2) = BoXi\lil-mlf{\cf\^ + \a 
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2^4 



(-1 + 7n1)^{-l + 14mf + rrig + 52mfm1 + + lAmfm 
+ ^(l-m^)5i?e(Cfff)C,^ff}. 



eff 1 2 



(C.17) 



C.3 The Functions , /3j , 7^ , (5^ 



The functions entering in the definition of the hadron moments in eq. (3.99) are given in this appendix. 
Note that the functions aj"''"^ and /J^"'"^^ multiply the Wilson coefficients jCy^^p and C^q, respec- 
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lively. The functions ^ (J^"'™) result from the Wilson coefficients C^^^ Re{Cg^^) , |Cg^^p, respec- 

tively. They cannot be given in a closed form since Cg^^ is an implicit function of a;o. 



The functions a^"'"*^ 



^(0,0) ^ 16 _ 2g^2 ^ ^^^4 ^ 22ml - H^?) + ?(-! + ^tfi^ + ^t) ln(4mf ) 
9 o 
64 



ym^(-9-2m2 + m^)ln(m,) , (C.18) 
(0.0) = l^o.o) , (C.19) 
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4°'°) = ^(-4 + 38m2 _ 42mi - 26m^ - 15mf ) + 16(-1 + 7fi1f{3 + 8m^ + 5mj) ln(4mf ) 

+ 32m^(-8 - 17m^ - 2m^ + 5m^) In(ms) , (C.20) 
a(°'^) = -(-41 - 49m^ + 256?fi^ - 128rfi^ - 43mf ) + — (-1 + m^)3(l + rhlf \n{Amf) 

+ ym^(3-m2-2m^)ln(m,) , (C.21) 
af''^ = , (C.22) 

4°'^) = ^(21 + 167m^ + 128m^ - 276m^ - 319m^) + y (-1 + mlf{3 + Uml + 21m^ 

32 

+ lOmf ) ln(4?fif ) + — m^(3 - 24mf - 18rfif + m^) ln(m,) , (C.23) 

O 

a(°'^) = —(-119 - 144m^ + 45m^ + 320m^ + 45m^) + -(-1 + m^f ln(4mf ) 

- ym6(8 + 3rfi2)ln(m,) , (C.24) 
aj°'^) = -^(-127 - 278m^ + 1075rrif - 800m^ + 49mf) + ^(1 - rn^f{-^ - '^^rh] 

- 5m^ + bml) ln(4?fif ) + ^m^(18 - 38m^ - 13m^) In(ms) , (C.25) 
4°'^) = ^(27 - 46m^ + 1681m^ - 688mf - 1189mf) + ^(-1 + m^)^(3 + 20m^ + 25m^) ln(4m;^) 

- ^m^(18 + 54m2 + 47m^) In(ms) , (C.26) 
4^'°^ = , (C.27) 

= -(-23 - 159mf - 112rnf + 304m^ - 45mf) (-1 + rn^t{l + mf) ln(4mf ) 

9 3 

+ — m^(-39 - 7ml + 6^s) In(ms) , (C.28) 



2 16 
4^'°) = -(-93 - Amml + 704mf - 127mf) + —{-! + mlf{?> + 8m^ + 5m^) In(4mf ) 

112 



.,(1.1) 



m^(3 + 3m^ + 2m^) ln(m5) , (C.29) 



, (C.30) 

2 8 
a(^'^) = —(-4 - 131?fi^ + 307mf - 416m^ + 178mf ) - -(-1 + mlf{l + 6m^ + 5rfi^) ln(4?fif) 

16 

+ -^"is(9 - 35m^ - 7m^) ln(m5) , (C.31) 



a(^'^) = ^(-60 - 185?fi^ + 173m^ + 160?fi^ + 70?fi^) + |(-1 + rhlf{l + rhlf{2, + bml) ln(4rnf ) 
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1 a 

+ — m^(3 - 21ml - l^m^) H^s) , (C.32) 
o 

4^'°^ = , (C.33) 
aJ^'O) = —(119 - 176ml - 1085?fi^ + 400mf + 835?fi^) + -^-(l - rhlf{l + ml) ln{Amf) 

- ^m^(28 + 5m2)ln(m,) , (C.34) 

4^'°^ = . (C.35) 



The functions /3, 



,(n,m) 



/?(°'°) = _(i _ Sm^ + 8m^ - - 24m^ ln(m,)) , (C.36) 

= l^r^ , (C.37) 
4°'°) = -3 + 8ml - 24mf + 24mf - 5rfif - 24m^ ln(m,) , (C.38) 
= -L(7 -25m^ +160m^ - 160mf + 25m^ -7mJ° +120m^ln(m,) + 120m^ ln(m,)) , (C.39) 
= /?f , (C.40) 
= ^ml{7 - 20ml + 20^^ - + 241n(ms) - 48m^ In(ms)) , (C.41) 

= A(2 _ Srn^ - 30m^ + 30mf + 3mJ° - 2ml'^ - 120m^ In(rfis)) , (C.42) 

= ;r^(43 - 135m^ + 1260m^ - 1440ra^ + 405m^ - 153m^° + 20m^^ + 1080m^ 111(7713) 
+ 840771^ ln(77is)) , (C.43) 
= ^(13 - 315771^ + 1500771^ - 1560771^ + 315771^ + 14777i^° - 100771^^ + 360771^ ln(7ns) 
+ 840771^ ln(77is)) , (C.44) 
= , (C.45) 
/?J^'°) = _L(i3 -135^2 _i607nf+3207nf -4577if + 7?7if -600777^ In (7775) - 120777^ In (771^)) , (C.46) 

= 1(3 _ 9^2 ^ 16777^ - 48771^ + 457nf - 77nf + 247n^ ln(7ns) - 727n^ In(7ns)) , (C.47) 
= , (C.48) 

P[^''^^ = ^(23 - 45777^ + 1080777^ - 1440777^ + 585777^ - 243777^° + 40777^2 + 1080777^ ln(?773) 

+ 600777^ ln(777s)) , (C.49) 

/^(i'^) = ^(13 + 45777^ - 120777^ - 45777^ + 147777^° - 40777^^ + SGOmj ln{rhs) 

- 600777^ ln(777s)) , (C.50) 

= , (C.51) 
/3f = _(_i + grhl - 45777^ + 457n^ - 9777^° + 777^^ - 120777^ ln(7hs)) , (C.52) 
/?f '"^ = . (C.53) 

The functions 7^"'"*^ 

Note that in the expressions given below Cf ^ = Cf ^(s = 1 — 2xo + ml). The lower and upper limits 
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of the XQ-integrals are: Xg^*"' = TUg and Xg"*^ = ^{1 + mj — Arhf). 



maa; 

7^°'°^ = 128 / ° dxoJxl - rhl{-2ml + xq + rhlxQ)Re{cf^) , (C.54) 
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7i°'^'' — ~^ I dxoJxQ — m'j{—4:rhg — lOm^xo + 2mgXo + 6xq + IGm^XQ — 5xq 

- 5m2a;3)i?e(Cff^) , (C.61) 
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7^^'^^ = — — / dxo — (4to^ + ArhgXo — 2m^xo — 3to^Xo — ITrn^Xg + wi^Xq + 7m^Xo 

3 . _ " 2 



'Xq - 

+ 3x^ + 13m^x^ - 5x1 - 5m^x^)i?e(C^ff) , (C.67) 
128 pxq^^"^ I 

7^^'^^ = — I dxQXoJxl - m2(-4m^-6m^+3xo-15m^xo +7x^ +15m^x^)i2e(Cf^^) ,(C.68) 
3 ixg^'" 

7^''°^ = , (C.69) 
7^ = / ° dxox/xg - m2(-2m^ + m^xo + m^xo + 2m2x^ - xg - 7nlxl)Re{cf^) , (C.70) 

7^^'°^ = . (C.71) 
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The functions 6^ 



(n,m) 



^(0,0) ^ 32 r^o ^1^2 _ rnl{-2ml + Sxq + Sm^^o - Axl)\Ct^\^ , (C.72) 

^f'^^ = ^4°'°^ , (C.73) 

^(0,0) _ / dxp (6r?7,g — xq — QrhgXo — 12rn^xo + Sm^Xg + 15xq + ISm^Xg 

•^-S-*" ^xl - ml 

- 2{)xt)\cf\\ (C.74) 

4°'^^ = ^ / dxoxojxl - ml{-2ml + Sxq + Srfi^xo - 4a;^)|(7fffp , (C.75) 

,5f = , (C.76) 

<^2°'^^ = — / dxp — (6mg — 3m^xo + 13m^a;o — 3xo — 30m3a;o — 57m^Xo + 3a;o 



+ 43m^a;i] + 48x^ + eOm^x^ - mx\)\Cf^f , (C.77) 
4°'^^ = ^ I dxQxlJxl - m1{-2ml + 3x0 + 3rfi^xo - 4xl)\cf^f , (C.78) 



,5}°'^) = / ° dxoJxl - m2(-4m^ - Gtti^xo + 6m^xo + 18xg + 20m^xg - 39xg - 15m^x^ 

+ 20x^)|CffV, (C.79) 

(^2^'^^ = — / (ixp — ^ (12mg — 6m^xo + Stti^xq — 3xq — 33m^Xo — 78m^XQ + 6x9 
"'^o ^^xf^-mi 

+ eSm^xg + 51x^ + Ihmlxl - 100a;^)|Cf^|^ , (C.80) 
4^'°^ = , (C.81) 

(1,0) 32 ^.4,0-^2 ,0-4 o-~22 q3 q-23 

0^ = -;r f dxp — (— zm^ + orngXp + onigXp — 2m gX^ — 6Xq — 6m gX^ 

3 ^Jxl-rhl 
+ 4x^)|Cfff|2, (C.82) 
'^^^'^^ - ' ° dxo\/x§ - m2(-6?7i^ + 3x0 + 5m^xo + 3xg + 15m^x^ - 20xg)|Cf^^|2 , (C.83) 

in. • 



4^'^^ = , (C.84) 
(5f '^^ = — I dxQ , =(4m^ - 6m^xo - 9m^xg - 15m^xg + 27m^x^ + 21m^xg 

+ 9xo - 9m^Xo - 27x^ - 15m^x^ + 20x^)|C|^p , (C.85) 

4^'^^ = — I dxQxJxl - m?g{-Qml + 3xo + 5m^xo + 3xg + 15m^xg -20xg)|Cf^^|2 , (C.86) 
3 Jxg**" 

4^'°^ = , (C.87) 

128 r^o''"'^ I 

ij(2>o) _ / dxo\ Xq — ml{—2m^g + 3777.5X0 + 3771^X0 — 2mlx^ — 3xo 

- ^rhixl + Axt)\cf\^ , (C.88) 

(5^ = . (C.89) 
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{xo) 



B_ 



+ 
+ 



+ 

+ 
+ 



9m| 



J 

L 



ims/niB 
■i(l+m2/m|) 

Ima/mB 



m. 



16 
3m^ 



TT 



:(l + 31n(4q-))C?" +^ 



rrii 



dxoi6^xlReicf)cf + ^(3 - 4xo)xg|cf ^p) 



m. 



(19 + 12 ln(4^))Cfff^ + [' dxo{^xo{-S - Qxg + 28xl)Re{cf^)Cl 
ttitd Jo o 



^eff 

'7 



xo{-3 + 3a;o + 48x^ - 80xl)\cf^\'^) 



m 



2 ' 



(C.90) 



45m^2 



:-119mi - 144m^V2 _ go(^^2 _ ^12) in(4^))ceff2^_l^(2m| - Sm^Cl 



+ 



+ 



+ 



+ 



i(l+m?/m|) 



64 



.effx^eff 



o?a;o£^aJo(-m^ - 4m^xo + 2m|a;o + 2m^Xo)-Re(C9")C7' 



ml 



i(l+m2/m|) 

dxo Q^2 ^2 «^2^2 _ e^2 ^3M/^eiI|2 



■a;o(— 3m^ + 6m^a;o + 6m^Xo — 8m^Xo)|Cf 



TT 



— + -irC'i 
niB o 



+ 



^(55 + 841„(4^,)C?ff%43|s 



(ia;o(^(6 - 5.xo)xi]i?e(Cfff)C7^ff + ^(18 " 39:^0 + 20xl)xl\cff) 
o y 



mi 



(11 + 41n(4^))Cf"' + 13^ + / ' (ixo(^x2(-3 - 6xo + 35xl)Re{Cf^)C. 



r2 



ml 



90 



64 



^effxrveff 



7 



16 



+ ^xl{-S + 6x0 + 51x1 -100xl)\Cl 



eff|2^ 



m 



2 ' 



(C.91) 



Do TT 



^(l + 31n(4f|-))C?" +23^ 



27 



ml 



r2 

270 



+ dxo{^{3 - 5xo)xlReicf)cf + ^(9 - 27xo + 20xl)xl\a 
Jo o 9 



eff|2-j 



+ 



-(5 + 31n(4^))C7eff' + 13.^ 



10 



m. 



90 



+ r dxoi^xliS + 7xo)Re{cf)cf + ^xg(3 + 3xo - 20^2) |C; 

Jo o 6 



m| 
(C.92) 

Ao 



nil 



£>0 



TT 
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+ 



+ 



+ 



dxo(128(l - xo)xlRe{cf^)cf^ + ^(3 - 7xo + Axl)xl\cf^f) ^ 



mj 



-(31 + 241n(4^))C7^ff + 



•-^10 



rrir 



(C.93) 



.1 

/ dxo{ 
Jo 



128 



32 



^ xli3 + 7xo)Re{cf)cf + ya;^(3 + 3a;o - 20x'o)\C^ 



eff|2^ 



mi 



B 



a. 



TT 



^(2,0)^2 + 



135 



(119 + 601n(4 



ml 



^^135 



(C.94) 



+ 



dx^{^xtRe{cf)cf + ^(-3 + ^xo)xl\C^^ 



effi2N 
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